\ 



■ ■ COMPUTER SCIENCES CORPORATION 



CSC/TR-75/6001 


PRELIMINARY DRAFT 

MULTISATELLITE ATTITUDE DETERMINATION/ 
OPTICAL ASPECT BIAS DETERMINATION (MSAD/OABIAS) 
• SYSTEM DESCRIPTION AND OPERATING GUIDE 

VOLUME 3 

, . _ -OPERATING GUIDE 

Prepared for 

GODDARD SPACE FLIGHT CENTER 
By 

COMPUTER SCIENCES CORPORATION 
Under 

Contract NAS 5-11999 

Task Assignment 396 - 

— ■ ' « 

Prepared by: Approved by: 



M. Joseph 

Date 

R. D. Cardwell 

Date 


J. Keat 

K. S. Liu 
M. E. Plett 


Technical Area Manager 



i>VL# oXjl(»«>LX 

T. Shinohara 
J. R. Wertz 


F. J. Douglas 

Quality Assurance Assistant 

Date 




R. L. Taylor 
Project Manager 

Date 



PRELIMINARY DRAFT 




Page intentionally left blank 


PRELIMINARY DRAFT 


PREFACE 


This volume is the third of four which make up the Multisatellite Attitude 
Determination/Optical Aspect Bias Determination (MSAD/OABIAS) System 
Description and Operating Guide. The volumes are 

Volume 1 - Introduction and Analysis 
Volume 2 - System Description 

Volume 3 - Operating Guide 

Volume 4 - Program Listing and Sample Execution 

This volume contains a complete description of all MSAD/OABIAS NAMELIST 
control parameters, a description and a sample of all printed output unique 
to OABIAS and of each IBM 2250 graphics display, an explanation of and user 
response for all error messages generated by the MSAD/OABIAS System, and 
a listing of the job control language (JCL) required to operate the system. 

Volume 1 contains an introductory exposition of the MSAD/OABIAS System 
and describes the analsrtic basis for the OABIAS aibsystem. This includes a 
detailed discussion of the recursive estimator algorithm, each of the 12 state 

I 

vector elements, and the 8 observation models used. 

Volume 2 describes the system flow and the components of the MSAD/OABIAS 
System, The table language descriptions provide detailed information relating 
the operational displays on the IBM 2250 display device to specific COMMON 
areas and subroutines within the MSAD/OABIAS System. 

Volume 4 contains the program listing with supplementary output and line- 
printer plots of all IBM 2250 displays occurring during a sample execution of 
the program. This volume preserves, in source form, the MSAD/OABIAS 
System as it is presented in this document. 
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ABSTRACT 


This document describes the Multisatellite Attitude Determination/Optical 
Aspect Bias Determination (MSAD/OABIAS) System, designed to determine 
spin axis orientation and biases in the alignment or performance of optical 
or infrared horizon sensors and Sun sensors used for spacecraft attitude de- 
termination. MSAD/OABIAS uses any combination of eight observation models 
to process data from a single onboard horizon sensor and Sun sensor to deter- 
mine simultaneously the two components of the attitude of the spacecraft, the 
initial phase of the Sun sensor, the spin rate, seven sensor biases, and the 
orbital in-track error associated with the spacecraft ephemeris information 
supplied to the system. In addition, the MSAD/OABIAS system provides a 
data simulator for system and performance testing, an independent determi- 
nistic attitude system for preprocessing and independent testing of biases de- 
termined, and a multipurpose data prediction and comparison system. 

MSAD/OABIAS has extensive capabilities for an interactive graphics mode and , 

makes use of the Graphics Executive Support ^stem (GESS), formerly known 

as the Multisatellite Attitude Determination System (MSAD), services. MS AD/ 

OABIAS is a multisatellite system capable of supporting, in its present form, 

3 

Small Scientific Satellite (S ), Interplanetary Monitoring Platform (IMP), At- 
mosphere Explorer (AE), and Synchronous Meteorolcgical Satellite (SMS) 
missions or any similar missions using optical or infrared horizon scanners 
and providing attitude data that can be read by the MSAD/OABIAS System. 


V 

PRELIMINARY DRAFT 



Page intentionally left blank 


TABLE OF CONTENTS 


VOLUTvIE 1— INTRODUCTTON AND ANALYSIS 
Section 1 - Introduction 


1.1 

1.2 

1.3 


1.-3.1 

1.3.2 

1.3.3 

1.3.4 

1.3.5 


1.4 

1.4.1 

1.4.2 


Motivation for Bias Detei’mination . . . . 

Historical Bacl'Cground 

System Overview . . . . 

OADRIV 

ODAP 

OASYS . 

OABIAS 

PLOTOC 

OABIAS Recursive Least-Squares Filter 

Introduction 

Observation Models 


Section 2 - Sensor Descriptions 

2.1 

2.2 

2.3 

2.4 

Section 3 - Analvsis 


Sun Sensor 

Horizon Detector 

AE-C Wheel-Mounted Horizon Sensors 
RAE-B Panoramic Attitude Sensor . . 


3. 1 Introduction 

3.2 . OASYS — Deterministic Attitude Determination 

Subsystem 

3. 2. 1 Input to Deterministic Process 

3.2.2 Design Assumptions 

3. 2. 3 Deterministic Logic 

3. 3 The Recursive Estimator Approach 

3.3.1 Comparison Between Recursive Processing and 

Batch Processing Estimation Methods 

3.3.2 The Basic Recursive Estimator Algorithm . . . 

3. 3. 3 Discussion of Weighting Factors 

3.4 Implementation of the Recursive Estimator 

Algorithm in OABI/iS 

3. 4. 1 Principal Inputs 

3. 4. 2 Observation Models 

3.4.3 State Vector Elements 


1-1 

1-1 

1-2 

1-4 

1-6 

1-7 

1-7 

1-8 

1-9 

1-9 

1 - 9 
1-11 

2 - 1 
2~i 

2 r^ 
“-0 

2-5 

2-6 

3-1 

3-1 

3-4 

3-4 

3-5 

3-6 

3-16 

3-16 

3-19 

3-25 


3-28 

3-28 

3-29 

3-30 



PRELiMlNARY DRAFT 

TABLE OF CONTENTS (Cont'd) 

Section 3 (Cont'd) 

3.4.4 Iterative Operation 3-32 

3.4.5 Single Observation Iteration 3-34 

3.5 Basic Geometr}’- 3-36 

3.5.1 Coordinate Frame Flow Diagram 3-36 

3.5.2 Coordinate Frames GI* and GI 3-38 

3.5.3 Coordinate Frames SI and SC and State Vector 

Elements Xg (ip) and Xg (cc) 3-40 

3.5.4 Sxm Sensor Geometry and State Vector 

Elements Xg (A/3) and ^10 3-43 

3. 5. 5 Horizon Detector Geometrv and State Vector 

Elements X 4 (Ay) , Xg ( 0 ^ , Xg (0§) , 

and Xj^j^ (Cjj) ............................. 3—43 

3. 5. 6 State Vector Element x <7 (Ap) 3-48 

3 . 5.7 State Vector Element Xj 2 (At) 3-50 

3. 5. 8 Transformation Matrix A and State Vector 

Elements (s^) and X 2 (S 2 ) 3-50 

3.5.9 Transformation Matrices Bj and Bq - 3-53 

3. 6 Observation Models 3-54 

3.6.1 Model 1 — Sun Angle Model 3-55 

3. 6.2 Model 2 — Sun Sighting Time Model 3-58 

3. 6 . 3 Model 3 — Nadir Vector Projection Model 3-61 

3. 6.4 Model 4 — Horizon Crossing Time Model 3-66 

3. 6 . 5 Model 5 — Sun to Earth- In and Sun to Earth-Out 

Dihedral Angle Model 3-70 

3.6.6 , Model 6 — Earth Width Model 3-74 

3.6.7 Model 7~Small Target Model 3-77 

3. 6 . 8 Model 8 — Sun to Earth Mid- Scan Dihedral Angle Model • • . 3-79 

3. 7 Computation of Central Body Angular Radius p^, 3-83 

3. 8 Computation of Horizon Crossing Vector D 3-87 

3. 9 Weighting Factors of the Observation Models 3-92 

3.10 Summary of Section 3 3-97 

Appendix A - Model Observation Partial Derivatives 

Appendix B - Derivation of Recursive Processing Algorithm 
Used in OABIAS 

Glossary 

viii 

PRELIMINARY DRAFT 



PRELIMiNARY DRAFT 


TABLE OF CONTENTS (Cont'd) 

VOLUME 2— SYSTEM DESCRIPTION 

Section 4 - System Flow 4-1 


4.1 Overview 4-1 

4. 2 External Interfaces 4-4 

4. 3 System Flow, Nongraphic Mode 4-5 

4.4 System Flow, Graphic Mode 4-6 

Section 5 - Baseline Diagram and Subroutine Descriptions 5-1 


5. 1 Baseline Diagram 5-1 

5.2 Subroutine Descriptions 5-1 

Section 6 - COMMON Area Descriptions 6-1 


Appendix C - Description of GESS Tables 

VOLUME 3— OPERATING GUIDE 

Section 7 - Operating Guide 7-1 


7.1 Introduction 7-1 

7.2 Resources 7-1 

7.3 Card Input 7-1 

7.3.1 . NAMELIST CONTRL . .' 7-4 

7. 3. 2 GESS Array Allocation Sizes 7-5 

7.3. 3 GESS Display Status Flags 7-5 

7.3.4 NAMELIST MAIN 7-6 

7.3.5 NAMELIST LIST 7-19 

7. 3. 6 NAMELIST OPMANl 7-31 

7.3.7 NAMELIST BIASNL 7-40 

7.3.8 Description of a Sample NAMELIST 7-49 

7.4 OA BIAS Data Set 7-60 

7.5 AE Data Set 7-63 

7. 6 JPL Lunar and Solar Ephemeris File 7-64 

7. 7 SUNRD Lunar and Solar Ephemeris File 7-64 

7. 8 GTDS Spacecraft Ephemeris File 7-65 

7. 9 DODS Spacecraft Ephemeris File 7-65 

7.10 Spacecraft Attitude Tape 7-66 

ix 

PRELIMINARY DRAFT 


TABLE OF CONTENTS (Cont'd) 


Section 7 (Cont'd) 

7.11 GESS Nonresident Tables 7-66 

7.12 Printed Output 7-66 

7. 12. 1 Printed Output From the GESS Executive 7-66 

7. 12. 2 Printed Output From OADRR^ 7-68 

7. 12. 3 Printed Output From the ODAP Subsystem 7-68 

7.12.4 Printed Output From the OASYS Subsystem 7-68 

7.12.5 Printed Output From the OABIAS Subsystem 7-68 

7.13 CalComp Plot Tape 7-93 

7. 14 GESS Displays 7-93 

7.14.1 Display Status Flags and Key Assignments 7-93 

7. 14. 2 Array Allocation Sizes . . 7-97 

7. 14. 3 Main Control Display 7-97 

7. 14.4 Options for Copying AE Data Set 7-97 

7. 14. 5 Simulation Options Display 7-97 

7. 14. 6 NAMELIST LIST Display 7-102 

7. 14. 7 Options for Reading Data 7-102 

7. 14. 8 Reader Record Displays ' 7-102 

7. 14. 9 Data Record Display 7-102 

7.14. 10 Data Selection Options for OASYS 7-103 

7.14.11 NAMELIST OPMANl Display 7-112 

7. 14. 12 OASYS Block Average Display 7-112 

7. 14. 13 OASYS Plots of Single Frame Results 7-112 

7. 14. 14 NAMELIST BIASNL Display 7-118 

7.14.15 Final Results From OABIAS 7-118 

7. 14. 16 , Table of Errors From OABIAS .......’ 7-118 

7. 14. 17 Table of Correlation Coefficients From OABIAS 7-126 

7.14. 18 State Component Plots From OABIAS 7-126 

7.14. 19 Plots of Uncertainties in State Components 7-126 

7.14.20 Residual Plots From OABIAS 7-126 

7.14.21 Options for Data Prediction 7-126 

7.14.22 Plot of Predicted and Observed Rotation Angles 7-126 

7.14.23 Plot of Predicted and Observed Earth Widths 7-134 

7.14. 24 Core Storage and Time Remaining Display 7-134 

7.15 Control With Interactive Graphics 7-138 

7. 15. 1 Programmed Function Keys 7-138 

7. 15. 2 Asynchronous Calls 7-140 

7. 16 Error Messages 7-142 

7.16. 1 Messages Displayed in a Graphic Mode 7-142 

7.16. 2 Printed Message From the ODAP Subsystem 7-149 

PRELIMINARY DRAFT 



TABLE OF CONTENTS (Cont'd) 


Section 7 (Cont'd) 

7.16.3 Printed Messages From the OASYS Subsystem 7-155 

7.16.4 Printed Error Message From OABIAS 7-156 

7. 17 Job Control Language 7-162 

VOLUME 4— PROGRAM LISTING AND SAMPLE EXECUTION 
Appendix D - Program Listing 
Appendix E - Sample Execution 


PRELIM INAb DRAFT 



PRELIMiNARY DRAFT 

LIST OF ILLUSTRATIONS 

Figure 


1- 1 MSAD/OABIAS Functional Subsystems 1-5 

2- 1 Sun Sensing Geometxy. . 2-2 

2-2 Horizon Detection Geometry . 2-4 

2- 3 RAE-B PAS Geometry 2-7 

3- 1 Coordinate Frame Flow . 3-37 

3-2 Criteria for Rotation of Frame GI 3-41 

3-3 Geometry of Frames GI, SI, and SC . 3-42 

3-4 Sun Sensor Geometry at Sun Sighting . 3-44 

3-5 Horizon Detector Geometry at Central-Body- In . 3-46 

Crossing . 3-46 

3-6 Generation of Earth Radius Bias Ap by Sensor 

Triggering Level Bias . 3-49 

3-7 Geometrj'’ for Model 1 — Sun Angle Model . 3-57 

3-8 Geometry for Model 3 — Nadir Vector Projection Model . . . 3-62 

3-9 Geometry for Model 4 — Horizon Crossing Time Model . . . . 3-67 

3-10 Geometry for Model 5 — Sun to Earth-In and Sun to 

Earth-Out Dihedral Angle Model ; . 3-71 

3-11 Geometry for Model 6 — Earth Width Model 3-76 

3-12 Geometry for Model 7 — Small Target Model 3-78 

3-13 Geometry for Model 8 — Sun to Earth Mid-Scan Dihedral 

Angle Model 3-80 

3-14 Pg Computation Geometry — Spherical Central Body ..... 3-84 

3- 15 Resolution of Sign Ambiguity of Horizon Crossing ' 

Vector D 3-90 

4- 1 External Interfaces for the MSAD/OABIAS System ...... 4-2 

5- 1 MSAD/OABIAS Baseline Diagram 5-2 

5-2 DHMOD Flowchart 5-59 

5-3 DIAMOD Flowchart 5-73 

5-4j FRAPRO Flowchart 5-103 

5-5 INITIL Flowchart 5-120 

5-6 LNMOD Flowchart 5-146 

5-7 LRMOD Flowchart 5-165 

5-8 OABIAS Flowchart 5-176 

5-9 OADRIV Flowchart 5-185 

5-10 PRINTB Flowchart 5-287 

5-11 SANMOD Flowchart 5-315 

5-12 SCBMOD Flowchart 5-329 

5-13 STMMOD Flowchart 5-343 


xii 

PRELIMINARY DRAFT 



LIST OF ILLUSTRATIONS (Cont'd) 


Figure 


7-1 

7-2 

7-3 

7-4 

7-5 

7-6 

7-7 

7-8 

7-9 

7-10 

7-11 

7-12 

7-13 

7-14 

7-15 

7-16 

7-17 

7-18 

7-19 

7-20 

7-21 

7-22 

7-23 

7-24 

7-25 

7-26 

7-27 

7-28 

7-29 

7-30 

7-31 

7-32 

7-33 

7-34 

7-35 

7-36 

7-37 

7-38 

7-39 


NAMELIST CONTROL 

NAMELIST MAIN . 

Input to O ALIAS Table 

OABIAS Attitude Summary — Heading Definition 

Attitude and State Vector Table From OABIAS 

Uncertainty Table From OABIAS . . 

Partial Derivative Table From OABIAS 

Recursive Estimator Gain Table From OABIAS . . 

Attitude, Observable, and Vector Table From OABIAS . . • 

Residual Table From OABIAS 

Error Messages From OABIAS 

Error Summary Count Table From OABIAS 

OABIAS Filter Response in Right Ascension of Spin Axis . . 

OABIAS Filter Response in Declination of Spin Axis 

Estimated Uncertainty in Right Ascension of Spin Axis. . . . 
Estimated Uncertainty in Declination of Spin Axis ....... 

Residuals of Sun Angle Model 

Residuals of Single Horizon Dihedral Angle Model. ...... 

Display Status Flags and Key Assignments 

Array Allocation Sizes . . 

Main Control Display . . 

Options for Copying AE Data Set 

Simulator Options 

NAMELIST LIST . . 

Options for Reading Data 

Header Record Displays 

Data Record Displays 

Data Selection Options for OASYS . . 

NAMELIST OPMANl 

OASYS Block Average Results 

Alphas Versus Frame Number, Including Rejected Points . 

NAMELIST BIASNL . . 

Final Results From OABIAS . . 

Table of Errors From OABIAS . . 

Correlation Coefficients . . 

Right Ascension Versus Frame Number . . 

Uncertainty in Right Ascension Versus Frame Number . • . 

Residuals From Sun Angle Model . . 

Options for Data Prediction . . 


7-67 

7-69 

7-71 

7-74 

7-75 

7-77 

7-79 

7-80 

7-81 

7-82 

7-84 

7-86 

7-87 

7-88 

7-89 

7-90 

7 - 91 , 

7-92 

7-94 

7-98 

7-99 

7-100 

7-101 

7-103 

7-105 

7-106 

7-109 

7-111 

7-113 

7-116 

7-119 

7-120 

7-123 

7-125 

7-127 

7-128 

7-129 

7-130 

7-131 


XI 11 


PRELIMINARY DRAFT 



LIST OF ILLUSTRATIONS rCont'd) 


Figure 


7-40 Predicted and Observed Rotation Angles Versus Time .... 7-135 

7-41 Predicted and Observed Earth Widths Versus Frame 


Number 7-136 

7-42 Core Storage and Time Remaining Display 7-137 

7-43 Error Messages From ODAP 7-150 

7-44 JCL to Compile, Link, and Execute MSAD/OABIAS ..... 7-163 

7-45 JCL to Execute MSAD/OABIAS 7-171 


LIST OF TABLES 

Table 


1-1 Dependence of Observation Models on State Vector 

Components ; 1-12 

3-1 Coordinate Frames, Axes, and Unit Vectors 3-39 

3-2 Weighting Factor Equation Used in OABIAS 3-95 

5-1 Identification of the Elements of the Arrays Belonging 

to the State Vector Class of Arrays 5-25 

5-2 Identification of the Elements of the Arrays Belonging 

to the Model Class of Arrays * 5-26 

5-3 Identification of the Elements of the Arrays Belonging 

to the Diagnostic Output Control Class of Arrays 5-27 

5-4 Interpretation of Diagnostic Output Control by Array 

IDIAGN in COMMON INBIAS 5-28 

7-1 OABIAS Error-Message Error - Table Correlation 

Matrix 7-158 


.xiv 

PRELIMiMARY DRAFT 



PRELIMINARY DRAFT 


SECTION 7 - OPEILA.TING GUIDE 


7.1 INTRODUCTION 

The purpose of Section 7 is to provide the user with the information required 
to operate the MSAD/OABIAS System in both graphic and nongraphic modes. 

The reader is assumed to be familiar with the material in Sections 1, 2, and 3, 
but not necessarily familiar with Sections 4, 5, and 6. 

7.2 RESOURCES 

The MSAD/OABIAS System was designed to run on the IBM System 360 
Model 95. If all system options are exercised and the suggested overlay struc- 
ture is used, approximately 450K bytes of core storage are required. 

A 2314 disk pack is required for the GESS nonresident tables and for the pro- 
gram load module. A 2250 or 2260 display device is required for operation in 
a graphic moc e. The OABIAS data set will normally reside on disk or tape; in 
a simulation mode a scratch disk or scratch tape can be used. The Atmosphere 
Explorer (AE) data set, if used, must reside on disk or tape. Control param- 
eters are normally read from cards and output routed to a lineprinter, but 

< 

job control language (JCL) modifications can alter this at the option of the user. 
The remainii^ data sets (attitude tape, CalComp plot tape, and ephemeris files) 
are optional and may reside on any stroage device (e.g. , tape, disk, or data 
cell). 

7. 3 CARD INPUT; NAMELIST PARAMETERS 

• 

All of the following card input is normally read from FORTRAN unit 5, with the 
exception of NAMELIST CONTRL and the optional GESS control cards, which 
are read from FORTRAN unit 97. However, in a nongraphic mode the user 
may specify the FORTRAN unit number for NAME LIST BIASNL using the pa- 
rameter REDUNT. The first NAMELIST BIASNL in each set must appear on 
unit 5. ~ 
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Card input should be ordered as follows for a job in a nongraphic mode: 
FORTRAN unit 97: 

1. NAMELIST CONTRL 

2. GESS array allocation sizes (optional; read only if IRDART = 1 
in NAMELIST CONTROL; normally omitted) 

3. GESS display status flags (optional; read only if IRDXST = 1 in 
NAMELIST CONTRL; normally omitted) 

FORTRAN unit 5: 

1. NAMELIST MAIN 

2. NAMELIST LIST (optional; read only if sumulation is specified 
(ISIM = 1) in NAMELIST 3VLAIN) 

3. NAMELIST OPMANl 

4. NAMELIST BIASNL 

5. Additional sets of NAMELIST parameters for NAMELIST BIASNL; 
read only if REDUNT (1) ^ 0 in the previous NAMEIJST BIASNL 

6. NAMELIST MAIN 

7. , NAMELIST LIST 

8. NAMELIST OPMANl 

9. NAMELIST BIASNL 

10. Additional NAMELIST sets, ordered the same as items 6-9, above 

For the graphic mode, the same set of card input may be used, except that all 
card input following item 4 will be ignored. 

Note that NAMELIST parameters are not reset to their default values after 
processing the first set. Thus, in a nongraphic mode, additional sets of 
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NAMELIST parameters need specify only those parameters which are to be 
changed from the previous set. A description of each NAMELIST follows. 

Each NAMELIST has been divided into logical subsets for purposes of descrip- 
tion. 

\ 
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CONTRL 


7.3.1 NAMELIST CONTRL 

NAMELIST CONTRL is read by the GESS Executive and is read only once in 
a job. 

Parameters in NAMELIST CONTRL are as follows: 


Name 

Type 

Default 

Description 

IFTABL 

1*4 

96 

FORTRAN unit number for GESS non- 
resident tables 

IFTUBE 

1*4 

23 

FORTRAN unit number for graphics de- 
vice 

IFTPRT 

1*4 

6 

FORTRAN unit number for GESS 
printout 

IRDART^ 

1*4 

0 

Array allocation sizes reading option: 
= 0, do not read array allocation sizes 
= 1, read array allocation sizes 

IRDXST^ 

1*4 

0 

Display status flags reading option: 
= 0, do not read display status flags 
= 1, read display status flags 

GSATID 

R*8 

»*» 

Character string for satellite identifica- 
tion 

GRUNID 

R*8 


Character string for run* identification 


IDIREC , 

IRDTPD I : 

NUMSUB I 
NUMCNC I 

NUMSCA y Not applicable to MSAD/OABIAS 

KOFFEE I 

IRYEAR I 

IRMON I 

IRDAY ' 


^See Section 7,3.2. 
2 

See Section 7.3.3. 
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CONTRL 


7.3.2 GESS Array Allocation Sizes 

GESS array allocation sizes may be read from cards to change the default ar- 
ray allocation sizes in a nongraphic mode. See Reference 13, page 4-25 for 
a description. This input is needed only if IRDART = 1 in NAMELIST CONTRL. 
These array allocation sizes are normally omitted. 

7.3.3 GESS Display Status Flags 

GESS display status flags may be read from cards to change the display status 
flags in a nongraphic mode. See Reference 13, page 4-25 for a description. 

This input is needed only if IRDXST = 1 in NAMELIST CNTRL. These display 
status flags are normally omitted. 


7-5 

PRELIMINARY DRAFT 



PRELIMINARY DRAFT 

MAIN 

Subsystem Selection 

7.3.4 NAMELIST MAIN 

NAMELIST IVIAIN is read by OADRIV, the main control module of MSAD/ 
OABIAS. 

The parameter for subsystem selection is as follows; 

Name Type Default Description 

ISIM 1*4 1 Subsystem selection indicator: 

= 1, simulate data 

= 2, read data from OABIAS data set 
= 3, display data already in core 
= 4, copy AE data set to OABIAS data 
set 

=-5, terminate 
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MAIN 

Simulation 


The following parameters apply only for simulation (ISIM = 1): 

Name 

Ty£e 

Default 

Description 

NFRAME 

1*4 

200 

Number of frames to simulate 

IREWND 

1*4 

1 

Rewind indicator for simulated data set; 
= 0, do not rewind before writing simu- 
lated data 

= 1, rewind before writing simulated 
data 

OMEGAl(2) 

R*4 

100., 0.0 

Spin rate coefficients (degrees per 
second, degrees per second^) 

TOMEGA(6) 

R*4 

6*0.0 

Reference time for spin rate (universal 
time (UT): year, month, day, hours, 
minutes, seconds). The spin rate at 
time t is computed as: 

OMEGAl(l) + (t - TOMEGA) * 
OMEGAl(2) 

IDISK 

1*4 

50 

FORTRAN unit number for simulated 


OABIAS data set 
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MAIN 

Copying AE Data Set 


The following parameters apply only to copying the AE data to the OABIAS data 


set aSIM = 4); 
Name 

Type 

Default 

Description 

IBLKl 

1*4 

1 

First block on AE data set to be copied 

IBLK2 

1*4 

99999 

Last block on AE data set to be copied 

lAEUN 

1*4 

49 

FORTRAN unit number of the AE data 
set 

ICSKIP 

1*4 

1 

Record skipping indicator. Copy every 
ith wheel horizon sensor frame on the 
AE data set, where i = ICSKIP 

AESENS(4) 

R*4 


Sensor mounting angles for AE (degrees) 

(1) 


104.0 

Wheel horizon sensor 1 

(2) 


118.0 

Wheel horizon sensor 2 

( 3 ) 


110.0 

Body horizon sensor 1 or 2 

( 4 ) 


110.0 

Not used 

IBSKIP 

1*4 

1 

Record skipping indicator for body hor- 
izon sensor data. Copy every ith body 
horizon sensor frame on the AE data 
set, where i = IBSICIP 

BHSREF 

R*4 

90.0 

Azimuth of AE body horizon sensor 
with respect to Sun sensor (degrees) 

ISUNSP 

1*4 

0 

Inertial spin period indicator: 

= 0, use spin periods reported with data 
= 1, use Sun times to generate inertial 
spin periods 

OLDRPM 

R*4 

0.0 

Initial inertial spin period value re- 
ported for first frames until valid 
inertial period available from Sun times. 
Valid only if ISUNSP = 1 

IBOLOF 

1*4 

0 

OABIAS AE data set indicator: 

= 0, do not write OABIAS AE data set 
= 1, write OABIAS AE data set on 

FORTRAN unit number specified by 
IBOLOU 

IBOLOU 

1*4 

49 

FORTRAN unit number for OABIAS AE 
data set 
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MAIN 

Reading Data 


The following parameters apply to reading the OABIAS data set (ISIM = 2 ): 


Name 

Type 

Default 

Description 

IRE Cl 

1 * 4 : 

1 

Record number of first record to read 

IREC2 

1*4 

99999 

Record number of last record to read 

lADD 

1*4 

0 

Option for combining newly read data 


with data already in core: 

= 0, replace old data 
= 1, add new data to old data (concate- 
nate) 

= 2, compress old data, deleting frames 
which are flagged; then add new 
data to old data 


ISKIP 

;i*4 

1 

Record skipping indicator. Select 
every ith record, where i = ISKIP 

IPLOT 

1*4 

0 

Rotation angle plot indicator; 

= 0, plot Earth-in and Earth-out on 
same plot 

= 1, plot Earth-in and Earth-out on 
separate plots only; used to avoid 
exceeding 2250 buffer size 

lEFMFG 

1*4 

0 

Ephemeris vector indicator; 


=r 0, do not preserve ephemeris vectors 
with data 


= 1, preserve ephemeris vectors with 
data ('input data records ' master 
number must be 200 or less; data 
sent to Optical Aspect Attitude 
Determination System (OASYS) 
must start with frame number 1) 

IHEADR 1*4 3 Selection indicator for header records 

and/or data records: 

^ 1, save data records only 
= 2, save header records only 
= 3, save both header records and data 
records 
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Reading Data 


Name 

Type 

Default 

Description 

ICBSEL{4) 

1*4 


Array specifying central body selection 


based on central body flag on OABIAS 
data set: 

= 0, do not select this central body 


(skip records with this central body) 
= 1, select this central body 


(1) 


1 

Earth 

(2) 


1 

Moon 

( 3 ) 


0 

Sun 

( 4 ) 


1 

Unidentified 

ICNSEL(8) 

1*4 

8*1 

Array specifying selection of configura- 
tion flag: 

ICNSEL(i) = 0, do not select configura- 
tion flag i (skip records 
with this configuration) 

= 1, select configuration flag i 
(the sensor configurations 
are described in Sec- 
tion 7.4) 

IDISK 

1*4 

50 

FORTRAN unit number for OABIAS data 
set 

Kxrxx 

1*4 

0 

Periodic reading indicator 
= 0, read all data 
= 1, read data periodically 
PERIOD, BNDWTH, and TREF6 apply 
only if IXXIXX = 1. If IXXIXX = 1, then 
only data in the intervals T = TREF6 + 
N*PERI0D to TREF6 + N*PERIOD + 
BNDWTH are read. N is any integer, 
positive or negative 

PERIOD 

R*4 

0.0 

Period for periodic reading (see above) 

BNDWTH 

R*4 

0.0 

Interval width for periodic reading (see 
above) 

TREF6(6) 

R*4 

0.0 

Reference time for periodic reading 
(j^ear, month, day, hour, minute, 
second) (see above) 
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OASYS Processing 


The following parameters apply only to selecting a block of data for processing 
by OASYS (ISIM = 3); 


Name 

Type 

Default 

Description 

NFRAMl 

1*4 

T 

First frame number of block of data to 
pass to OASYS 

NFRAM2 

1*4 

200 

Last frame number of block of data to 
pass to OASYS 

lEPHEM 

1*4 

1 

Option for using ephemeris data from 
OABIAS data set: 

= 0, do not use ephemeris data from 
OABIAS data set 

= 1, use ephemeris data from OABIAS 
data set (data set must have been 
read with lEFhlFG = 1) 

ITERM 

1*4 

1 

Option for using terminator flags on 
OABIAS data set: 

= 0, do not use terminator flags 
= 1, use terminator flags (reject all 
triggerings identified as termina- 
tors) 

ICB 

1*4 

1 

Not used 

TADJ 

R*4 

0.0 

Time adjustment (seconds). Added to 
all times to be passed to' OASYS 

INOISE 

1*4 

3 

Option for adding noise to data to be 


passed to OASYS: 

= 1, process data twice, once without 
noise and once with noise (applies 
in nongraphic mode only) 

= 2, process data once, with noise 
= 3, process data once, without noise 
NOTE: If INOISE = 3, then the re- 
maining parameters in this sec- 
tion do not apply 

STDV(6) R*4 6*0. 0 Standai*d deviation of Gaussian noise 

to be applied to data 

(1) Telemetry frame time (seconds) (not 

used) 

(2) Sun time (seconds) 

(3) ' Time in (seconds) 
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MAIN 

OASYS Processing 


Name 

Type 

Default 

Description 

STDV(4) 

(5) 

(6) 

RM 

6*0.0 

Time out (seconds) 

Spin rate (degrees per second) 
Sun angle (degrees) 

BIAS(6) 

R*4 

6*0.0 

Biases to be added to data (same order 
and units as for STDV) 

PROB(6) 

R*4 

6*0.0 

Probability of telemetry error (bit 
change); same order as for STDV 

IRAND 

1*4 

123456789 

Initial number for random number gen- 
erator (any odd integer) 

BGAM 

RM 

0.0 

Bias to be added to each sensor mounting 
angle (degrees) 

SQUANT 

R*4 

0.0 

Sun angle quantization (degrees): 

= 0. 0, no effect 

^ 0. 0, each Sun angle is quantized as 
follows: 

= OS/SQUANT) * SQUANT + 
SQUANT/2. 0 

where /3 is the am angle, and 
the quantity in parentheses is 
truncated to an integer 

ABETAl 

R*4 

0.0 

Sun ai^e coefficient for linear fit 
(degrees) (see note below) 

ABETA2 

R*4 

0.0 

Rate of chaise of Sun angle (degrees 
per hour) (see note below) 

TBETA(6) 

R*4 

6*0.0 

Reference time for Sun angle fit 


(UT: year, month, day, hours, 
minutes, seconds) 

NOTE: If ABETAl = 0. 0, and 


ABETA2 = 0. 0, then the Sun 
angle from the OABIAS data set 
is used, and the Sun angle fit 
parameters (ABETAl, ABETA2, 
TBETA) are ignored. 

If ABETAl 5 ^ 0.0, or 
ABETA2 7 ^ 0.0, then the Sun 
angle computed as follows 
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MAIN 

OASYS Processing 


Name Type Default 


Description 


TBETA(6) R’*=4 6*0. 0 replaces the Sun angle from the 

(Cont'd) OABIAS data set: 

/3(t) = ABETAl + ABETA2* 

(t - TBETA) 

where jS(t) = the Sim angle at 
time t , 

If TBETA is earlier than 
September 1, 1957 (e. g. , if 
TBETA = 6*0.0), then the in- 
itial time from the block of 
data is used in place of TBETA 
in the above expression. 


OMEGl 

R*4 

0.0 

Spin rate coefficient for linear fit 
(degrees per second) (see note below) 

OMEG2 

R*4 

0.0 

Rate of change of spin rate (degrees 
per second^) (see note below) 

TOMEG (6) 

R*4 

6*0.0 

Reference time for spin rate (UT: 
year, month, day, hours, minutes, 


seconds) 

NOTE: If OMEGl = 0.0 and OMEG2 = 
0.0, the spin rates from the 
OABIAS data set are used. If 
OMEGl 0. 0 or OMEG2 ^ 0. 0, 
then the spin rate computed 
as follows replaces the spin 
rate from the OABIAS data set; 
oi (t) = OMEGl + OMEG2* 

(t - TOMEG) 

where a:(t ) = spin rate at time t 
If TOMEG is earlier than 
September 1, 1957, then the 
initial time from the block of 
data is used in place of TOMEG. 
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MAIN 

Data Prediction 


The following parameters apply only to the data prediction subsystem, PLOTOC. 
Parameters whose description begins with an asterisk (*) have three compo- 
nents which refer to plot 1, 2, and 3, respectively. 


Name Type Default 


Description 


I STATE (3) 1*4 1,2,0 *State selection indicator for plot of 

predicted and observed rotation angles 
and scan widths: 

= 0, no plot 

= 1, use specified values of ALF, DEL 
(see note below) and biases 
= 2, use specified values for BETA, 

PHI (see note below) and biases 
= 3, use OABIAS final state for attitude 
and biases 

= 4, use OASYS block average results 
for attitude and the input biases 
used in OASYS 
>4, repeat old plot 
NOTE; When returning to the options 

for data prediction display from 
within PLOTOC, the subroutine . 
automatically returns 11, 22, 

33, or 44 to the "Use State" 
(=ISTATE) for each plot to in- 
dicate whether the previous plot 
was 1, 2, 3, or 4, respectively. 
The only user options are to re- 
display the arrays previously 
displayed in that plot or to gen- 
erate a new plot or to generate 
no plot. 


ALF(3) 

R*4 

0.0 

♦Right ascension (geocentric inertial 
(G.I. )) (degrees) 

DEL(3) 

R*4 

0.0 

♦Declination (G.I.) (degrees) 

BETA (3) 

R*4 

0.0 

♦Sun angle (degrees) between Sun vector 
and spin axis at block start time 

PHI (3) 

R*4 

0.0 

♦Phase of attitude on Sun cone (degrees) 
at block start time; defined as the di- 


hedral angle from the plane defined 
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Data Prediction 


Name 

Type 

Default 

Description 

PHI (3) 
(Cont'd) 

R*4 

0.0 

by the Sun vector and the Z-axis (G.L ) 
to the plane defined by the Sun vector 
and the spin axis 

BSIGMA (3) 

R*4 

0.0 

’•'Bias to be added to sensor mounting 
angle (degrees) 

BRHO(3) 

R*4 

0.0 

*Bias to be added to the angular radius 
of the central body (degrees) 

BAl (3) 

R*4 

0.0 

*Bias to be added to observed in- 
triggering (degrees) 

BA2(3) 

R*4 

0.0 

’•’Bias to be added to observed out- 
triggering (degrees) 

BTIME (3) 

R*4 

0.0 

’•‘Time adjustment to be added to telem- 


etry times (seconds) 

NOTE: For all values of ISTATE, the 
program fills in the values ap- 
propriately for ALF, DEL, 
BETA, PHI, BSIGMA, BRHO, 
BAl, BA2, and BTIME, and 
these values can be examined 
the next time this display 
appears. For example, if 
ISTATE = 1, then the attitude 
specified by AL!f and DEL is 
converted to the coordinate 
system for BETA and PHI, and 
BETA and PHI are set. If 
ISTATE = 2, then BETA and PHI 
are converted to ALF and DEL. 
If ISTATE = 3, then the OABIAS 
final state is used to set ALF 
and DEL, and this attitude is 
converted to set BETA and 
PHI. The biases (BSIGMA, 
BRHO, BAl, BA2, and BTIME) 
are set to the corresponding 
bias values from the OABIAS 
final state. 
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MAIN 

^ ' Data Prediction 

Name 

Type 

Default 

Description 

ICENT(3) 

P4 

1 

♦Central body indicator: 

= 1, Earth is central body 
= 2, Moon is central body 

RADE(3) 

R*4 

6378.16 

♦Effective radius of Earth (kilometers) 

RADM(3) 

R*4 

1738. 

♦Effective radius of Moon (kilometers) 

DANG E (3) 

R+4 

89.7 

♦Dark angle of Earth (degrees) 

DANGM(3) 

R*4 

89.7 

♦Dark angle of Moon (degrees) 

I0BLAT{3) 

1*4 

1 

♦Oblateness indicator: 

= 0, use spherical Earth with radius 
RADE + HT 
= 1, use oblate Earth 

HT(3) 

R*4 

0.0 

♦Height of atmosphere above Earth 
surface (kilometers) 

IPSKIP(3) 

1*4 

1 

♦Frame skipping indicator. Generate 
two frames of predicted data for every 
(i 4- 2)th frame, where i = IPSKIP (used 
to reduce execution time) 

I WIDTH 

1*4 

0 

Scan width indicator: ^ 

= 0, plot both rotation angles and scan 
widths 

= 1, plot scan width onlj^ 

ITIME 

1*4 

0 

Parameter for selecting equal time in- 
crement plotting (applies to all plots): 

= 0, plot at the frame times 
= 1, plot at NPTS equal frame points 
between TSTAET and TEND (see 
note below) 

NOTE: NPTS, TSTART, and TEND 
apply only if ITIME = 1 

NPTS 

1*4 

200 

Number of points for equal time in- 
crement plotting; must be 200 or less. 
If NPTS = 0, the number of frames of 
observed data is used 




7-16 , . 


PRELIMINARY DRAFT 



PRELIMI^3ARY DRAFT 

MAIN 
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Name Type Default Description 


TSTART{6) 

R*4 

0.0 

Start time for data prediction (UT; 
year, month, day, hours, minutes, 
seconds). If T START is earlier than 
September 1, 1957, then the start time 
of the observed data is used 

TEND (6) 

R*4 

0.0 

Final time for data prediction (UT: 
year, month, day, hours, minutes, 
seconds). If TEND is earlier than 
September 1, 1957, then the end time 
of the observed data is used 

IORBC(3) 

1*4 

1 

Orbital motion correction indicator: 

= 0, do not correct for orbital motion 
within a spin period. Compute 
prediction based on spacecraft 
position at the in-triggering time 
= 1, correct for orbital motion during 
a spin period. Compute prediction 
for in-triggering based on space- 
craft position at the in-triggering 
time and prediction for out- 
triggering based on spacecraft 
position at the out-triggering time 




NOTE: Execution time for predictions 
is approximately doubled if 
lORBC = 1 

ISAME 

1*4 

1 

Duplicate parameters indicator: 

= 0, use specified parameters for each 
plot 

= 1, use plot 1 parameters (ICENT on 
only) for plot 2 and plot 3 

IDATA 

1*4 

1 

New data plot parameter: 

= 0, use old data plot 
= 1, generate new data plot 

IVSTAT(3) 

1*4 

1,2.0 

*IVSTAT(i) retains the value of ISTATE (i) 
used to generate the current plot i 

ISAVTM 

1*4 

0 

ISAVTM retains the value of I TIME used 


to generate the current plots 
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Data Prediction 

Name Type Default Description 

OMEGAP R*4 0 Spin rate used to generate time for 

predicted data for flagged frames. 
Automatically set to spin rate calcu- 
lated by OABIAS (degrees per second) 

NRDF 1*4 1 Number of revolutions per frame of data; 

(used to generate time for predicted 
data for flagged frames) 
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LIST 

Attitude 


7.3.5, NAMELIST LIST 


NAMELIST LIST is read by subroutine ODAPIN and contains parameters spec- 
ifying the simulation conditions. 

Spacecraft attitude parameters are as follows: 


Name 

Type 

Default 

Description 

ALPHA 

RM 

0.0 

Right ascension (G. I. ) of the spin axis 
(degrees). Valid range: 0 ^ ALPHA < 
360. This parameter does not apply if 
lATAPE = 0 

DELTA 

R*4 

0.0 

Declination (G. I. ) of the spin axis (de- 
grees). Valid range: -90 ^ DELTA ^ 
90. This parameter does not apply if 
LATAPE = 0 

lATAPE 

1*4 

0 

Attitude tape indicator: 

= 0, get attitude from attitude tape 
= 1, use fixed ALPHA and DELTA read 
from cards 
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Name 
TO (6) 


TF(6) 
BELT AT 

INTER P 
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LIST 

Time 


time parameters are as follows: 


Type Default 


Description 


R*4 6*0.0 


R*4 6*0.0 

R*4 60.0 


1*4 0 


Starting time for data simulation (UT: 
year, month, day, hours, minutes, 
seconds). If TO is not specified and 
lATAPE = 0, then the starting time is 
taken as the start time of the attitude 
tape. If LATAPE = 1, TO must be 
specified 

Not used; should not be specified 

Time interval at which simulated data 
will be generated (seconds). Valid 
range; DELTAT > 0 . This parameter 
does not apply if lATAPE = 0 and 
INTERP = 0 

Interpolation indicator for attitude tape: 

= 0, use the time spacing on the attitude 
tape, without interpolation (i. e. , 
start at the first time on the atti- 
tude tape greater than or equal to 
TO and use each suceeding time). 

Note that the time spacing on the 
attitude tape need not be uniform 
= 1, use linear interpolation with the 
attitude tape to obtain the attitude 
at the desired times, using the 
specified value of DELTAT. 

INTERP does not apply if lATAPE = 1 
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Orbit 

Simulation orbital elements (in the arrays dimensioned (2), subscript = 1 ap- 
plies to spacecraft orbit around Earth or Moon, subscript = 2 applies to orbit 
of Moon around Earth) are as follows: 

Name Type Default Description 

A(2) R*4 2*0.0 Semimajor axis (kilometers) 

E(2) R*4 2*0.0 Eccentricity (unitless) 

EYE (2) R*4 2*0.0 Inclination (degrees) 

EMO(2) R*4 2*0.0 Mean anomaly at epoch time (degrees) 

WO(2) R*4 2*0.0 Argument of perigee (degrees) 

RANODE(2) R*4 2*0.0 Right ascension of ascending node (de- 

grees) 

ORBITE(6) 

ORBITM(6) 

ISDN 


ISPC 1*4 2 Spacecraft ephemeris indicator: 

= 1, use ORBGEIN to generate space- 
craft orbit, using orbital elements 
for an orbit around the Earth (see 
module description of ORBGEN) 
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R*4 6*0.0 Epoch time of orbital elements for 

spacecraft (UT: year, month, day, 
hours, minutes, seconds) 

R*4 6*0. 0 Epoch time of orbital elements for 

Moon (UT: year, month, day, hours, 
minutes, seconds) 

1*4 1 Sun ephemeris indicator: 

= 1, use SUNl to obtain Sun position 
(equinox of date) (see Reference 6) 
= 2, use SUNRD to obtain Sun position 
(uses a direct-access data set on 
FORTRAN unit 14; either equinox 
of date or equinox of 1950. 0 de- 
pending on data set accessed) (see 
Appendix B) 

= 3, use RJPLT to obtain Sun position 
(accesses a Jet Propulsion Labo- 
ratory (JPL) ephemeris data set 
on FORTRAN unit NRJPLT; equi- 
nox of date) (see Reference 5) 
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LIST 

Orbit 



Name 


Type 


Default 



ISPC 1*4 2 =2, use DTAPRE to obtain spacecraft 

position (accesses a standard Defi- 
nitive Orbit Determination System 
(DODS) EPHEM data set on 
FORTRAN unit NORBl; must be 
Earth-centered ephemeris) (see 
Reference 3) 

= 3, use GE THDR/ GE TV C T to obtain 
spacecraft position (accesses a 
Goddard Trajectory Determination 
System (GTDS) ephemeris file, 
either direct-access or sequential, 
depending on the value of LEVEL 
(see below); the ephemeris may be 
either Earth -centered or Moon- 
centered; the header record of the 
data set specifies the origin for the 
coordinates) (see Reference 4) 

= 4, use ORBGEN to generate spacecraft 
orW t, using orbital elements for an 
orbit around the Moon (see module 
description of ORBGEN) 

= 5, use DTAPRE to obtain spacecraft 
position (accesses a standard DODS 
EPHEM data set on FORTRAN unit 
NORBl; must be Moon-centered 
ephemeris) (see Reference 3) 

IMOON 1*4 0 Moon ephemeris indicator: 

= 0, do not obtain Moon position. This 
option must not be chosen if ICB >2, 
or if a Moon-centered ephemeris 
file is used for GETHDR/GETVCT, 
or if ORBGEN is used for a space- 
craft orbit around the Moon 
(ISPC = 4), or if DTAPRE is used 
with a Moon-centered EPHEM tape 
(ISPC = 5) 

= 1, use ORBGEN to generate orbit of 
Moon using orbital elements for an 
orbit around the Earth (see module 
description of ORBGEN) 
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Orbit 


Name 

Type 

Default 

Description 

IMOON 

1*4 

0 

= 2, use SUNRD to obtain Moon position 


(accesses a direct-access data set 
on FORTRAN unit 14; either equinox 
of date or equinox of 1950. 0, de- 
pending on data set accessed) (see 
Appendix B) 

= 3, use RJPLT to obtain Moon position 
(accesses a JPL ephemeris data 
set on FORTRAN unit NRJPLT; 
equinox of date) (see Reference 5 ) 

NOTE; The orbital elements and epoch 
time for the spacecraft must be 
specified if ISPC = 1 or 4. Or- 
bital elements and epoch time 
for the Moon must be specified 
if IMOON = 1. 


NRJPLT 

1*4 

28 

FORTRAN unit number for JPL ephem- 
eris data set, read by RJPLT 

NORBl 

1*4 

30 

FORTRAN unit number for DODS EPHEM 
data set, read by DTAPRE 

NGTDS 

1*4 

29 

FORTRAN unit number for GTDS ORBIT 
file, read by GETHDR/GETVCT (only if 
a sequential file is used, .LEVEL = 0; 
see below) 

LEVEL 

1*4 

0 

Level number for GETHDR/GETVCT: 

= 0, .use sequential orbit file, FORTRAN 
unit number NGTDS 


> 0, level number on direct-access orbit 
file, FORTRAN unit number 31 
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Name 

IPRINT 


IREPOR 
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Printout 


which control simulation printout are as follows: 

Type Default Description 

1*4 0 Print indicator: 

= 0, print every output line (i. e. , do 
not suppress printing) 

= 1, suppress printout of an output line 
if the data condition flag for this 
line is the same as the data flag 
for the previous line (i.e,, print 
an output line only for those times 
at which the data flag changes). 

The output line for the first time 
in the prediction interval is always 
printed, regardless of the value of 
the data flag. The output line for 
the final time is also printed, un- 
less this problem terminates ab- 
normally 

1*4 0 Report generation indicator: 

= 0, generate final report summary 
= 1, do not generate final report 
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Control 


Simulation control parameters are as follows: 


Name 

Type 

Default 

Description 

SIGMA 

R*4 

90.0 

Initial value for the sensor mounting 

— 



angle between the spin axis and the 


optical axis of the horizon telescope 
(degrees). Valid range: 0 ^ SIGMA 
^ 180. The sensor mounting angle is 
incremented by the value STEP after 
each frame or after each spin period 
if no sensor triggering occurs 


STEP 

R*4 

0.703125 

Step size for panoramic scanner (de- 
grees). If STEP = 0.0, the sensor 
mounting angle will remain constant 
at the value SIGMA 

BIASRE 

R*4 

0.0 

Bias to be added to angular radius of 
Earth (degrees) 

BIASRM 

R*4 

0.0 

Bias to be added to angular radius of 
Moon (degrees) 

lOBLAT 

1*4 

1 

Model for Earth radius 


= 0, use spherical Earth, of radius = 


RADE+HT 

= 1, use oblate Earth, with atmosphere 
height HT above surface 


RADE 

R*4 

6378.16 

Earth radius (kilometers); used only 
if lOBLAT = 0 

HT 

R*4 

0.0 

Height of atmosphere layer above 
surface (kilometers) 

ICB 

1*4 

3 

Central body indicator: 

= 1-2, consider Earth only as central 


body 

= 3-6, consider both Earth and Moon as 
central bodies. Note that there 
is no option provided to consider 
the Moon as the only central 
body; however, if data for the 
Earth as central body are not 
desired, the dark angle of the 
Earth, DANGE, may be set to 
360. 0 degrees, causing the 
Earth to be ignored 


— 7-25 

PRELIMINARY DRAFT 



Name 

OBLATl 

0BLAT2 


EPS 

THETAC 

DANGE 

DANGM 
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Control 


Type Default 


Description 


R*4 

R*4 


R*4 


R*4 


R*4 


R*4 


0.0033528 Earth oblateness parameter (see note 
below) 


0.0 


0.0 


45 


90.0 


90.0 


Earth oblateness parameter for asym- 
metric Earth (see note below) 

NOTE; The Earth radius is computed 
using the formula 
R = 6378. 16 (1. 0 - OBLATl 
* sin^ 0 + OBLAT2 * sin 0) 

+ HT 

where R = effective Earth radius 
in kilometers 
0 = geocentric latitude 

Angular field of view of sensor (degrees). 
Valid range: 0 ^ EPS ^ 180. This pa- 
rameter is used to determine if the sen- 
sor field covers the entire view of the 
central body or if part of the sensor field 
hits the central body 

Sensor cutoff angle (degrees). The sen- 
sor is disabled when the dihedral angle 
from the plane defined by the spin axis 
and the Sun vector to the plane defined 
by the spin axis and the optical axis of 
the horizon telescope is less than THETAC 
or greater than 360 - THETAC. Valid 
range: 0 ^ THETAC ^ 180 

Dark angle of the Earth, in degrees 
measured between the vector to the 
terminator and the negative of the Sun 
vector 

Dark angle of the Moon (degrees). (Note 
that the dark angle is approximately 
90 degrees - hg, where hs = half the 
angle subtended by the Sun = 0.27 degree. 
However, the dark angle may depend on 
the extent of the central body's atmos- 
phere and the sensitivity of the sensor. 

For an infrared sensor, both DANGE 
and DANGM should be set to -90) 
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Name 


DELPHI 


PHITOL 


ISKIP 


Type 

R*4 


R*4 


1*4 


Default Description 


0. 5 Angular increment (degrees) used when 

searching for pairs of terminator cross- 
ir^s. If there are two terminator cross- 
ings for which the spacecraft rotation 
angles are closer then DELPHI, these 
crossings may be overlooked. Valid 
range: DELPHI > 0. Values of 
DELPHI less than 0.40 should not be 
used, unless the routine TER MIN is 
modified to allow more than 1000 steps 
in the terminator search. Larger values 
of DELPHI may be used to reduce ex- 
ecution time 

0.01 Maximum allowable error in the phase 
angle of a terminator crossing. The 
half-interval search is terminated when 
the interval becomes smaller than 
PHITOL. Valid range: PHITOL > 0 . 
Values of PHITOL less than 0. 001 
should not be used, due to the limited 
precision of REAL *4 arithmetic 
(degrees) 

0 Data condition skipping indicator: 

= 0, simulate one frame of data for 
each time point; do hot skip time 
points 

= 1, skip any time point if the data con- 
dition flag for this time point is 
blank (i. e, , do not store the data 
for this time point in the simulated 
data arrays). Continue increment- 
ing the time until a frame is ob- 
tained for which the data condition 
flag is not blank. (This option 
must be used in conjunction with 
a value for IFLAG. See below) 
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Description 

Flag indicator; 

= 1, normal flag function 
= 2, set flag to blanks if no sensor 
triggering occurs for this frame; 
otheiwise, use normal flag func- 
tion 

= 3, set flag to blanks if no sunlit hor- 
izon triggerings occurred; other- 
wise, use normal flag function 
= 4, set flag to blanks if a triggering 
was caused by a terminator cross- 
ing or if no triggerings occurred; 
otherwise, use normal flag function 
= 5, use normal flag function if at least 
one triggering occurred at a ter- 
minator; otherwise, set flag to 
blanks 


7-28 

PRELIMINARY DRAFT 



PRELIMINARY DRAFT 

LIST 

Uncertainties 

Simulation uncertainty parameters are used mainly for simulation summary 
printout to indicate the expected quality of an attitude solution based on a given 
frame of data given the following uncertainties. If all uncertainties equal 0. 0, 
no attitude uncertainty is computed. Simulation uncertainty parameters are as 


follows; 




Name 

Type 

Default 

Description 

ERRBET 

R*4 

0.0 

Uncertainty in Sun angle (degrees) 

ERRGAM 

R*4 

0.0 

Uncertainty in sensor mounting angle 
(degrees) 

ERRA 

R*4 

0.0 

Uncertainty in rotation angle of space- 
craft from Sun crossing to horizon 
crossing (degrees) 

ERR AD 

R*4 

0.0 

Uncertainty in angular radius of the 
central body (degrees) 

ERRTIM 

R*4 

0.0 

Uncertainty in absolute time of horizon 
crossings, or uncertainty in ephemeris 
information (seconds) 

IDE BUG 

1*4 

0 

Debug indicator for uncertainty printout: 
= 0, no debug printout 
> 0, generate debug printout on 

FORTRAN data set reference num- 
ber IDEBUG 
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Orbital Motion Correction 


Parameters controlling corrections for orbital motion in data simulation are 

as follows: 

Name Type Default Description 

TOL R*4 0.0001 Desired tolerance for sensor triggering 

time (seconds). If no iteration is de- 
sired, a large value may be used (e.g., 
l.OE+50, ) so that convergence will 
always occur on the first iteration 

MAXIT 1*4 10 Maximum allowed number of iterations 
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7.3.6 NAMELIST OPMANl 

NAMELIST OPMANl is read by subroutine OPINIT and contains parameters 
specifying the processing options for OASYS. 

Orbital elements (in arrays dimensioned (2), subscript = 1 applies to spacecraft 
orbit around Earth or Moon, subscript = 2 applies to orbit of Moon around 
Earth) are as follows; 


Name 

Type 

Default 

Description 

A(2) 

R*4 

2*0. 

Semimajor axis (kilometers) 

E(2) 

R*4 

2*0. 

Eccentricity (unitless) 

EYE (2) 

R*4 

2*0. 

Inclination (degrees) 

EMO(2) 

R*4 

2*0. 

Mean anomaly (degrees) 

WO(2) 

R*4 

2*0. 

Argument of perigee (degrees) 

RANODE(2) 

R*4 

2*0. 

Right ascension of ascending node (de- 
grees) 

ISUN 

1*4 

1 

Sun ephemeris indicator; 

= 1, use SUNl to obtain Sun position 
(equinox of date) (see Reference 6; 
= 2, use SUNRD to obtain Sun position 


(uses a direct-access data set on 
FORTRAN unit 14; either equinox 
of date or equinox of 1950. 0, de- 
pending on data set accessed) (see 
Appendix B) 

= 3, use RJPLT to obtain Sun position 
(accesses a JPL ephemeris data 
set on FORTRAN unit N RJPLT; 
equinox of date) (see Reference 5) 

ISPC 1*4 2 Spacecraft ephemeris indicator: 

= 1, use ORBGEN to generate space- 
craft orbit, using orbital elements 
for an orbit around the Earth (see 
module description of ORBGEN) 
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Type Default Description 

1*4 2 =2, use DTAPRE to obtain spacecraft 

position (accesses a standard 
DODS EPHEM data set on FORTRAN 
unit NORBl; must be Earth- 
centered ephemeris) (see Refer- 
ence 3) 

= 3, use GETHDR/GETVCT to obtain 
spacecraft position (accesses a 
GTDS ephemeris file, either direct 
access or sequential, depending on 
the value of LEVEL (see below); 
the ephemeris may be either Earth- 
centered or Moon-centered; the 
header record of the data set spec- 
ifies the origin for the coordinates) 
(see Reference 4) 

= 4, use ORBGEN to generate space- 
craft orbit, using orbital elements 
for an orbit around the Moon (see 
module description of ORBGEN) 

= 5, use DTAPRE to obtain spacecraft 
position (accesses a standard DODS 
EPHEM data set on FORTRAN unit 
NORBl; must be Moon-centered 
ephemeris) (see Reference 3) 

1*4 0 Moon ephemeris indicator: 

= 0, do not obtain Moon position. This 
. option must not be chosen if 

CBFLAG > 2, or if a Moon-centered 
ephemeris file is used for GETHDR/ 
GETVCT, or if ORBGEN is used 
for a spacecraft orbit around the 
Moon (ISPC = 4), or if DTAPRE 
is used with a Moon-centered 
EPHEM tape (ISPC = 5) 

= 1, use ORBGEN to generate orbit of 
Moon using orbital elements for 
an orbit around the Earth 
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Name Type Default 


Description 


IMOON 1*4 0 =2, use SUNRD to obtain Moon position 

(Cont'd) (accesses a direct access data set 

on FORTRAN unit 14; either equi- 
nox of date or equinox of 1950. 0, 
depending on data set accessed) 

(see Appendix B) 

= 3, use RJPLT to obtain Moon position 
(accesses a JPL ephemeris data 
set on FORTRAN unit NRJPLT; 
equinox of date) (see Reference 5) 

NOTE: The orbital elements and epoch 
time for the spacecraft must 
be specified if ISPC = 1 or 4. 
Orbital elements and epoch time 
for the Moon must be specified 
if IMOON = 1. 


NRJPLT 

1*4 

28 

FORTRAN unit number for JPL ephem- 
eris data set, read by RJPLT 

NORBl 

1*4 

30 

FORTRAN unit number for DODS 
EPHEM data set, read by DTAPRE 

NGTDS 

1*4 

29 

FORTRAN unit number for GTDS 
ORBIT file, read by GETHDR/GETVCT 
(only if a sequential file is used, 

LEVEL = 0; see below) 

LEVEL 

1*4 

0 

Level number for GETHDR/GETVCT: 


= 0, use sequential orbit file, FORTRAN 
imit number NGTDS 
> 0, level number on direct-access 
orbit file, FORTRAN unit number 
31 


TORBIT(6) 

R*4 

6*0. 

Epoch time of orbital elements for 
■ spacecraft (year, month, day, hour. 




minutes, seconds) 

TMOON(6) 

R*4 

6*0. 

Epoch time of orbital elements for 
Moon (year, month, day, hour, minutes, 
seconds) 
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Uncertainty parameters for OASYS data weighting are as follows: 


Name 

Type 

Default 

Description 

ERRBET 

R*4 

0.0 

Uncertainty in Sun angle (degrees) 

ERRGAM 

R*4 

0.0 

Uncertainty in sensor mounting angle 
(degrees) 

ERRA 

R*4 

0.0 

Uncertainty in rotation angle of space- 
craft from Sun crossing to horizon 
crossing (degrees) 

ERRAD 

R*4 

0.0 

Uncertainty in angular radius of the 
central body (degrees) 

ERRTIM 

R*4 

0.0 

Uncertainty in absolute time of horizon 
crossing, or uncertainty in ephemeris 
information (seconds) 

IDEBUG 

1*4 

0 

Debug indicator for subroutines UNCERT, 
UNCDBL, and UNCDH: 

= 0, no debug printout 
> 0, generate debug printout on 

FORTRAN data set reference num- 
ber IDEBUG 


NOTE: If it is desired to utilize the data weighting and uncertainty options of 
the system, the following parameters must be specified: 

• The NAMELIST parameter lAPIOR must be 4 or 8. The 
NAMELIST parameters APRA and APDEC must be specified. 

(See OPMANl Control.) 

• One or more of the uncertainty parameters (ERRBET, ERRGAM, 
ERRA, ERRAD, and ERRTIM) must have a nonzero value. 
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Bias parametei’s for OASYS are as follows: 


Name 

Type 

Default 

Description 

lOBLAT 

1*4 

1 

Kartli radius model indicator: 

= 0, use spherical Earth model, with 
radius = (RADE + IIT) 

= 1, use oblate Earth model, with at- 
mosphere layer at height HT above 
Earth surface (the value of RADE 
is ignored) 

HT 

R*4 

0.0 

Height of atmosphere layer above Earth 
surface (Idlometers) 

ABIASl 

R*4 

0.0 

Bias ;iddcd to the central body in- 
triggering (degrees) 

ABIAS2 

R*4 

0.0 

Bias added to central body out- 
triggering (degrees) 

BIASRE 

R*4 

0.0 

Bias added to angular radius of Earth 
(dfgrees) 

BIASRM 

R*4 

0.0 

Bias :idded to angular radius of Moon 
(degrees) 

RADE 

R*4 

6378. IG 

Effective optical aspect radius of Earth 
(kilometers) 

RADM 

R*4 

1738.0 

Effective optical aspect i^adius of Moon 
(kilometers) 

DANGE 

R*4 

89.7 

Dark :uigle of Earth (degrees) 

DANGM 

R*4 

89.7 

!!ar!c angle of Moon (degrees) 

BGAMMA 

R*4 

0.0 

Bias on sensor mounting angle (degrees) 

OB L ATI 

R*4 

0.0033328 

Earth obhitcness parameter (see note 
below ) 

OBLAT2 

R*4 

0.0 

Earth oblatcness parameter for asym- 
metric I kirtli (see note below) 


The effective Earth radius, R, 
in kilometers, is computed as 
U ^ (k'lTS. IG (1.0 - OBLYTl 
* sin“ 0 + OBLAT2 * sin 0) 
-» HT 

where 0 = geocentric latitude 
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OASYS printout control parameters are as follows: 

Name Type Default Description ■ 

IPLOT 1*4 5 Lineprinter plot level: 

= 0, no lineprinter plots 
s 1, plots of alpha versus frame number 
and delta versus frame number 
(including rejected points) 

^ 2, plots of alpha versus frame number 
and delta versus frame number 
(omitting rejected points) 
i .3, plots of arc length uncertainty 
versus frame number 

- 4, plots of nadir angle versus frame 

number 

- 5, plots of dihedral angle versus 

frame number 

DEBUG 1*4 8 Printout level indicator: 

= 0, no printout 

- 1, print error messages 
2, print averages 

Is print input parameters 
«• 4, print single frame results 
.■), print pointer frames 

- (i, i)rint input telemetry frames 

- 7, print Sun vectors 

s, print ephemeris vectors 
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OASYS control parameters 

are as 

follows: 

Name 

Type 

Default 

Description 

lAPIOR 

1*4 

0 

A priori attitude indicator: 

= 0, no a priori attitude available 
= 4, use a priori attitude only if deter- 
ministic logic fails 

= 8, use a priori attitude first to resolve 
all ambiguities 

APRA 

R*4 

0.0 

A priori right ascension (G. I.) (degrees) 

APDEC 

R*4 

0.0 

A priori declination (G. I. ) (degrees) 

EPSILN 

R*4 

50.0 

Sensor cutoff angle (degrees) 

EPS 

R*4 

1.4 

Angular field of view of sensor (degrees) 

ATTOL 

R*4 

360.0 

Error limit on a priori attitude (degrees) 

TRMCHK 

1*4 

0 

Terminator rejection flag: 

= 0, no restriction 

= 4, reject data if terminator is visible 

ITMAX 

1*4 

10 

Maximu m number of iterations allowed 
to resolve ambiguities in SPINA V (see 
Reference 8) 

SPNSIG 

R*4 

3. 

Rejection tolerance. An attitude is 
rejected if its deviation from the block 
average is greater than SPNSIG*STDV, 
where STDV is the computed standard 
deviation 

SPNTOL 

R*4 

20. 

Rejection tolerance. An attitude is 
rejected if its deviation from the block 
average is greater than SPNTOL (de- 
grees) 

SRLOW 

R*4 

0. 

Lower bound on spin rate (degrees per 
second) 

SRHIGH 

R*4 

1.E6 

Upper bound on spin rate (degrees per 
second) 

SUN LOW 

R*4 

0. 

Lower bound on Sun angle (degrees) 

SUNHI 

R*4 

180. 

Upper bound on Sun angle (degrees) 
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Name 

Type 

Default 

Description 

CBFLAG 

1*4 

2 

Central body indicator: 

= 1-2, consider Earth only as central 
body 

= 3-6, consider both Earth and Moon 
as central bodies 

= 7-10, consider Moon only as central 
body 

DELPHI 

R*4 

0.5 

Terminator search increment, used by 
predictor. A pair of terminator 
crossings may be overlooked if their 
angular separation is less than DELPHI 
(degrees) 

PHITOL 

R*4 

0.01 

Terminator angle tolerance. Used to 
terminate half-interval search in predic- 
tor (degrees) 

DIHTOL 

R*4 

10. 

Maximum expected error in degrees in 
dihedral angle computed from double 
horizon scan midtime. Used to resolve 
ambiguities in double horizon width 
method 

IDISK 

1*4 

0 

OABIAS data archive control parameter: 
^ 0, do not write archive data set 
^ 0, write archive data set on FORTRAN 
unit number IDISKi See Section 7.4 
for a description of this data set 

BDYCHK 

1*4 

0 

Double horizon crossing flag: 

= 0, use double horizon methods only 
if central body is fully sunlit 
= 1, use double horizon methods only 
if a scan at the a priori attitude 
does not cross the terminator 
= 2, use double horizon methods regard- 
less of lighting conditions on cen- 
tral body 

MIN 

1*4 

1 

Single horizon-in flag: 


= 0, do not process with single horizon- 
in method 


= 1, process with single horizoa-ia 
method 
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Name 

Type 

Default 

Description 

MOUT 

1*4 

1 

Single horizon-out flag: 

= 0, do not process with single horizon- 
out method 

= 1, process with single horizon-out 
method 

MDBL 

1*4 

1 

Double horizon width flag: 

= 0, do not process with double horizon 
width method 

= 1, process with double horizon width 
method 

MDIH 

1*4 

1 

Double horizon dihedral flag: 

= 0, do not process with double horizon 
dihedral method 

= 1, process with double horizon dihedral 
method 

SUNSIG 

R*4. 

3. 

Rejection tolerance. A Sun angle is 
rejected if it differs from its smoothed 
value by more than SUNSIG*STDV, 
where STDV = the computed standard 
deviation 

SUNTOL 

R*4 

20. 

Rejection tolerance. A Sun angle is 

rejected if it differs from its smoothed 

value by more than SUNTOL (degrees) 

1 

ISNPRO 

1*4 

0 

Sun angle processing indicator: 

= 0, normal Sun angle processing (linear 
fit) 

= 1, do not smooth Sun angles 

SUNIN 

R*4 

0.0 

Sun angle input; 

= 0. 0, no effect 

4 0.0, use the value SUNIN for every 
Sun angle; ignore the Sun angles 
in the telemetrj% and ignore die 
value of ISNPRO (degrees) 


NOTE: The following parameters in NAME LIST OPMANl are not applicable for 
MSAD/OABIAS ; see Reference 8 for a discussion of these parameters: 


SENANG STIMEl SRSIG WTNAD 

APPER FTIilEO SRTOL WTDIII 

MIDTBI FTIiNIEl ISRPRO WTNAD2 

TWTDTH IGCORE WTSUN WTD1H2 

STIIvIEO 
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7.3.7 NAMELIST BIASNL 

NAMELIST BIASNL is read by subroutine INITIL and contains parameters 
specifying processing options for the OABIAS Subsystem. 

Parameters defining the initial state are as follows: 


Name Tvpe Default Description 


XO(l) 

R*8 

0.0 

Initial estimate of ri^t ascension (G. I. ) 
(degrees) 

XO(2) 

R*8 

0.0 

Initial estimate of declination (G. I. ) 
(degrees) 

XO(3) 

R*8 

0.0 

Initial estimate of phase (degrees) 

= 0. 0, program will compute the initial 
phase 

XO(4) 

R*8 

0.0 

Initial estimate of bias in horizon 
scanner mounting angle from spin axis 
(degrees) 

XO(5) 

R*8 

0.0 

Initial estimate of borizon scanner 
azimuth relative to the Sun sensor on 
central body in-triggering (degrees) 

XO(6) 

R*8 

0.0 

Initial estimate of horizon scanner 
azimuth relative to the Sun sensor on 
central body out-triggering (degrees) 

XO(7) 

R*8 

0.0 

Initial estimate of bias on the angular 
radius of the central body (degrees) 

XO(8) 

R*8 

0.0 

Initial estimate of bias in Sun ar^le 
(degrees) 

XO(9) 

R*8 

0.0 

Initial estimate of spin rate (revolutions 
per minute); 

= 0. 0, compute initial estimate using 

average of spin rates in INDATA 

XO(IO) 

R*8 

0.0 

Initial estimate of Sun sensor plane tilt 
(degrees) 

XO(ll) 

R*8 

0.0 

Initial estimate of horizon sensor plane 
tilt (degrees) 

XO(12) 

R*8 

0.0 

Initial estimate of orbit time adjustment 
(seconds) 
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Name 

Type 

Default 

Description 

PO(12) 

R*8 

12*0.0 

Estimate of the error in the initial state 
(units same as XO) 

ATTOPT 

1*4 

2 

Initial attitude options: 

= 1, use OASYS block average results 


for both attitude and uncertainties; 
ignore XO(l), XO(2), PO(l), PO(2) 
= 2 , use XO(l) and XO(2) for initial 
attitude; use PO(l) and PO(2) for 
initial uncertainties 
= 3, use state vector and uncertainties 
from previous OABIAS run; ignore 
XO and PO entirely 
= 4, use state vector from previous 
OABIAS run, and input values of 
PO; ignore XO entirely 
= 5, use OASYS block average results 
for attitude, and input values for 
PO; ignore XO(l) and XO(2) 
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Parameters specifying OABIAS models are as follows: 


Name 

Type 

Default 

Description 

MODEL(IO) 

L*4 

10*F 

Observation model selector. If true^ 
the selected model is processed 

(1) 



Sun angle 

(2) 



Sun time 

(3) 



Nadir vector projection 

(4) 



Horizon time 

(5) 



Single horizon dihedral angle 

(6) 



Double horizon scan width 

(7) 



Small target nadir angle 

(8) 



Scan midtime dihedral angle 

(9) 



Unused 

(10) 



Unused 

RESMOD(IO) 

L*4 

10 *F 

Residual calculation selector. If true, 
the residual for that model is computed, 
RESMOD(i) is defined in the same way 
as MODEL(i) 
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Parameters defining processing options are as follows: 


Name 


RE J OPT . 


Type 

P4 


Default 


lUPDAT 


1*4 


NUMITR 


1*4 


CONVRG 


R*8 


11*1. D-5 


REDUNT(IO) 1*4 


5,9*0 


INTOUT 


lOBLAT 


ITER 


1*4 


1*4 


1*4 


Description 


Data rejection option: 

= 0, process all data 
= 1, reject an observation if it was re- 
jected by the corresponding model 
in OASYS 

State vector updating sequence: 

= 0, update after every model process- 
ing 

= N, update at every Nth frame 

Iteration controller; 

= 0, no iteration 

= N, allow up to N iterations on each 
state update 

Array of convergence criteria for itera- 
tion loop. If A < CONVRG, the 
iteration loop is completed (units the 
same as XO) 

In a nongraphic mode, if REDUNT(l) ^ 0, 
another set of parameters for NAMELIST 
BIASNL will be read from Idle FORTRAN 
unit number specified by REDUNT(l). 

If REDUNT(l) = 0, OABIAS will return 
after completing the processing re- 
quested in this NAMELIST. In a graph- 
ics mode, REDUNT has no effect. 
REDUNT(2) throu^ REDUNT (10) are 
not applicable in MSAD/ OABIAS 

Results printout frequency controller; 

= 0, print after'each model processing 
= N, print after each Nth frame 

Oblateness model indicator; 

= 0, use spherical Earth model 
= 1, use oblate ilarth model 

Maximum number of iterations for 
linearity fix: 

= 1, do not use linearity fix 
= N, maximum of N iterations 


iM 
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Type Default Description 

1*4 0 Spin rate processing option: 

= 0, use block average spin rate for 
processing every frame 
= 1, use spin rate from each telemetry 
frame for processing that frame 
ISPINR does not apply if PO(9) 5 ^ 0.0 
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Parameters controlling OABIAS printout are as follows; 


Name 

Type 

1*4 

Default 

Description 

ERUNIT 

19 

Unit number for printout of error mes- 
sages and diagnostic output 

IDUMPL 

1*4 

0 

Diagnostic output level controller; 

= 0, print only summary of errors 
== 4, print error messages and summary 

IPLOT (12) 

1*4 

12*3 

Plotting option for printer plots of state 
vector and state uncertainties: 


= 0, do not plot state or uncertainty 
= 1, plot state only 
= 2, plot uncertainty only 
= 3, plot both state and uncertainty 

There is one indicator for each state 
component as listed for XO. Note that 
if PO(I) = 0. 0 for any state component, 
plots will not be generated for that state 
component, regardless of the value of 
IPLOT 


TDIAGN{32) 1*4 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

(6) 
< 7 ) 
( 8 ) 

(9) 

(10) 
( 11 ) 


32*0 Diagnostic printout flag for each sub~ 

routine: 

= 0, print nothing 
s 1, print header and trailer 
^ 2, input variables 
2 4, output variables 
^ 6, internal variables before and after 
they are passed via an internal 
subroutine call 
2 8, other internal variables 

AMATRX 

A PARTS 

BIASER 

DIAFUN 

DIAMOD 

FRAPRO 

INITIL 

LCOMP 

LNFUN 

LNMOD 

LPARTS 
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Name Type Default 


Description 


(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27)-(32) 

IPRINT 1*4 20 


ILEVEL 1*4 4 


LRFUN 

LRMOD 

OABIAS 

PRINT 

PSIPHA 

RECURS 

ROTATE 

SANE UN 

SANMOD 

STMFUN 

STMMOD 

DHFUN 

DHMOD 

SCBFUN 

SCBMOD 

Unused 

FORTRAN unit number for results 
printout. Results are printed on units 
IPRINT, IPRINT + 1, IPRINT + 2, 
IPRINT + 3, IPRINT + 4, and IPRINT + 5 

Results printout level controller: 

= 0, print nothing 

= 2, print time, state vector, a , 6 , 
frame number, observation type 
= 4, print above plus covariance matrix 
= 6, print above plus gain matrix 
= 8, print above plus matrix of parti als, 
weights, residuals 

= 10, print above plus observation and 
calculated observation 
= 12, print above plus inertial vectors 
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BIASNL 
Data Weighting- 


Parameters defining observation weights are as follows: 


Name 

Type 

Default 

Description 

SANGRN 

RM 

l.E-4 

Sun angle granularity (degrees) 

STMGRN 

R*4 

l.E-4 

Sun time granularity (seconds) 

OAGRAN 

R*4 

l.E-4 

Timing granularity of the optical aspect 
(OA) telescope (seconds) 

RHOGRN 

R*4 

1.0 

Granularity in central body angular 
radius for model 7 (degrees) 
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Parameters defining OABIAS calculation tolerances are as follows: 


Name 

Type 

Default 

Description 

TOLDEN 

R*4 

l.E-4 

Mimmum value allowed for S 2 R 2 “ S2Rt 
in LCOMP subroutine 

TOLDET 

R*4 

l.E-4 

Determinant singularity tolerance in 
LPARTS subroutine 

BELT AT 

R*4 

10.0 

Time increment used for computing 
spacecraft velocity by numerical differ- 
entiation (seconds) 

ROTLIM 

R=*=4 

20.0 

Declination tolerance for coordinate 
transformation (degrees). A rotation 
is performed if the absolute value of 
the initial declination is < ROTLIM 
or >90 - ROTLIM 

UNTOL(12) 

R*8 

12*1. D-5 

Tolerances for linearity fix (units the 
same as for XO) 

GAMDEG 

R*8 

— 

Not used 

IDISK 

1*4 

— 

Not used 
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7.3.S Description of a Sample NA M KfJ.ST 

This subsection describes a sample set of input cards. The input NAMELIST 
parameters together with the default values arc designed to first simulate a 
spacecraft with an infrared (IR) sensor in an elliptic orbit. The simulated data 
span covers G minutes about 1 hour after perigee. The simulated data is first 
processed by OASYS and then by OABIAS. This sample problem provides the 
printout and graphs which are described in Sections 7. 12 and 7. 14. 


The sample NAME LIST s together with the default values are given in Sec- 
tion 7.3. 8. 1. A description of tlie effects of each NAMELIST is given sepa- 
rately in the following subsections. 


7. 3. 8. 1 Sample NAMELIST 
Card Input 


Column 2 
&CONTRL 


&END 

&JIAIN 


NFRAME - 25 


Default Values Assumed 


IFTABL = 96 
IFTUBE = 23 
IFTPRT = 6 
IRDART = 0 
IRDXST = 0 ' 
GSATID = 
GRUNID = 


ISIM = 1 

IREWND = 1 
OMEGAl = 100., 0. 
IRECl = 1 
IREC2 = 99999 
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Card Input 


IDISK = 51 


NFRAM2 = 25 
lEPHEM = 0 
ITERM = 0 
ICB = 0 


INOISE = 2 

STDV = 4* 0. 005, 
2 * 0.0 


SQUANT = 0.5 


&END 

&LIST ALPHA = 351. 

DELTA = -20. 
lATAPE = 1 


Default Values Assumed 

lADD = 0 
ISKIP = 1 
IHEADR = 3 
ICBSEL = 1, 1 , 0 , 1 
ICNSEL = 8*1 

NFRAMl = 1 


TADJ = 0.0 


BIAS = 6 *0.0 
PROB = 6 * 0. 0 
IRAND = 1^3456789 
BGAM =0.0 


ABETAl = 0.0 
ABETA2 = 0.0 
OMEGl = 0.0 
OMEG2 = 0.0 
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Card Input 

TO = 1974., 1., 4., 

23., 4., 0. 

DELTAT = 15. 

A = 24548., 0. 

E = 0. 73264, 0. 

EYE = 28.3, 0. 

WO = 180. , 0. 
RANODE = 260. , 0. 

ORBITE = 1974. , 1. , 

4. , 22. , 

51., 26. 


ISPC = 1 


ISKIP = 1 
IFLAG = 2 
SIGMA = 86. 
STEP = 0.0 


lOBLAT = 0 


ICB = 1 


THETAC = 0.0 
DANGE = -90. 


Default Values Assumed 


EMO = 0. , 0. 


ISUN = 1 

IMOON = 0 
IPRINT = 0 
IREPOR = 0 


BIASRE = 0. 0 
RADE = 6378.16 

HT = 0. 0 


EPS = 0.0 
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&END 

&OPMAN1 


Card Input Default Values Assumed 

DELPHI = 0.5 
PHITOL = 0, 01 

ERRBET = 0.25 
ERRGAM = 0. 1 
ERRA = 0. 1 
ERRAD = 0.5 
ERRTIM = 30. 

IDEBUG = 0 
TOL = 0.0001 
MAXIT = 10 

ISUN = 1 

ISPC = 1 

IMOON = 0 

TORBIT = 1974. , 1. , 

22., 51., 

26. 

A = 24548., 0. 

E = 0.73264, 0. 

EYE =28.3, 0. 

EMO = 0. 0, 0. 

W.O = 180. , 0. 

RANODE = 260. , 0. 

ERRBET = 0.25 
ERRGAM = 0. 1 

ERRA = 0. 1 ■ - 

ERRAD = 0.5 
ERRTIM = 30. 
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Card Input 


lOBLAT = 0 


DANCE = -90. 0 


lAPIOR = 8 
APRA = 350. 
APDEC = -19. 
EPSILN = 0.0 


CBFLAG = 1 


Default Values Assumed 


IDEBUG = 0 


HT = 0. 0 
ABIASl = 0.0 
ABIAS2 = 0. 0 
BIASRE =0.0 
RADE = 6378.16 


BGAMMA = 0.0 
IPLOT = 5 
DEBUG = 8 


EPS = 1.4 
ATTOL = 360. 
TRMCHK = 0 
ITMAX = 10 
SPNSIG = 3. 
SPNTOL = 20. 
SRLOW = 0. 
SRHIGH = 1. E6 
SUNLOW = 0. 
SDNHI = 180. 

DELPHI = 0. 5 
PHITOL = 0.01 
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Card Input 


ISNPRO = 1 


&END 

&BIASNL 


PO = 3 * 100. , 

9 * 10. , 

ATTOPT = 5 

MODEL = 6 ^ T, F, 
T, 2 * F 


RESMOD = 6 * T 
F, T 
2 * F 

REJOPT = 1 


REDUNT(l) = 0 


ERUNIT = 12 
- IDUMPL = 4 


Default Values Assumed 

DIHTOL = 10. 

IDISK = 0 
MIN = 1 
MOUT = 1 
MDBL = 1 
MDIH = 1 
SUNSIG = 3. 

SUNTOL = 20. 


SUNBSr = 0.0 


XO = 12 ♦ 0.0 


lUPDAT = 0 
NUMITR = 0 

INTOUT = 0 
lOBLAT = 0 
ITER = 1 
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Card Input Default Values Assumed 

IPLOT = 12 3 

IDIAGN = 32 * 0 

IPRINT = 6 
ILEVEL = 12 
SANGRN = 30. 

STMGRN = 0. 1 
OAGRAN = 0.5 

TOLDEN = l.E-4 
TOLDET = l.E-4 
ROTLIM = 20. 

&END 

7. 3. 8. 2 NAMELIST CONTRL 

The GESS nonresident tables will be read on FCiRTRAN unit number 96, and 
the graphic device specified on FORTRAN unit number 23 will be used (default; 
IFTABL = 96; default: IFTUBE = 23). GESS output will be printed on 
FORTRAN unit number 6 (default: IFTPRT = 6), 

I 

7. 3 . 8 . 3 NA ME LIST MAIN 

Twenty -five frames of data (N FRAME = 25) will be simulated (default: 

ISIM = 1) using the constant spin rate of 100 degrees per second (default: 
OMEGAl = 100. , 0. ). The data will be written on a device specified by 
FORTRAN unit number 51 (IDISK = 51). The device will be rewound prior to 
being written on (default: IREWND = 1). 

The ephemeris data, terminator flags, and central body flags on the OABLAS 
data set will not be used in the OASYS and OABIAS runs (lEPHEM = 0, 

ITERM = 0, ICB = 0). These options will let OASYS obtain the ephemeris data 
and identify the central bod 3 ’- and terminator crossings. The data will be 
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processed once with noise (IN OISE = 2) generated the standard deviations 
specified by STDV(6). No biases will be added to the simulated data (default: 
BIAS = 6 * 0. 0; default: BGAM = 0.0; default: TADJ = 0. 0) , but the Sun 
angles will be quantized (SQUANT = 0. 5). 

7. 3. 8. 4 NAMELIST nST 

This NAMELIST contains parameters specifying the simulation conditions. The 
data will be simulated every 15 seconds (DELTAT = 15.), starting at 2304 hours, 
January 4, 1974 (TO = 1974. , 1. , 4. , 23. , 4. , 0. ), based on the attitude 
(ALPHA = 351. , DELTA = 20. ) read from cards (lATAPE = 1). 

The Earth is the only central body considered here (ICB = 1), and the Moon po- 
sition is not required (default: IMOON = 0). The spacecraft orbit will be gen- 
erated by subroutine ORBGEN (ISPC = 1) using the orbit parameters A, E, 

EYE, EMO, WO, RfYNODE, and ORBITE. The Sun position will be obtained 
by subroutine SUNl (default: ISUN = 1). 

It is assumed that the horizon sensor has a mounting angle of 86. 0 degrees 
(SIGMA = 86. 0) and an angular field of view of 0. 0 degree and that there is no 
cutoff angle (THETAC = 0. 0). Because the sensor mounting angle is constant, 
STEP must be set to zero (STEP = 0. 0). No bias on the angular radius of 
Earth is considered (default: BIASRE = 0. 0). Note that all biases except 
BIASRE are in NAMELIST MAIN. This is because the bias on the angular 
radius, unlike the other biases, cannot be simply added to unbiased data. 

A spherical Earth model is used (lOBLAT = 0). In this example an IR sensor is 
assumed so that the terminator has no influence. In order to ensure both horizon 
in- and out-triggerings, the dark angle of Earth must be set to -90. 0 degrees 
(DANGE = -90). The IR sensor is assumed to trigger at the surface of Earth; 
i.e. , the effective height of atmosphere layer is zero (default: HT = 0. 0). 

The expected qualitj’- of the attitude solutions obtained by OASYS will be printed 
out in terms of the uncertainty in the attitude caused by uncertainties in Sun 
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angle, sensor mounting angle, rotation angle, angular radius of Earth, and 
absolute time of horizon crossing (ERRBET = 0.25, ERRGAM = 0. 1, 

ERRA = 0. 1, ERRAD = 0. 5, ERRTIM = 30. ). 

7. 3. 8. 5 NAMELIST 0PM AN 1 

This NAMELIST contains parameters specifying the processing options for the 
OASYS Subsystem. 

The spacecraft orbit and the Sun position will be obtained in the same way as 
in the simulation (default: ISUN = 1, ISPC = 1; default: IMOON = 0). The or- 
bit parameters. A, E, EYE, EMU, WU, RAN ODE, and TORBIT, should be 
the same as those in NAMELIST LIST. Note that the epoch time here is de- 
noted by TORBIT while ORBITE is used in LIST. 

The computation of the uncertainty in the determined attitude will be based on 
the uncertainty in the Sun angle, sensor mounting at^le, rotation angle, angular 
radius of the Earth, and absolute time of the horizon crossing (ERRBET = 0.25, 
ERRGAM = 0. 1, ERRA = 0. 1, ERRAD = 0. 5, ERRTIM = 30. ). 

As in the simulation, the spherical Earth model with zero effective atmosphere 
layer and a dark angle of -90. 0 degrees will be used (lOBLAT 0; default: 

HT = 0. 0; DANGE = -90, ), No bias adjustment will be made in determining 
the attitude (default: ABIAS2 = 0.0; default: BIASRE = 0.0, default: 

BGAMMA = 0. 0). 

As in the simulation, only the Earth will be considered as a central body 
(CBFLAG = 1). The Sun angles will not be smoothed (ISNPRO = 1). Note that 
the Sun angles in this example problem are quantized as described in LIST. 

All four attitude determination methods, the single horizon-in and horizon-out, 
double horizon width, and double horizon dihedral method, will be processed 
to determine the attitude (default; MIN = 1; default; MOUT = 1; default: 

MDBL = 1; default: MDIH = 1). The a priori attitude is used first to resolve 
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all ambiguities (APRA = 350. , APDKC = 19,, lAPIOR = 8). All solutions will 
be included prior to the block averaging (default: ATTOL = 360. ). An attitude 
will be rejected if it deviates from the block average by more than 20 degrees 
or by more than 3 times the computed standard deviation (default: SPNTOL = 20 
default: SPNSIG = 3.). The attitudes, uncertainties, nadir angles, and dihe- 
dral angles obtained from the four attitude methods will be plotted by the line- 
printer (IPLOT = 5). 

7. 3. 8. 6 NAMELIST BIASNL 

This NAMELIST contains parameters specifying processing options for the 
OABIAS Subsystem. In this sample problem, the user is assumed to have 
no a priori estimate of the attitude and biases. The best choice, then, is 
ATTOPT = 5, which allows OABIAS to use the block average attitude from 
OASYS as the initial estimate. Zero bias will be used as the initial estimate 
(default: XO = 12 * 0. 0). The initial estimate of the uncertainties will be 
given by PO. All 12 elements of the state vector will be corrected. Note that 
PO(J) = 0. 0 restricts Jth component of the state vector to its initial value. 

All models but the small target nadir angle model will be processed, and the 
residuals for each model will be computed (MODEL = 6 * T, F, T, 2 * F, 
RESMOD = 6 * T, F, T, 2 * F). Generally, updatii^ the state vector after 
processing each observation (default: lUPDAT = 0) with no iteration (default: 
NUMITR = 0) gives the best results. Since the spherical Earth model is used 
in simulation, the use of that model is appropriate (default; lOBLAT = 0 ). 
Usually, the attitude estimated by OASYS is sufficiently close to the true state, 
unless unexpectedly large biases are present. The linearity' fix, then, is not 
necessary (default: ITER = 1). The data rejected in each model by the block 
averaging in OASYS will be rejected in processing of lie iu rresponding model 
(REJOPT = 1) in OABLAS. 
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The printout will include tlie state vector, uncertainties, partial derivatives 
vector, gain vectors, observation and weight, and the statistics of the resid- 
uals (ILEVEL = 12) on FORTRAN units numbered 6, 7, 8, 9, 10, 11, respec- 
tively (IPRINT = 6). Both the state vector and uncertainties will be plotted 
(default: IPLOT = 12 * 3) on FORTRAN unit number 6 (IPRINT = 6). Note 
that the residuals for a model will be plotted on FORTRAN unit number 6 
(IPRINT = 6) if the corresponding array element of RESMOD is true. A line 
of printout will occur for every observation, and a point will be plotted at every 
frame (default: INTOUT = 0), The error messages and error summary table 
will be printed (IDUMPL = 4) on FORTRAN unit number 12 (ERUNIT = 12). 

No diagnostic printout will be generated (default: IDIAGN = 32 * 0). 


The granularities required for each model are shown below: 


Model 


Granularity Required 


Sim angle model 

Sun time model 

Nadir vector projection model 

Horizon time model 

Single horizon dihedral angle model 

Double horizon scan width model 

Small target nadir angle model 

Scan midtime dihedral angle model 


SANGRN 

STMGRN 

OAGRAN 

OAGRAN 

STMGRN, OAGRAN 

OAGRAN 

RHOGRN 

STMGRN, OAGRAN 


Since all models but the small target nadir angle model will be processed, all 
granularities but RHOGRN are required (SANGRN = 30. , STMGRN = 0. 1, 
OAGRAN = 0. 5). 

/N /N 

When any two vectors of the horizon vector, ; spin axis vector, S ; and 
spacecraft position vector, R , are either nearly parallel or nearly antiparal- 
lei, a singularity in computing the partial derivatives of is encountered* 
Whenever the estimate and geometry result in a matrix whose determinant is 
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less than TOLDET, the processing of that observation is terminated (default: 
TOLDET = 1. E-4). TOLDET is required for the horizon time model, single 
horizon dihedral angle model, and double horizon scan with model. 

The G. I. coordinate system will be rotated if the absolute value of the initial 
estimate of the declination of the spin axis is less than 20 degrees or between 
70 and 90 degrees (default: ROTLIM = 20. ). 

The NAMELIST parameters that have no effect in this sample problem are 
CONVRG, RHOGRN, DELTAT, and LINTOL. If an iteration on each state up- 
date is needed (NUMITR ^ 0), CONVRG is required. If the linearity fix is at- 
tempted (ITER >1), LINTOL is required. If the correction of the in-track time 
error is required and the central body is the moon, DELTAT must be specified. 

7.4 OABIAS DATA SET 

The OABIAS data set is written by the Radio Astronomy Explorer-B (RAE-B) 
Attitude Determination System or written by the OASYS Subsystem running 
under the SSS-A, SMS-A, or IMP-J Attitude Determination Systems. In a sim- 
ulation mode, the simulator within MSAD/OABIAS also writes a data set in this 
format. For the AE spacecraft, the subroutine AECOPY copies the AE data 
set into this format. 

The data set consists of one or more blocks of data, with each block preceded 
by a pair of header records. Each block of data consists of one or more 
frames; each frame consists of two records, an in-triggering followed by an 
out-triggering. 

The data set is a FORTRAN A-format data set (i. e. , all parameters are read 
and written in an A format); it resides on FORTRAN unit IDISK, where IDISK 
is a parameter in NAMELIST MAIN (default value = 50). 
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The two 50-byte header records which precede each block of data are as 
follows: 

Record 1: 


Name 

Description 

Type 

IHFIAG 

Header flag = -999 

1*2 

TREF 

Current time of job running (local time) 

(year, month, day, hours, minutes, seconds) 

R*8 

TFIRST 

Time of first sample (UT; seconds from 
0 hours, September 1, 1957) 

R*8 

THE^ (1) 

Upper SAS zenith reference angle from 
+Z (degrees) 

R*4 

(2) 

Lower SAS zenith reference angle from 
+Z (degrees) 

R*4 

(3) 

PASl null zenith reference angle from 
+Z (degrees) 

R*4 

(4) 

PAS2 null zenith reference angle from 
+Z (degrees) 

R*4 

PAD (1) 

Padding (ignored) 

R*8 

(2) 

Padding (ignored) 

R*8 

Record 2: 



Name 

1 

Description 

Type 

NSAMP 

Number of samples = number of data records 
following this header 

1*2 

PHI^ (1) 

Upper SAS azimuth reference angle from 
+X (degrees) 

R*4 

(2) 

Lower SAS azimuth reference angle from 
+X (degrees) 

R*4 

(3) 

PASl null azimuth reference angle from 
+X (degrees) 

R*4 

(4) 

PAS2 null azimuth reference angle from 
+X (degrees) 

R*4 

^Applies to R\E' 

-B spacecraft only. 
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PRELIMIMARY DRAFT 

Description 

Type 

PHl^ (5) 

PASl Sun azimuth reference angle from 
+X (degrees) 

R*4 

(6) 

PAS2 Sun azimuth reference angle from 
+X (degrees) 

R*4 

RHOEB 

Bias on Earth angular radius (degrees) 

R*4 

RHOMB 

Bias on Moon angular radius (degrees) 

R*4 

PAD (1) 

Padding (ignored) 

R*8 

(2) 

Padding (ignored) 

R*8 

The 50-byte data records (always occurring in pairs) which constitute each 
frame on a block of data are as follows: 

Name 

Description 

Type 

INVERT^ 

Inversion flag: 
= 0, normal 
= 1, inverted 

1*2 

TTRIG 

Time of horizon triggering (UP: seconds from 
0 hours, September 1, 1957) 

R*8 

BETA 

Sun angle (degrees) (0-180) 

R*4 

GAM 

Sensor mounting angle (degrees) (0-180) 

R*4 

RPM 

Spin rate (revolutions per minute) 

R*4 

A 

Dihedral angle from Sun to horizon (degrees) 
(0-360) 

R*4 

PASPHI^ 

Phase angle used in computing A (degrees) 
(0-360) 

R*4 

RMINUS (3) 

Unit vector from spacecraft to center of cen- 
tral body, inertial coordinates 

R*4 

RMAG 

Distance from spacecraft to center of central 
body (Idlometers) 

R*4 

KCNFLG 

Configuration flag (see note below) 

I ^2 


Applies to RAE-B spaceci-aft only. 
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Name 


Description 


Type 


CBPAS Central body indicator: L*1 

= E, Eartli 
= M, Moon 

= S, Sun (central body vector ignored) 

= U, unlcnown (central body vector ignored) 

HORTER Horizon terminator indicator: L*1 

= H, horizon 
= T, terminator 
= U, unknown 


NOTE: For the RAE-B spacecraft, the configuration flag is defined as follows: 


KCNFLG 

Sun Angle 

Sun Sensor Used for 
Sun Triggering 

Horizon Sensor 
Used 

= 1 

<90 

SAS 

PASl 

= 2 

<90 

SAS 

PAS2 

= 3 

<90 

PASl 

PASl 

= 4 

<90 

PAS2 

PAS2 

= 5 

>90 

SAS 

PASl 

= 6 

>90 

SAS 

PAS2 

= 7 

>90 

PASl 

PASl 

= 8 

>90 

PAS2 

PAS2 


For the AE spacecraft, the configuration flag is defined as follows : 


= 1, wheel horizon sensor 1, gimbal Sun sensor • 

= 2, wheel horizon sensor 2, gimbal Sun sensor 
= 3, body horizon sensor 1 or 2, digital Sun sensor 

For all other spacecraft, the configuration flag is fixed at 1. 

7. 5 AE DATA SET 

The AE input data set is written by the AE-C Attitude Determination System. 
This data set can be accessed only through the use of subroutine RDBUFF from 
the AE libraries. See Reference 14 for a description of this data set. The 
AE input data set is read from FORTRAN unit number lAEUN, a parameter in 
NAMELIST MAIN. 
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Not all of the data brought in to subroutine AECOPY is used in OABIAS proc- 
essing. Since some of the data which is not used is of interest in modeling the 
performance of the bolometer for the AE wheel mounted sensors, the option 
exits to write wheel horizon sensor data on FORTRAN unit number IBOLOU. 
IBOLOU is a parameter in NAMELIST MAIN, The AE bolometer data set con- 
sists of one or more frames of data with each frame in the following format; 


Name 

Description 


WHSTIM 

Wheel horizon sensor tri^ering time for 
each frame (UT: seconds from 0 hour, 
September 1, 1957) 

R*8 

EWl 

Earth width times for wheel horizon sensor 1 
(seconds) 

R*4 

EW2 

Earth width times for wheel horizon sensor 2 
(seconds) 

R*4 

SI 

Split-to-index pulse times (seconds) 

R*4 

WHSSP 

Wheel spin period for each frame (seconds) 

R*4 

WHSFLG 

Wheel horizon sensor flag; 
= 0, (EBCDIC) good frame 
^ 0, bad frame 

L*1 


7. 6 JPL LUNAR AND SOLAR EPHEMERIS FILE 

The JPL ephemeris file contains lunar and solar ephemeris information ac- 
cessed by subroutine RJPLT. The file is required if the parameter ISUN = 3 
or IMOON = 3 in NAMELIST LIST or OPMANl. (See Sections 7.3. 5 and 
7.3.6.) The FORTRAN unit number is specified by the parameter NRJPLT in 
NAMELIST LIST or OPMANl (default value = 28), 

See Reference 5 for a description of this data set. 

7.7 SL’NRD LUNAR AND SOLAR EPHEMERIS FILE 

The SUNRD ephemeris file contains lunar and solar ephemeris data accessed 
by subroutine SUNRD. The file is required if the parameter ISUN = 2 or 
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IMOON = 2 in NAMELIST LIST or 0PM AN 1. (See Sections 7. 3. 5 and 7. 3. 6. ) 
The FORTRAN unit number for the file is 14. 

See Appendix B in Volume 1 for a description of this data set. 

7.8 GTDS SPACECRAFT EPHEMERIS FILE 

The GTDS ephemeris file contains spacecraft ephemeris data accessed by sub- 
routine GETHDR/GETVCT. The file is required if the parameter ISPC = 3 in 
NAMELIST LIST or OPMANl. (See Sections 7. 3. 5 and 7.3. 6. ) If the param- 
eter LEVEL = 0, then the FORTRAN unit number is specified by the parameter 
NGTDS (default value = 29); if LEVEL ^ 0, the FORTRAN unit number is 31. 

See Reference 4 for a description of this data set. 

7.9 DODS SPACECRAFT EPHEMERIS FILE 

The DODS ephemeris file contains spacecraft ephemeris data accessed by 
DTAPRE (an entry point of ROITAP). The file is required if ISPC = 2 or 5 
in NAMELIST LIST or OPMANl. (See Sections 7.3.5 and 7. 3. 6. ) The 
FORTRAN unit number is specified by the parameter NORBl (default value = 30). 

See Reference 3 for a description of this data set. 

I 

7.10 SPACECRAFT ATTITUDE TAPE 

The spacecraft attitude tape is a standard Multisatellite Attitude Prediction 
(MSAP) attitude data set which is used by the simulator to simulate a space- 
craft with a slowly varying attitude. This data set is required if the parameter 
lATAPE = 0 in NAMELIST LIST. The FORTRAN unit taimber for the attitude 
tape is 12. 

See Reference 20 for a description of this data set. 


7-65 

PRELIMINARY DRAFT 



PRELIMINARY DRAFT 

7. 11 GESS NONRESIDENT TABLES 

The GESS nonresident tables reside on a partitioned data set and contain the 
display tables which describe each GESS display in the system. The tables 
are accessed one at a time by the GESS Executive as required. The FORTRAN 
unit number is specified by parameter IFTABL in NAMELIST CONTRL (default 
value = 96). 

See Reference 13 for a description of this data set. 

7. 12 PRINTED OUTPUT 

Printed output is generated by the GESS Executive, the OADRIV control mod- 
ule, the Optical Aspect Determination Prediction (ODAP) Subsystem, the 
OASYS Subsystem, and the OABIAS Subsystem. 

Printed output may be generated on FORTRAN unit numbers 6, 9, 10, 94, and 
95, plus the unit numbers specified by the parameters IFTPRT (in NAMELIST 
CONTRL), ERUNIT (in NAMELIST BIASNL), and IPRINT through IPRINT + 5 
(in NAMELIST BIASNL). Thus if IFTPRT = ERUNIT = IPRINT = 6, then 
units 6, 7, 8, 9, 10, 11, 94, and 95 may be used. The following data definition 
(DD) card is required for each unit: 

//FT XX FOOl DD SYSOUT = A, DOB = (RECFM = VBA, 

// LRECL = 137, BLKSIZE = 7265, BUFNO = 1) 

The block size for each data set should be chosen to optimize input/output 
(I/O) time and core requirements based on the amount of printout expected for 
each unit. A DUMMY DD card should be provided for any unit number from 
which printout is not desired. 

7, 12. 1 Printed Output From the GESS Executive 

Printed output from the GESS Executive consists of a listing of the parameters 
in NAMELIST CONTRL (see Figure 7-1) on FORTRAN unit IFTPRT. The 
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SAMEUIST/CONTPL/ 

IFTABL * 96 FQPTRAN UNIT NUMBER FOR DISPLAY TABLES 

IFTUDE a 23 FORTRAN UNIT NUMBER FOR THE GRAPHICS DEVICE 

IFTPRT = 6 FORTRAN UNIT NUMBER FOR THE PRINTER 

IROART =0 EQ !• READ IN MASTER NUMBERS 

NE 1. DO NOT READ IN MASTER NUMBERS 

IRDXST =0 EQ !• READ IN CISPLAY STATUS FLAGS 

NE 1. DO NOT READ IN DISPLAY STATUS FLAGS 

IDIREC =0 EQ 0. DO NOT READ IN DIRECTORY ARRAYS 

NE 0. READ IN DIRECTORY ARRAYS 

IRDTPO =0 EQ 0. DO NOT READ NAMEL IST/TPDSET/ 

NE 0. READ NAMELIST/TPDSET/ 

-0 NUMSUa * 3 number of major SUBSYSTEMS 

I 

o 

-i NUMCNC * 0 NUMBER CF CONVERSION AND CORRECTION ROUTINES 

NUMSCA = 0 NUMBER OF- SPECIAL CAPABILITY ROUTINES 

KOFFEE =0 EQ 0. TERMINATE AFTER FINE ATTITUDE 

EQ I. RECYCLE TO TELEMETRY PROCESSING 

GSATIO = ♦♦♦♦♦♦♦♦ SATELLITE IDENTIFICATION 

CRUNIO = «*♦♦*♦♦* RUN IDENTIFICATION 

I * 

IHMON » 9 MONTH OF REFERENCE DATE 

IROAY = I CAY OF REFERENCE DATE 

IRYEAR =< I9S7 - YEAR CF REFERENCE DATE 


Figure 7-1. NAMELIST CONTRL 
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only other printout from GESS consists of printed copies of character displays; 
see Section 7. 14 for examples of these displays. 

In addition, the GESS Executive will generate a core dump on the data set 
SYSPRINT if the MSAD/ OABIAS System abends in a graphics mode and the 
GESS Executive is able to intercept the abend. A DUMMY DD card can be 
provided for this data set to minimize the time required to recover from an 
abend in the graphic mode. 

7.12.2 Printed Output From OADRIV 

The only printed output from OADRIV consists of a listing of the parameters in 
NAMELIST MAIN on FORTRAN unit 6 (see Figure 7-2). In a nongraphic mode, 
this printout is generated following the reading of each new set of NAMELIST 
parameters for NAMELIST MAIN; in a graphic mode this printout is generated 
following each display at which the operator could have changed any parameters 
in NAMELIST MAIN. 

7.12.3 Printed Output From the ODAP Subsystem 

Printed output from ODAP includes all printout described in Reference 20. 

7. 12. 4 Printed Output From the OASYS Subsystem 

Printed output from OASYS includes all printout described in Reference 8, with 
the exception of the printout from GCONES. 

7. 12. 5 Printed Output From the OABIAS Subsystem 

The printed output of OABIAS results consists of tables and plots of the state 
vector, other quantities relevant to the recursive estimator, and the statistics. 
The printout is controlled by the NAME LIST parameters, ILEVEL, IPRINT, 
IDUMPL, ERUNIT, INTOUT, and IPLOT(12). The printout generated by 
OABIAS in the sample problem (Section 7, 3. 8) is shown in Figures 7-3 through 
7-18. 
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Figure 7-2. NAMELIST MAIN 
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7. 12. 5. 1 Input to OABIAS Table 

The input to OABIAS is printed on FORTRAN unit number IPRINT. At the 
beginning of each pass through OABIAS, all necessary parameters in BIASNL 
are printed. This printout can be eliminated by setting ILEVEL to zero. 

Figure 7-3 shows the input to OABIAS table generated by the sample problem. 

7. 12. 5.2 OABIAS Attitude Summary — Heading Definition 

The definition of the column headings of subsequent tables is printed on 
FORTRAN unit number IPRINT when ILEVEL is greater than 0 . Figure 7-4 
shows an example of this table. 

7.12.5.3 Attitude and State Vector Table 

The attitude and state vector table includes the attitude and biases estimated 

at each frame and the final estimate of the uncertainties. This table is printed 

on FORTRAN unit number IPRINT when ILEVEL is greater than or equal to 

2. The first and last pages of this table generated in the sample problem are 

shown in Figure 7-5. Note that the right ascension and declination of the spin 

axis are printed rather than the state vector components, S and S„ . 

X u 

7.12.5.4 Uncertainty Table 

The uncertainty table includes the square root of each diagonal element of flie 
covariance matrix, P (uncertainty), estimated at each frame. This table is 
printed on FORTRAN unit number IPRINT + 4 when ILEVEL is greater than 
or equal to 4. Figure 7-6 shows the uncertainty table generated in the sample 
problem. Note that the uncertainty in the right ascension and declination are 
printed rather than those in the state vector components, and S^. 

7. 12.5.5 Partial Derivative Table 

The partial derivative table includes tiie partial derivatives, G , of tlie observa- 
tion with respect to each element of the state vector evaluated at each frame. 
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Figure 7-3. Input to OABIA.S Table (1 of 3) 
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Figure 7-3. Input to OABIAS Table (2 of 3) 
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PRELIMINARY DRAFT 


Figure 7-3. Input to OABIAS Table (3 of 3) 
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Figure 7-4. OABIAS Attitude Summary — Heading Definition 



PRELIMINARY DRAFT 


i 


©H: 


2 


a it 

o a 

ti. H • 

<J u 


»• tfl 

D t- < 

o < a. 

Im 

> lA 

a. • 
< o 

X X. 

X < X 
3 O 
u> oj a 
3 ^ 
!A 3 01 


0 — 
0 — 


X X 
— s 
*- £ 


9000000000000000000000000900090009000009CIOOOOOOOOOOO 

9999900999999099909999999999099999999909099999999990 

J 





o 



N 


N 



■4 

^ 9 

9 

V 

9 

a 

a 

H 

o 

a 

a 

a 




y 

y 

y 




9 


9 

9 

9 

• 90 

> 

9 

9 

9 

9 

y 

a 

a 

1 

A 

IX 

3 




o o 

9 

9 


o o 

O 

o 

O 



9 

O 

9 

9 

O 

9 

o 

9 

9 















9 

9 O 

9 

9 

O 

9 

9 

o 

9 

9 

O 

o 

+5 

o 

o 

o 

9 

O 

o 

o 

O 

9 

9 

o 

o 

o 

O 

9 

9 

9 

O 

9 

O 

9 

9 

9 

9 

o 

9 

9 

O 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 9 

O 

O 

9 

O 

9 

o 

O 

9 

O 

o 

© 

O 


o 

o 

0 

1 

O 

0 

1 

0 

1 

0 

1 

o 

{ 

O 

9 

1 

9 

1 

0 

1 

9 

1 

O 

1 

9 

1 

9 

1 

9 

1 

9 

1 

9 

1 

9 

9 

1 

0 9 

1 t 

9 

t 

0 

1 

0 

1 

9 

1 

9 

1 

9 

0 

1 

9 

1 

9 

1 

O 

9 

1 

O 

1 

9 

1 

9 

1 

0 

1 

6 o 
1 1 

9 

t 

O 

1 

0 

1 

9 

1 

O 

t 

9 

1 

9 

1 

9 

1 

o 

} 

9 

i 

© 

s 

Q 

1 



o 

♦ 

9 


<t 

« 

* 

m 


n 


O 


•» 

O 

m 

9 

K 

♦ 

« 

a 

<r 

y 

y 

•n 

a 

0 

a 


9 



♦ 

n 

a 

a 

A 

A 9 

9 

9 


•> 

9 

# 

a 

a 

9 

9 





o 

o 

o 

o 






9 

N 

n 

n 

<y 

y 

cy 

ty 



O 



♦ 

a 


♦ 

♦ 



a 


n 

f> 

y 

y 





y 

y 

y 

y 

y 

y 

y 

y 

a 

n 



o 

o 

o 

o 

o 

o 

o 

o 

9 

O 

9 

9 

9 

9 

9 

9 

o 

9 

9 

9 

o 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

a 

9 9 

9 

o 

9 

9 

9 O 

9 

9 

9 

o 

9 

9 

o 

9 

9 

9 

© 

O 

'U 

a i 

o 

o 

o 

9 

o 

o 

o 

o 

O 

9 

9 

O 

O 

9 

9 

9 

o 

O 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

O 

9 

9 

9 

9 

9 

9 

9 

o 

O 

9 9 

O 

O 

O 

o 

9 

9 

9 

9 

9 

9 

9 

9 




CO 

9 

9 

9 


9 

9 

9 

9 

9 

9 


n 

*n 

n 

m 

<y 

cy 

<y 

y 

y 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

pw 

f- 




f. a 


a 

a 


a 

pik 

fk 

pw 

a 



Pk 


<0 


<0 

<0 


c 

9 

<0 

•0 

O 

0 

O 

4 

S 


Pk» 




fy 


h> 

r* 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

■a 

a 

a A 


a 

a 


3 

a 

a 

a 

a 

Q 

a 

a 


<o 


■Q 

9 


o 

0 


o 


<e 


O 

0 

9 

•0 

o 

•0 

•O 

a 

a 

a 

o 

a 

■0 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

9 

a a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 




9 

9 

o 



o 

V 

'if 

■o 

'<1 

o 

9 

O 

•tf 

o 

<u 

•c 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a 

a a 

a 

a 

V 

a 

a 

a 

a 

V 

a 

V 

a 

a 

tt 1 


•• 



•• 



■* 

■* 


*• 





** 


■* 

•* 






"* 

*“ 

* 

“* 

** 


■" 



*• 


•• 


* 

• •• 

«• 

■“ 



“ 



•• 



"* 

•" 


o 

o 


N 

9 

in 

n 

M 


9 


o 

0 

<0 

<y 

M 


•> 

<y 


9 

9 

9 

9 

9 

9 

a 

y 

y 



a 

9 9 

9 

y 


9 

o a 

a 

a 

9 

y 

N 

y 

N 


9 

_ 

y 

y 


o 

o o 

O 

O 





o 

o 

.•y 

•y 

N 

y 

N 




o 

9 



4t 

♦ 

A 

m 

m 

UP 

<n 

<r 

♦ 

n 

*r 

4*» 

n 

fM 

y 

y y 

y 

y 

y 

a 

fl 


fl 

a 

y 





a 

o 

o 

9 

o 

o 

o 

9 

9 

9 

9 

o 

9 

9 

Vf 

9 


9 

9 

9 

9 

9 

V 

9 

O 


V 

w 

u 

9 

9 

9 

u 

o 

9 

9 

O 

9 

9 9 

9 

o 

o 

u 

O 

u 

o 

9 

9 

9 

o 

o 

0 

a 

0 

1 

CJ 

0 

1 

0 

1 

o 

♦ 

0 

1 

0 

1 

0 

1 

9 

1 

O 

1 

O 

9 

I 

O 

t 

9 

1 

9 

1 

O 

1 

9 

1 

9 

I 

0 

1 

9 

t 

O 

1 

9 

1 

9 

1 

9 

1 

o 

o 

9 

1 

0 

1 

o 

e 

I 

9 

1 

O 

» 

0 

1 

o 

9 

1 

c 

1 

0 

1 

9 

1 

9 O 
1 1 

9 

1 

9 

1 

0 

1 

0 

1 

0 

1 

o 

t 

9 

1 

9 

1 

O 

1 

O 

f 

0 

1 

0 

1 




9 

9 

■0 ^ 

N 


9 

<y 

-1 

9 

•n 

tn 

O 

9 

9 

K 

a 

<y 


a 

'0 

9 

9 

9 

a 

0 

a 

A 



n 

9 

'll 


a 

a 

(0 O 


a 

a 

9 


a 

A 



9 

K 

fk 




o 



o 

o 

O 

O 

o 

9 




n 

r» 

<y 

ty 



9 



y 

y 


n 

n 



f| 


y 

n 


y 



*• y 


y 


a 

a 

y 

y 







o 

o 

o 

O 

o o 

o 

o 

O 

9 

O 

9 

o 

o 

9 

o 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

o 

9 

9 

o 

9 

O 

9 

9 

9 

9 

9 

o 

9 O 

9 

o 

o 

9 

o 

o 

9 

9 

o 

9 

o 

o 

u; 

3 

o 

o 

0 

1 

0 

1 

0 

1 

0 

1 

0 

1 

9 

0 

1 

0 

1 

9 9 
1 

9 

9 

9 

9 

9 

9 

9 

9 

9 

9 

O 

9 

O 

9 

9 

9 

9 

o 

9 

9 

9 

O 

9 

9 

9 

9 

o O 

a 

O 

o 

9 

9 

9 

O 

O 

o 

9 

9 

9 




o 



< 

# 

!•! 




•J» 

o 

O 

<y 

•y 

a 

a 

O 

a 

<Ji 

yy 

a 


y 

»n 

y 

■o 

<S 

a 

a 

a 


« 

♦ 

9 

y 

o 

o a 

■c 

a 

o 

9 

0* 

y 

y 

9 

A 

9 


y 




o 



nl 


py 

N 

M 


m 

n 

•n 

•n 


<y 

y 

o 

o 

O 



a 

■n 

U) 

a 

•A 

A 

a 

u> 

A 

•A 

a 

O 

a 


«' 


9 

a 

a 


•o 

a 


a 

A 

P« 

© 

•35 


o 

o 

o 

o 

o 

o 

o 

ij 

O 

o 

9 

9 

o 

9 

9 

9 

o 

o 

9 

o 

9 

9 

9 

9 

9 

9 

o 

9 

9 

Q 

a 

9 

O 

9 

9 

9 

o 

9 

9 9 

a 

o 

o 

u 

9 

9 

9 

o 

O 

o 

9 

O 

U 1 

o 

o 

o 

a 

o 

o 

o 

o 

o 

o 

O 

9 

o 

9 

9 

O 

o 

9 

9 

9 

9 

1 

9 

O 

9 

o 

9 

o 

9 

9 

O 

o 

9 

o 

9 

O 

O 

9 

O 

9 9 

o 

9 

o 

9 

9 

9 

9 

9 

9 

9 

O 

© 

















9 

9 

M 

<n 



_ 


9 

1. 

9 

o 



9 

y 

y 

n 

y 

a 

*<• 

a 

PW 

IX Ml 




9 

9 

•n 

A 

A 


y 

y 

y 

















9 

o 

9 

o 

9 

O 



N 

y 

y 

in 

n 

n 

M 

y 

y 

9 

o 


p« 

o 

9 


o 

9 


O 

9 





y 

y 

y 


o 

o 

o 

o 

o 

o 

9 

o 

9 

o 

9 

C 

c 

9 

o 

O 

o 

9 

9 

9 

c 

9 

O 

O 

o 

o 

a 

9 

o 

9 

9 

o 

o 

o 

9 

O 

o 

O 

9 O 

o 

O 


9 

9 



O 


o 








« 










• 










• 



• 


• 

t 




• 




• * 














o 

o 

o 

o 

o 

9 

O 

o 

O 

o 

9 

9 

9 

o 

9 

o 

9 

9 

9 

O 

1 

9 

1 

0 

1 

0 

1 

o 

o 

9 

9 

o 

9 

O 

O 

9 

9 

9 

9 

1 

9 

1 

9 

1 

0 

1 

9 9 
1 1 

0 

1 

9 

I 

9 

1 

9 

1 

9 

1 

9 

1 

9 

1 

0 

1 

O 

i 

0 

1 

0 

1 

o 




9 

9 





9 




■0 

(D 

N. 

K 

A 

tf> 


a 

O 

a 

9 


N> 

0 

y 




9 

9 

y 

a 



•n 

a 

If 

a 

y 


y 




a 

9 







o 



O 

o 

o 

O 

o 

e 




<y 

«y 

y 

fy 



O 



y 

y 

y 


n 

n 

f» 

y 

n 

y 

y 


y 





M 


y 


y 








o 

o 

o 

o 

o 

O 

o 

o 

O 

o 

9 

9 

O 

O 

o 

9 

9 

9 

9 

9 

9 

o 

9 

o 

o 

o 

9 

9 

9 

9 

9 

9 

o 

o 

O 

o 

9 

o 

9 9 

o 

9 

o 

9 

O 

o 

9 

9 

O 

© 

9 

o 

uj 

3 

o 

o 

o 

o 

o 

o 

o 

o 

f 

9 

1 

9 

9 

9 

1 

0 

1 

9 

1 

0 

1 

0 

1 

9 

t 

9 

1 

9 

9 

1 

9 

1 

9 

1 

O 

1 

0 

1 

9 

1 

9 

1 

9 

1 

0 

1 

O 

1 

9 

9 

1 

9 

1 

9 

1 

9 9 
1 1 

9 

!• 

0 

1 

o 

» 

9 9 
1 1 

0 

1 

9 

1 

9 

1 

9 

1 

9 

1 

9 

1 

0 

1 

9 

1 

O 

H 

0 

1 

O 

J 

© 

t 



n 

fS 





> 

9 


9 

♦ 


9 

•A 

♦ 

O 

a 

a 


K 

K 

P* 




9 



•a 

f» 

>A a 

y 

a 

a 

9 

pn 

a a 

4T 

a 



A 



9 

9 

fk» 

9 

9 


<« 

■* 

O 




M 

M 



CM 


a 

a 

*0 


9 

> 

cy 

N 

y 

a 

a 

y 

y 

y 


a 

9 

A 


a 

a 



9 

N 

N 



o 




A 


9 






o 

o 

o 




9 

IW 




<0 


n 

o 

c 

“• 

o 

* 

— 

o 


9 

a 

a 

a 

a 

a 

a 

>» 

<r 

« 

a 

a 

<« 

A 

A 

<9 

A a 

9 

9 

9 

a 


« 

9 

9 

9 

n 

a 

a 

El 

N 

«M 

cy 

A 


■M 






-• 

m 


9 

9 

o 

o 

9 

O 

9 

'I* 

> 

y 

y 

y 

y 

y 

y 

y 

y 

y 

y 

y 

y 

y 

y 

y 

y y 

y 


y 



y 

y 

y 


y 

y 


1 

lA 

>/) 

lA 


A 


'0 

<0 

rf» 



<n 



lO 

A 

a 

a 

a 

a 


4T 

n 

fl 

A 

■A 

tf) 

«) 

A 

lA 

9 

in 

(A 

•A 

A 

A 

lO 

A A 

A 

A 

A 

A' 

'a 

A 

A 

A 

A 

A 

A 

A 

1 

o 

O 


•0 


0 

9 


> 



o 

<> 

n 

9 

9 

« 

<f 

a 

m 

m 

0 

>n 

a 

_ 

a 

K 

Q 

-T 

♦ 

■n 

♦ 

fi 

a 

K 

a 

a 

_ 



O 



9 

a 


y 

A 

9 

> 

a 




*• 






-• 


N 

«• 





> 


o 

9 

y 



o 

9 

a 




C 

o 

a 

0 

y 

y 

a 

h- 

fl* 

s a 

a 

fl 

a 


9 

y 

a 

A 

9 

f- 


Ik. 

■ 

9 

y 

n 

nj 


»n 

(M 

fy 

n 

•n 

« 

<0 

-0 

<0 

O 

■0 

a 

0 

a 

a 

9 

■“ 

-* 

- 

— 

9 

9 

> 

O' 

o 

9 

> 

A 

O' 

3* 

a 



p* fk 


X 

h* 

9 

« 

a 

9 

9 

9 

Ok 

0* 

9 

aJ 

•> 

9 

T» 


> 


> 

9 


9 

9 



3k 

J» 

> 

> 

9 


o 

O 

9 

9 

9 

9 

9 

9 

> 

> 

o 

e 

> 

■> 

9 


9 

> 

9 

a* » 


9 9 

9 


> 




s-> 


9 

a 

1 

1 

1 

7 

1 

1 

1 

i 

1 

1 

1 

»- 

1 

1 

1 

1 

1 

i 

\ 

» 

y 

1 

y 

1 

y 

1 

y 

1 

y 

t 

y 

1 

1 

• 

1 

y 

y 

1 

♦ 

1 

1 

1 

1 

i 

1 

1 i 

1 

1 


1 

1 

1 

7 

1 


1 

1 



tfl 

lit 


■# 

PV 

n 

O 

9 

9 


V 


o 

♦ 


a 

•s 

ry 

r» 


K 

_ 

<y 

y 

9 

a 

9 


u* 

n 

9 


n 

•A 

9 

a 




p>» 

a 




r* 


9 







♦ 

9 

O 

N> 


IN 

-• 

(\J 

n 

a 

9 

9 

O 

9 

u 

<n 



■3* 

a 

9 

TJ 

M 


o 



f» 

w 

9 

9 

o» 


■J* 

pk. 

A 

A 

•-?> ♦ 

<r 


a 

9 


a 



■0 

Ok 


9 



-• 

9 

9 

9 

> 

9 

9 

> 

9 

9 

o 

9 

9 

O 

c 

9 

9 

o 

a 


•• 



•« 

«• 

9 

O 

o 

c 

o 

U» 

o 

9 

O 

9 

9 

O 

O 9 

o 

o 

o 

9 

9 

o 

O 

9 

o 

o 

o 

© 





















































LJ 



O 

O 

9 

c 

O 

O 

C* 

O 

O 









































3 


in 

in 

n 

•n 


(0 


n 

A 


K) 

wi 

n 

.1 

•ji 

A 

a 

a 

m 

a 

a 

yy 

a 

A 

141 

41 

(A 


A 

A 

A 

4» 

.1 

41 

U) 

A 

•41 

A A 


A 


A 

A 

A 


14) 

A 


A 

A 




•> 


•A 


-» 


9 

-> 


■n 


m 

•n 

<n 

•J 


•n 

«n 


)•> 

m 

n 

fl 

fi 

n 

n 

fl 

fl 

f| 

n 

-I 

n 

n 

n 

n 

m 

a a 

a 

a 

a 

a 

a 

a 

a 

f) 

n 

a 

a 

a 


o 

5 


n*n<n'nr»n 4»>4>a>'uv«i4>nirvi 


9 -f < tt ♦ 

\\NSNX'VVNX 

♦ ♦ ♦ ♦ ■ 


,NXV^.XNS^SN 


S N X N "V N 

n n n •* — 

m oj N 


xxxnxxssnnxnsx’vsxvnxnnnxxnsxnxnnxnsv 

0«(M.NM<MNOI'4n4N^(VMMNN0tN<V(VNNMAin|('4r«N(^(g9N«<|O«<^M 




oS o 

3f - 


7-75 

PRELIMINARY DRAFT 


Figure 7-5. Attitude and State Vector Table From OABIAS (1 of 2) 
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Figure 7-5. Attitude and State Vector Table From OABIAS (2 of 2) 
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This table is printed on FORTRAN unit number (IPRINT +2) when ILEVEL is 
greater than or equal to 6. Figure 7-7 shows the partial derivative table 
generated in the sample problem. 

7. 12.5.6 Recursive Estimator Gain Table 

The recursive estimator gain table includes the elements of the gain vector, K , 
for each element of the state vector estimated at each frame. This table is 
printed on FORTRAN unit number (IPRINT +1) when ILEVEL is greater than 
or equal to 8. Figure 7-8 shows the recursive estimator gain table generated 
in the sample problem. 

7.12.5.7 Attitude, Observable, and Vector Table 

The attitude, observable, and vector table consists of the attitude, the residual, 
the inverse of the weight, the observation, the computed observation, and the 
three components of the Sun and position vector at each frame. This table is 
printed on FOR'^RAN unit number (IPRINT +3) when ILEVEL is greater than 
or equal to 10, Figure 7-9 shows the attitude, observation, and vector table 
generated in the sample problem. 

7. 12. 5. 8 Residual Table 

« 

The residual table includes the mean residual and the standard deviation of the 
residuals computed at each frame. This table is printed on FORTRAN unit 
number (IPRINT +5) when ILEVEL is greater than or equal to 12. The first 
and last page of the residual table generated in the sample problem are 
shown in Figure 7-10. 

7.12.5.9 Error Messages 

The error messages are printed on FORTRAN unit number ERUNIT when 
IDUMPL is greater than or equal to 4. Some of the messages generated in the 
sample problem are shown in Figure 7-11. See Section 7.16.4 for a detailed 
description of the error messages. 
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Figure 7-8. llecursive Estimator Gain Table From OABIAS 
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Figure 7-9. Attitudes Observables and Vector Table From OABIAS 
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Figure 7-10. Residual Table From OABIAS (1 of 2) 
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Figure 7-11. Error Messages From OABIAS 
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7. 12. 5. 10 Error Summary Count Table 

The error summary count table contains the number of errors that occurred in 
each subroutine during the pass. This table is printed on FORTRAN unit 
number ERUNIT. If IDUMPL is set to zero, only this table is printed and the 
error messages shown in Figure 7-11 are not printed. Figure 7-12 shows the 
error summary count table generated in the sample problem. See Sec- 
tion 7. 16.4 for a description of the error messages corresponding to the loca- 
tions in the error summary count table. 

7.12.5.11 State Vector Plot 

Each element of the state vector is plotted on FORTRAN unit number IPRINT 
if IPLOT for the corresponding element is 1 or 3 . As examples, Figures 7-13 
and 7-14 show the plots of the right ascension and declination of the spin axis 
generated in the sample problem. 

7. 12. 5. 12 Uncertainty Plot 

The square root of each diagonal element of the covariance matrix (uncertainty) 
is plotted on FORTRAN unit number IPRINT if IPLOT for the corresponding 
element is 2 or 3 . As examples, the uncertainty plots for the right ascen- 

4 

Sion and declination of the spin axis generated in the sample problem are 
shown in Figures 7-15 and 7-16, respectively. 

7. 12. 5. 13 Residual Plot 

The residual for each model is plotted on FORTRAN unit number IPRINT if 
RESMODfor the corresponding model is .TRUE. . 

As examples, tlie residual plots for the Sun angle model and single horizon di- 
hedral angle model are shown in Figures 7-17 and 7-18, respectively. The 
mean residual and standard deviation at the end of the block are also shown in 
these plots. 
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Figure 7-12. Error Summary Count Table From OABIAS 
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Figure 7-13. OABIAS Filter Response in Right Ascension of Spin Axis 
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Figure 7-14. OABIAS Filter Response in Declination of Spin Axis 
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Figure 7-15. Estimated Uncertainty in Right Ascension of Spin Axis 
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Figure 7-16. Estimated Uncertainty in Declination of Spin Axis 
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Figure 7-18. Residuals of Single Horizon Dihedral Angle Model 
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7. 13 CALCOMP PLOT TAPE 

The GESS Executive can be used to generate a CalComp plot tape of the GESS 
plot displays. If CalComp plots are desired, a DD card must be included for 
the tape. See Section 7. 18 for a description of the DD card. 

Following termination of the job, the tape must be taken to the CalComp plotter 
in the basement of Building 1 (GSFC). The user should specify plain white pa- 
per, 12 inches wide, for the plots. The number of CalComp files on the tape 
should be specified as 999, CalComp plotting time will be approximately 5 min- 
utes per plot. A maximum of about 200 plots will fit on 1 reel of tape. 

7. 14 GESS DISPLAYS 

This subsection describes each GESS display and provides a hard-copy example. 
For each control parameter display (e.g. , NAMELIST display), the hard-copy 
example is annotated with the variable name associated with each parameter, 
unless the descriptor line on the display already includes the variable name. 
Note that certain integer control flags are displayed as alphanumeric fields 
rather than integers (e, g. , the parameter lOBLAT is displayed as YES or NO 
rather than 1 and 0). In each case where this occurs, the descriptor line 

I 

on the display lists all possible alphanumeric keys which are valid. The hard- 
copy is amiotated with the corresponding parameter values in the same order. 

The displays are described in the order in which they normally appear. This 
order should be noted so that the operator can move efficiently forward or 
backward to the desired display in a graphic mode. 

« 

7. 14. 1 Display Status Flags and Key Assignments 

Figure 7-19 shows the display status flags and key assignments display gen- 
erated from the GESS Executive. The user may alter the status of any control 
point from this display, as described in Reference 13. The key assignments. 
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Figure 7-19. Display Status Flags and Key Assignments (1 of 3) 
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shown on tlie second page of Uiis display, are for information only. The user 
may move asynchronous Ij’- to this display at any time by depressing key 30 . 

7.14.2 Array Allocation Sizes 

Figure 7-20 shows the array allocation sizes display. The user may change 
array allocation sizes from this display. The values of 200 for OABIAS "Data 
Records" may be decreased, if desired, to save core but should not be in- 
creased to more than 200. The value of 300 for "Input Data Records" may be 
increased to allow additional input data records for data selection purposes. 
Approximately 1300 input data records is the maximum that can be displayed 
by the 2250 display device. The value of 20 for RAE-B "Header Records" may 
be increased, if desired, to allow saving more than 20 header records. The 
user may return to this display at any time by depressing key 0. The other 
array sizes should not be changed by the operator. 

7. 14. 3 Main Control Display 

Figure 7-21 shows the main control display generated from OADRIV, This 
display contains the parameter ISIM from NAMELIST MAIN described in 
Section 7.3.4. 

4 

The user may return to this display at any time by depressing key 1 . 

7.14.4 Options for Copying AE Data Set 

Figure 7-22 shows the display of options for copyii^ the AE data set to the 
OABIAS data set format. This display includes parameters from NAMELIST 
MAIN, described in Section 7.3.4. 

7. 14. 5 Simulator Options Display 

Figure 7-23 shows the simulator options display. This display contains simula- 
tion parameters in NAMELIST MAIN, described in Section 7.3.4. 
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Figure 7-20. Array Allocation Sizes 
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Figure 7-21. Main Control Display 
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Figure 7-22. Options for Copying AE Data Set 
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Figure 7-23. Simulator Options 
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7.14.6 NAMELIST LIST Display 

Figure 7-24 shows the display of parameters ia NAMELIST LIST. See Sec- 
tion 7.3.5 for a description of these parameters. 

7. 14. 7 Options for Reading Data 

Figure 7-25 shows the display of options for reading the OABIAS data set. 
This display includes parameters from NAMELIST MAIN, described in Sec- 
tion 7.3.4. 


7. 14. 8 Header Record Displays 

Figure 7-26 shows the three displays which contain data from the OABIAS 
header records. For a description of the data on the header records, see 
Section 7.4. 


7.14.9 Data Record Display 


Figure 7-27(1) shows the OABIAS data records display. The data on this dis 
play are described in Section 7.4, with the following exceptions: 


Heading 


Description 


REC. NO. 


Record number on data set 


FRM NO. 

FLAG 

First character 


Second character 


Third character 
Fourth character 


Frame number, the index in the internal 
arrays (1-200) 

A four-character flag; 

= I, inversion case (out-triggering 
occurred before in- triggering) 

= Blank, normal case 
= X, inversion flag invalid (not 0 or 1) 

= 1-8, configuration flag 
= X, configuration flag invalid (not 1-8) 

= E, M, S, U, central body flag 

= H, T, U, terminator/horizon flag 
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Heading Description 

R Data rejection flag (set automatically by 

OADRIV, and may be changed by operator): 

= valid frame 

^ flagged frame 

The following character is used by OADRIV 
for flagging data: 

= X, configuration flag invalid (not 1-8 or 
central body flag invalid (not E, S, 
or U) 

Following the data records display are six sets of plots of the data. Each set 
consists of a plot of data versus time followed by a plot of data versus frame 
number. The data plotted are the same as those tabulated in the data records 
display. Points may be flagged from the plots by use of the GESS FLAGPT 
options. 

The light pen is used to touch the FLAGPT option, each point to be flagged, 
and the END option. The six plot displays are as follows: 

1. Sun angle (degrees) 

2. Spin rate (degrees per second) 

3. Rotation angle. Earth-in (degrees) and Earth-out (degrees) 

4. Earth-in (degrees) 

5. Earth-out (degrees) 

6. Earth width (degrees) = Earth-out minus Earth -in 

Figure 7-27(2) shows the plot of rotation angles versus time. The other plots 
are similar. 

7. 14. 10 Data Selection Options for OASYS 

Figure 7-28 shows the display of data selection options for OASYS. These 
parameters are in NAMELIST MAIN described in Section 7.3.4. 
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Figure 7-25. Options for Reading Data 


PRELIMINARY DRAFT 



PRELliVliNARY DRAFT 


»»»»»» -n ^ A O * * 

0 1 S P L AY *444444*4*4 73 • 29 1 • tO • * t • 52 4*** 


44 FAEI-DI <4EA0ER RECORDS* DISPLAY t OF 3 

44 

44 FIRST LASf START TIME JOB TIME 

- 4* REC ' RS: YYPMOO*MHMMSS YYMHDD.MMMMSS 

44 I 740104*230400 730102*030405 

44 


-3 

I 


O 

C5 


- 4* 
44 

44 

-44— 

44 

44 

44 

-44— 

44 

44 

44 

44 

44 

4* 

44 

44 

44 

4 4 

4* 

44 

44 

4* 

44 

-44 - 
44 
44 
44 

4* 

4* ' 

44 

44 


44 

44 


U9CI ENC CF FILE ON O^ilAS DATA SET 
CFCINI^FAOER WHAT NJ« 


CALL OISPLAT 


OISP I OF 


4444 
44 44 


444 44-444444444*4444*44444 4 44444 
444*4#4444444444*44*444*4444444 


M 3 A O 
I S P L A Y 


44 

4 4 

4 4 44444444444444444444444*44444444 
****•44 44 444* 44 444 4444*44444444444 
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Figure 7-26. Header Record Displays (2 of 3) 
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Figure 7-26. Header Record Displays (3 of 3) 
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Figure 7-28. Data Selection Options for OASYS 
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7.14.11 NAMELIST QPMANl Display 

Figure 7-29 (1 and 2) shows the NAME LIST OPMANl display. These param- 
eters are described in Section 7.3.6. Figure 7-29(3) shows the source of the 
orbit to be used for processing (ephemeris tape, GTDS file, or orbit generator) 
and the orbital elements from either the header record of ephemeris tapes or 
GTDS files or the input parameters to the orbit generator. 

7.14.12 OASYS Block Average Display 

Figure 7-30 shows the display of block average results from OASYS. These 
values are stored in COMMON AVRAGE. 

7. 14. 13 OASYS Plots of Single Frame Results 

There are seven plot displays of the single frame results from OASYS. Each 
plot includes up to four curves, corresponding to the four different attitude 
computation methods in OASYS. If the GESS IDENTIFY option is used (see 
Reference 13), a number will appear on each curve to identify it as follows; 

= 1, single horizon method, in 
= 2, single horizon method, out 
= 3, double horizon width method 
= 4, double horizon dihedral method 

The seven plots are as follows: 

1. Alpha, including rejected points — This plot shows right ascension 
alpha, versus frame number, plotting the selected attitude from 
each ambiguous pair and including attitudes rejected in block aver- 
aging. 

2. Alpha, omitting rejected points — This plot is the same as plot 1, 
but attitudes rejected in the block average are omitted. 
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Figure 7-29. NAMELIST OPMANl (2 of 3) 
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Figure 7-30. OASYS Block Average Results (1 of 2) 
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Figure 7-30. OASYS Block Average Results (2 of 2) 
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3. Delta, including rejected points — This plot is the same as plot 1 
but is for declination. 

4. Delta, omitting rejected points — This plot is the same as plot 2 
but is for declination. 

5. Arc length uncertainty — This plot shows arc length uncertainty in 
attitude versus frame number. Points for which the uncertainty is 
undefined are omitted, but points rejected on the block average are 
included. 

6. Nadir angle plot — This plot shows nadir angle versus frame num- 
ber, including points rejected in the block average. 

7. Dihedral angle plot — This plot is the same as plot 6 but is for di- 
hedral angles. 

Figure 7-31 shows an example of plot 1. The other plots are similar. 

7.14.14 NAMELIST BIASNL Display 

Figure 7-32 shows the NAME LIST BIASNL display. These parameters are 

described in Section 7. 3. 7. 

7.14.15 Final Results From OABIAS 

Figure 7-33 shows the display of final results from OABIAS. These param- 
eters are stored in COMMON TABCOM. 

7. 14. 16 Table of Errors From OABIAS 

Figure 7-34 shows the display of errors from the OABIAS Subsystem, The 

errors are described in Table 7-1, Section 7.16.4. 
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Figure 7-31. Alphas Versus Frame Number, Including Rejected Points 
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Figure 7-32. NAMELIST BIASNL (1 of 3) 
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Figure 7-32 
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Figure 7-32. NAMELIST BIASNL {3 of 3) 
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Figure 7-33. Final Results From OABIAS (1 of 2) 
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Figure 7-33. Final Results From OABIAS (2 of 2) 
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7. 14. 17 Table of Correlation Coefficients from OABIAS 

Figure 7-35 shows the display of correlation coefficients from the OABIAS 
Subsystem. Each value in the table is computed as follows; 


A.. 

13 



where P is the final covariance matrix. Since the matrix is symmetric , the 
lower diagonal portion is filled with zeroes on the display. 

7. 14. 18 State Component Plots From OABIAS 

Figure 7-36 shows one example of a state component plot. The state compo- 
nent right ascension is plotted versus frame number. ■ 

7.14.19 Plots of Uncertainties in State Com p jnents 

Figure 7-37 shows one example of a plot of the imcertainty on the state compo- 
nents. The imcertainty in right ascension is plotted versus frame number. 

7.14.20 Residual Plots From OABIAS 

Figure 7-38 shows one example of a residual plot. The residuals from the Sun 
angle model are plotted versus fx’ame number. 

7.14.21 Options for Data Prediction 

Figure 7-39 shows the display of options for data prediction. These param- 
eters appear in NAMELIST MAIN and are described in Section 7. 3.4. 

7. 14. 22 Plot of Predicted and Observed Rotation Angles 

Figure 7-40 is an example of the plot of predicted and observed rotation angles 
versus time. The X-axis is in units of minutes from the start time of the pre- 
diction. The plot shows dihedral angles from the Sun to Earth-in and the Sun 
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Figure 7-35. Correlation Coefficients 
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Figure 7-39. Options for Data Pre(3iction (1 of 3) 
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Figure 7-39. Options for Data Prediction (3 of 3) 
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to Earth-out. The observed rotation angles appear as points. The predicted 
rotation angles appear as lines. There are up to four segments for each pre- 
dicted plot; first horizon, first terminator, second terminator, and second 
horizon. There are up to three predicted plots, corresponding to the three sets 
of inputs in the plot option table shown in Section 7. 14.21. The three predicted 
plots are identified by the numbers 1, 2, and 3 when the GESS IDENTIFY is 
used to identify the curves. Numbers 4 and 5 in the identify option referred to 
observed Earth-in data and Earth-out data, respectively. The plot of predicted 
and observ'^ed rotation angles versus frame number is identical to the plot in 
Figure 7-40 except that the X-axis Is frame number. The frame number plot 
is only available if ITIME = 0 (see Section 7. 14. 21). 

7. 14. 23 Plot of Predicted and Observed Earth Widths 

Figure 7-41 is an example of the plot of predicted and observed earth widths 
versus frame number. The three predicted plots and the observed plot corre- 
spond to the same data as the rotation angle plots (see Section 7. 14. 22). The 
Earth width equals the Earth-out dihedral angle minus the Earth-in dihedral 
angle. The plot of predicted and observ^ed Earth width versus time is identical 
to the frame number plot, except that the X-axis is measured in minutes from 
the start time of the prediction. 

7. 14. 24 Core Storage and Time Remaining Display 

Figure 7-42 shows an example of the core storage and time remaining dis- 
play. This display is generated whenever the operator enters the command 
'CORTIM' into tlae WHAT NOW field of any display. The display indicates the 
amount of free core storage within the user's region, the extent of core frag- 
mentation, and the computer processing unit (CPU) and I/O time remaining in 
the job. 
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PTflNGF PREDICTED AND OBSERVED ROTATION ANGLES VS. FRAME NO. 
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FLAGPT 
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Mgure 7-40. Predicted and Observed Rotation Angles Versus Time 
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7. 15 CONTROL WITH INTERACTIVE GILVPIIICS 

This subsection describes the ways in which the operator can modify program 
flow in a graphic mode. The reader is assumed to have a knowledge of tlie gen- 
eral capabilities of GESS graphics systems (see Reference 13, pages 6-14 ff. ). 

7. 15. 1 Programmed Function Keys 

The MSAD/OABIAS System uses 20 programmed function keys, in addition to 
keys 0, 30, and 31 which are normally provided in any MSAD system. The 
functions of tliese keys are listed below; 

Key 

0 

1 

4 

5 

6 

7 

8 
9 

10 
11 
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Function 

Causes move to array allocation sizes display. 
Allows the operator to change array allocation 
sizes. Always a backward move to the beg innin g 
of the system 

Causes move to main control display in OADRIV 

Causes move to tabular display of input data 

Causes move to plot of Sun angle data versus time 

Causes move to plot of spin period data versus 
time 

Causes move to plot of rotation angle data versus 
time. (NOTE: If IPLOT = 1 in the data reading 
options displaj', then key 7 will cause a move to 
the Earth-in data versus time plot. See Sec- 
tion 7. 3.4. ) 

Causes move to plot of Earth width data versus 
time 

Causes move to display of options for selecting 
data for OASYS 

Causes move to NAMELIST OPMANl display in 
the OASYS Subsystem 

Causes move to character display of block aver- 
age results from OASYS 
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Key 

Function 

12 

Causes move to plot of alpha (including rejected 
points) from OASYS 

13 

Causes move to plot of delta (including rejected 
points) from OASYS 

14 

Causes move to plot of arc length uncertainties 
from OASYS 

15 

Causes move to plot of nadir angles from OASYS 

16 

Causes move to NAMELIST BIASNL display in the 
OABIAS Subsystem 

17 

Causes move to character display of final results 
from OABIAS 

18 

Causes move to plot display of first state compo- 
nent in OABIAS 

19 

Causes move to plot display of uncertainty in first 
state component in OABIAS 

20 

Causes move to plot display of residuals from 
first model in OABIAS 

22 

Causes move to display of options for data pre- 
diction in PLOTOC Subsystem (NAMELIST MAIN) 

23 

Causes move to plot of predicted and observed 
rotation angles versus time 

24 

Causes move to plot of predicted and observed 
Earth-widths versus time 

30 

Causes move to GESS XSTOPS display used to 
control graphics displays and check key assign- 
ments 

31 

Causes move to GESS PRINT; prints lineprinter 
plot of character displays and CalComp plot of 
plot displays 

NOTE; Once the OABIAS Subsystem has been entered, the operator can leave 
this subsystem only by making a baclovard move (usually either key 1 
or kej’- 10). In a graphic mode, the OABL\S Subsystem is in an infinite 
loop, always returning to the NAMELIST BIASNL display to accept 
new parameters for reprocessing. 
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The programmed function keys are assigned for convenient use with the 
2250 graphic display device. The four long rows on the key console corre- 
spond to the four major subsystems of OABIAS — OADRIV, OASYS, OABIAS, 
and PLOTOC. The keys are normally used as follows: if the operator desired 
to return to the main control display in OADRIV from any point in the system, 
(e.g. , to read new data or terminate), key 1 is depressed. Within OADRIV 
the first long row is used to move rapidly back and forth to examine data 
tables and displays, flag data, and select data for further processing. Within 
the first long row, key 4 shows the tabular data display, keys 5 through 8 show 
plots of the input data, and key 8 shows the options for adding noise or biases 
to the data and selecting data for further processing. The second long row 
controls OASYS processing and result displays. The first key in the row shows 
the OASYS processing parameters. The second key is depressed to obtain the 
tabular results of the OASYS processing. The remainder of the keys are used 
to move quickly back and forth among the plots of OASYS results. The third 
row is the OABIAS Subsystem and is similar in structure to the second row. 

The first two keys are processing param eters and final results table, respec- 
tively. The next three keys show plots of OABIAS results. The fourth long 
row controls the PLOTOC Subsystem. The first key of the row ‘shows the op- 
tions for data prediction. The next two keys show the predicted and observed 
rotation angle and Earth -width plots. Finally the last two keys on the console 
bring up the XSTOPS display and provide the PRINT option. 

See Section 2. 4 for a detailed example of controlling system flow in a graphic 
mode. 

7.15.2 Asynchronous Calls 

The following control point names may be used for asynchronous calls. For an 
asynchronous call, the control point name is entered into the WHAT NOW field 
of any display; the specified control point is then invoked. Following execution 
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of the actions specified at that control point, control returns to the display 
from which the asynchronous call was made. 


Control Point Name 

" Action 

HEADER 

Generate displays of OABIAS header 
records (three displays) 

DATDIS 

Generate tabular display of OABIAS 
telemetry data 

BETA 

Generate plot display of Sun angle versus 
time and versus frame number 

OMEGA 

Generate plot displa}’^ of spin rate versus 
time and versus frame number 

ROANGL 

Generate plot display of rotation angles 
versus time and versus frame number 

EWIDTH 

Generate plot displa}!' of Earth width 
versus time and versus frame number 

CORTIM 

Generate core storage and time remain- 
ing display 

ol point HEADER would be 

invoked asjmchronously only if the operator 


desired to look at the header record information. Control points DATDIS, 
BETA, OMEGA, ROANGL, and EWIDTH could be invoked asynchronously to 

I 

allow the operator to examine the data or to change the rejection flags. The 
rejection flags from the data display are re-examined each time processing 
is performed in OASYS or OABIAS, and the rejection status of the data is 
modified appropriately. 


Control point CORTIM can be invoked whenever the operator wishes to examine 
core fragmentation or the amount of CPU or I/O time remaining in the job. 
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7.16 ERROR MESSAGES 


7.16.1 Messages Displayed ia a Graphic Mode 

Two types of error messages appear in a graphic mode; 


appear on a single line at the bottom of the 
on a separate display. For messages not listed below. 


ADDMSG Messages are as follows: 

Message 

U700 REQUESTED IN- 
TERVAL COPIED. 

NOR]\IAL COMPLETION 
OF AECOPY 


U701 END OF FILE ON ^ 
AE DATA SET. NORMAL 
COMPLETION OF AECOPY 


U702 END OF FILE ON 
AE DATA SET. RE- 
QUESTED INTERVAL NOT 
FOUND 


U703 NO RECORDS ON 
AE DATA SET 


U801 END OF FILE OR 
I/O ERROR. NAMELIST 
/BIASNL/ 


ADDMSG messages 
display; MESAGE messages appear 
see Reference 13. 


Description 


The requested block(s) fi'om the 
AE data set have been copied to 
the OABIAS data set. At least one 
frame was copied 

User response: None 

An end-of-file was encountered on the 
AE data set before the end of the re- 
quested interval was encountered. At 
least one frame was copied, 

User response: None 

An end-of-file was encountered on 
the AE data set before the beginning 
of the specified interval was found 
No records were copied 

User response: Specify a valid interval, 
or an interval beginning with block 1 

An end-of-file was encountered 
while attempting to read the first 
block from the AE data set. The 
AE data set is empty. No records 
were copied 

User response: Obtain a valid AE 
data set 

An end-of-file or I/O error was 
encoimtered while reading 
NAMELIST BIASNL in INITIL. 

See description of U920, below 
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Message 


Description 


An end-of-file was encountered on 
the OABIAS data set before reading 
of the requested interval was com- 
pleted 

User response: None 

Either the header record arrays 
or the data record arrays were 
filled before reading of the re- 
quested interval was completed 

User response: None 

The data interval specified by 
IRECl and IREC2 in NAMELIST 
MAIN (for a read request) or 
NFRAMl and NFRAM2 in NAMELIST 
MAIN (for a processing request) is 
invalid 

For a read request: 

IRECl > IREC2 or 
IRECl ^ 0 

For specifying data to be passed 
to OASYS: 

NFRAMK 1 or 
NFRAMl > NFRAM2 or ‘ 

NFRAM2 > the number of frames 
I of data currently in core 

^ NFRAM2 - NFRAMl >200 

User response; Correct the invalid 
parameter, and repeat the request 

U904 NO DATA REC- Following processing of the read 

ORDS IN CORE request, there are still no data 

records in core because tlie user 
requested reading of header rec- 
ords only, or the specified data in- 
terval is not on the data set 

User response: Read some data 
records into core by specifying the 
correct parameters 


U901 END OF FILE ON 
OABIAS DATA SET 


U902 ARRAYS FILLED 


U903 INVALID DATA 
INTERVAL SELECTED 
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Message 


Description 


U905 ERROR RETURN xxx 
FROM ODAP 


U906 START TIME = 
XXXXXX.XXXXXX 


U907 FIRST FRAME NO. 
TO BE PROCESSED 
MUST BE 1 IF EPHEM 
VECTORS SAVED 


U908 TO SAVE EPHEM 
VECTORS SET INPUT 
DATA RECORDS MASTER 
NO. .LE. 200 


U909 EPHEM VECTORS 
NOT PRESERVED 


Subroutine ODAP returned before com- 
pleting simulation of the requested 
block, probably because invalid pa- 
rameters were input to NAMELIST 
LIST. See the description of 
error code xxx in Section 7. 16.2 

User response: Repeat tlie simula- 
tion request, and specify valid pa- 
rameters in NAMELIST LIST 

This message always appears on plot 
displays of data versus time. START 
TIME is the time of the first unflagged 
frame of data (expressed as 
YYMMDD. HHMMSS) and is the zero 
point for the X-axis coordinate 

User response: None 

An incorrect inters'- al for passing data 
to OASYS has been requested. If 
ephemeris vectors are to be saved with 
the data, the data interval passed to 
OASYS must be 1 to N. 

N^ 200 

User response: Set NFRAMl = 1 on 
’'Select Data for OASYS" display or 
lEFMFG = 0 on "Options for Reading 
Data" display 

An attempt was made to save the 
ephemeris vectors passed with the 
data with an invalid "Input Data Rec- 
ords" master number 

User response: Set lEFMFG = 0 in 
"Options for Reading Data" display or 
return to array allocation sizes dis- 
play (key 0) and set "Input Data Rec- 
ords" master number to 200 or less 

An attempt was made to use the 
ephemeris information on the OABIAS 
file when that information was not 
preserved 


7-144 

PREL!/*.1INARY DRAFT 



PRELiA/HMARY DRAFT 


Message 


U920 END OF FILE 
NAMELIST/LIST/ 


U930 END OF FILE 
NAME LIST/OPMANl/ 


U947 ERROR 
RETURN XXX FROM 
EPHEMG 


Description 


User response: Set lEPHEM = 0 
on "Select Data for OASYS" display 
or re-read data with lEFMFG = 1 . 

(NOTE: See restrictions on pre- 
serving ephemeris vectors in error 
messages U907 and U908) 

An end-of-file was encountered while 
reading NAMELIST LIST in ODAPIN. 

The NAMELIST cards are probably 
missing or in the wrong order 

User response: Supply the necessary 
NAMELIST parameters via the display 

Ah end-of-file was encountered while 
reading NAMELIST OPMANl in 
OPINIT. See description of U920, above 

Subroutine EPHEMG returned an error 
code while attempting to access 
ephemeris data for the plot of pre- 
dicted and observed rotation angles. 

The error codes from EPHEMG are 
as follows: 

= 0, normal return 

= 1, time before start of tape (DTAPRE) 

= 2, time after end of tape (DTAPRE) 

= 3, I/O error on tape (DTAPRE) 

= 4, interpolation error (DTAPRE) 

= 5, too many tapes (DTAPRE) 

= 6, time before start of tape (RJPLT) 

= 7, time after end of tape (RJPLT) 

= 8, requested body not on file (SUNRD) 

= 9, I/O error (SUNRD) 

= 10, orbit level not found (GETHDR) 

= 11, I/O error reading header (GETHDR) 
= 12, time not on file (GETVCT) 

= 13, I/O error reading data record 
(GETVCT) 

= 14, input code is invalid (ISUN, ISPC, 
or IMOON, in COMMON/ORBITl/) 
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Message 


Description 


= 15, central body flag returned by 
GETHDR is not 1 or 2 
= 16, exceeded maximum iterations in 
solution of Kepler's equation 
(ORBGEN) 

User response: Depress key 2 to return 
to the NAMELIST OPMANl display 
and correct the ephemeris option 


U948 SPACECRAFT 
WITHIN CENTRAL 
BODY-PLOTOC 


The ephemeris data places the space- 
craft within the central body 

User response: Correct the ephemeris 
options on tlie NAMELIST OPMANl 
display or correct the effective cen- 
tral body radius on the display of 
options for data prediction 


U949 TO = XXXXXX.XXXXXX 
STDU = X.XXX, X.XXX, 
X.XXX. RESD = X.XXX, 
X.XXX, X.XXX 


This message always appears on the 
plot of predicted and obser\’’ed rotation 
angles versus time. TO is the start 
time of the prediction in the form 
YYMMDD. HHMMSS. STDV and RESD 
are the standard deviation and mean 
residual, respectively, of the predicted 
rotation angles minus the observed 
rotation angles. The three values are 
for plots 1, 2, and 3, respectively 


User response: None 


U950 PLOTS CALL FOR 
DIFFERENT INDE- 
PENDENT VARIABLE 


An attempt was made to obtain plots 
with different X-axis coordinates, such 
as plotting at data times combined with 
equal time increment plotting 


User response; Regenerate all plots 
by settir^ ISTATE to 4 or less for 
each plot 


U951 INDEPENDENT 
VARIABLES DISAGREE 
AT XXXTH ELEMENT 


Two plots have different X-coordinate 
elements for the XXXth point 

User response; Check source of X- 
array for probable error 



Message 


Description 


U952 TO = XXXXXX.XXXXXX This message always appears on plot 
STDV = X.XXX, X.XXX, of predicted and observed Eartli width 

X.XXX, RESD = X.XXX, versus time. Same interpretation as 

X.XXX, X.XXX U949 above with Earth width substituted 

for rotation angle 

User response: None 


U953 ILLEGAL USE 
STATE AT PLOT X HAS 
BEEN CORRECTED. 
PROCEED. 


An invalid use state was requested for 
plot X. The use state has been changed 
to reflect the current contents of that 
plot array 


User response: Proceed 


MESAGE MESSAGES: 


END OF FILE 
NAIME LIST/MAIN/ 


TERMINATION 

REQUESTED 


HEADER RECORD AT 
ODD POSITION ON DATA 
SET 


ODD NUMBER OF 
RECORDS ON DATA SET 


An end of file was encountered while 
reading NAMELIST MAIN in OADRIV. 
The NAMELIST cards for MAIN are 
probably missing 

User response: Continue. Supply the 
necessary parameters via the display. 
All other NAMELISTS in the deck will 
have been skipped over also 

The operator has requested termi- 
nation of the job 

User response: If termination is 
desired, continue. If termination 
is not desired, depress key 1 

A header record was encountered 
on the OABIAS data set as the sec- 
ond record of a pair of records 

User response: Continue. The 
first record in the pair will be 
ignored. Correct the program 
which operated the OABIAS data 
set to generate a valid data set 

An end-of-file was encountered 
while attempting to read the second 
record in a pair of records 
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Message 


Descriptioa 


V 


INVALID RECORD. 
NSAMP . LE. 0 


DATA RECORD 
PRECEEDS RECORD 
HEADER 


ARRAYS NOT ALLOCATED 
FOR AE COPY 


OABIAS 

AECOPY IN PROGRESS 
COPYING BLOCK 
NO. XX 

OABIAS 

DATA READ IN 
PROGRESS 

OABIAS 

OASYS PROCESSING 
IN PROGRESS 


User response: Continue. Correct 
the program which generated the 
OABIAS data set to generate a 
valid data set 

The parameter NSAMP on the 
OABIAS data set is ^ 0 

User response: Same as for pre- 
vious message 

The first record on the OABIAS 
data set is not a header record 

User response: Correct proc- 
essing of data records is not possi- 
ble unless the header record is 
available. Correct the program 
which generated the OABIAS data 
set to create a valid data set 

One or more of the arrays passed 
to AECOPY has not been allo- 
cated by the GESS Executive (i, e, , 
the array allocation size is zero). 
Return to die array allocation 
sizes display via key 0 and correct 
the array allocation sizes 

Copying of the AE data set is pro- 
ceeding. The XXth block has been 
processed 

User response: None 

Reading of the data into core is pro- 
ceeding 

User response: None 

Processing of the data by OASYS 
is proceeding 

User response: None 
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Message 


Description 


OABIAS processing is proceeding. 
After XX frames of data have been 
processed, the uncertainties in 
right ascension and declination 
are as given 

User response: None 

PLOTOC data prediction is proceed- 
ing. Plot number X has been 
generated 

User response: None 

7.16.2 Printed Messages From the ODAP Subsystem 

All error messages are written on FORTRAN unit number 6, as shown in 
Figure 7-43. Each error message provides the following information: 

Heading Description 


OABIAS 

BIAS DETERMINATION 
IN PROGRESS 
XX FRAMES PROCESSED 
XXX. XX UNCERTAINTY 
IN ALPHA = XXX. XX DEG. 
UNCERTAINTY IN DELTA = 
XXX. XX DEG 

OABIAS 

DATA PREDICTION IN 
PROGRESS PLOT NO. X 


FRAME NUMBER 


ERROR CODE 
CENTRAL BODY 


MESSAGE 


A number corresponding to the frame 
number which appears in the detailed 
printout. With this number, a user 
can determine at what time in the 
simulation the error occurred. If the 
frame number = 0, then the error 
occurred in the initialization phase of 
the program 

A number identifying the error type 

Either EARTH, MOON, or blanks. If 
this error occurred while processing 
for a particular central body, then the 
central body is indicated. If this 
column is blank, it does not apply 

A message describing the error 


The following is a list of all error messages produced by the program. Those 
messages with an error code less than 100 are warning messages only; an 
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Figure 7-43. Error Message From ODAP 



appropriate fixup is taken, and execution of the problem continues. Those 
messages with an error code greater tlian or equal to 100 result in termination 
of tlie simulation, with the additional message: 

EXECUTION OF THIS PROBLEM TERMINATED DUE TO SEVERE 
ERROR 

Control then returns to OADRIV, and message U905 is displayed. (See Sec- 


tion 7. 16. 1, above. ) 

Error 

Code Message 


1 CONES DO NOT INTER- 

SECT--ODAP1 


2 DIHEDRAL ANGLE 

UNDEFINED— ODAPl 


3 MAX ITERATIONS IN 

TERMIN— DELPHI 
TOO SMALL 


4 POINT DOES NOT LIE 

ON CENTRAL BODY— 
TERMIN 


Meaning 

The CONES routine indicates that 
the sensor scan does not intersect 
the central body, although a pre- 
vious computation indicated that the 
scan should intersect the central 
body. This error may occasionally 
result due to truncation error with 
REAL* 4 arithmetic 

User response: None 

The spin axis is parallel to the 
sun vector or the sensor mount- 
ing angle is zero or. 180 degrees, 
so rotation angles cannot be 
computed. 

User response: None 

The terminator search exceeded 
1000 steps. This error cannot 
occur if the default value of 
DELPHI is used 

User response; Increase DELPHI 

In the search for a terminator 
crossing, a point in the scan failed 
to lie on the central body, although 
a previous computation indicated 
that all points in this segment of 
the scan should lie on the central 



Message 


MAX ITERA.TIONS IN 
SEARCH--PHITOL TOO 
SMALL 


POINT DOES NOT LIE 
ON CENTRAL BODY- 
SEARCH 

ORDER OF CROSSINGS 
UNDEFINED 


INVALID INPUT TO 
SUBROUTINE xxxxxx 


INVALID INPUT TO 
SUBROUTINE xxxxxx 


Meaning 

body. This erroi* will occur oc- 
casionally due to the inaccuracy 
of REAL *4 arithmetic 

User response: None 

The half-interval search for a 
terminator crossing exceeded 
30 steps. This cannot occur if 
the default value of PHITOL is 
used 

User response: Increase PHITOL 
See Error Code 4 


The vector to the horizon crossing 
point is parallel to the nadir vector 
or the spin axis; the order of the 
horizon ci'osssings cannot be de- 
termined 

User response: None 

xxxxxx will be the name of one of 
the uncertainty routines, UNCERT, 
UNCDBL, or UNCDH. No solution 
was obtained for attitude, using the 
imbiased simulated data. This will 
occasionally result due to trunca- 
tion error with REAL *4 arithmetic 

User response: None 

xxxxxx will be the name of one of 
the uncertainty routines UNCERT 
or UNCDBL. The spacecraft po- 
sition obtained from EPHEMG is 
within the central body 

User response: Correct the orbit 
information being used by EPHEMG. 
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Code 


Meaning 
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Message 


31-46 

INVALID INPUT TO 
' SUBROUTINE xxxxxx 

70 

EOF ON ATTITUDE 
TAPE— NORMAL 
TERMINATION 

80 

EXCEEDED MAX ITER- 
ATIONS IN ODAP 

101-116 

ERROR IN EPHEM 
ROUTINE 

130 

STOP TIME COMES 
BEFORE START TIME 

150 

SPACECRAFT IS 
WITHIN EARTH 

160 

SPACECRAFT IS 
WITHIN MOON 

200 

EOF ON ATTITUDE 
TAPE WHILE READING 
HEADER RECORD 


xxxxxx will be tlie name of one of 
the uncertaintj’- routines, UNCERT, 
UNCDBL, or UNCDH. An error 
return was encountered from 
EPHEMG, witli error code 1-16, 
respectively. See description of 
EPHEMG 

An end-of-file was encountered 
while attempting to read a data 
record on the attitude tape 

User response: None 

The maximum number of iterations 
specified by the user to be used in 
the iterative process to correct for 
orbital motion has been exceeded 

User response: None 

Eiror return 1-16 respectively 
from subroutine EPHEMG. 

See description of EPHEMG 

The specified value of TF is earlier 
than the specified value of TO. 
(NAMELIST LIST) ‘ 

User response: Do not specify TF 

The spacecraft coordinates place it 
within the Earth 

User response: Correct the orbital 
elements or orbit tape 

See Error Code 150 

An end-of-file occurs within the 
header record, or there is an il- 
legal, missing, or DUMMY DD 
card for the attitude tape 

User response: Check the attitude 
tape DD card 
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Error 

Code Message 

210 ATTITUDE TAPE DOES 

NOT COVER ANY OF 
THE DESIRED TmE 
INTERVAL 


220 TIME REVERSAL- 

BACKSPACING IS NOT 
PROVIDED FOR ATTI- 
TUDE TAPE 


Meaning 

The range covered bj’- the attitude 
tape, as indicated on the header 
record, does not include any of the 
requested time interval, as defined 
by TO and TF (NAMELIST LIST); 
or the attitude tape has fewer than 
two data records; or the header 
record does not correctly indicate 
the start and stop times of the data 
records 

User response : Check the printout 
of the header record to see if the 
tape covered the requested interval 

Interpolation has been specified 
(INTERP = 1, NAMELIST LIST), 
and the requested time is earlier 
than the time interval covered by 
the two records which have been 
saved. This error should not occur 

User response: Correct the program 
so that the error does not occur 


The following unnumbered message is provided: 


Message 


Meaning 


END OF FILE ON 
UNIT 5— NAMELIST/ 
LIST/NOT FOUND 


An end-of-file was encountered on 
FORTRAN data set reference num- 
ber 5, while attempting to read 
NAMELIST LIST for the first 
problem. This message is always 
accompanied by the graphics mes- 
sage U920 described in Section 
7.16.1 

User response: Supply the 
NAMELIST parameters via the 
display device, if this is a graphics 
run 
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7.16.3 Printed Messages From the OASYS Subsystem 


All error messages from the OASYS Subsystem are written on FORTRAN unit 
The following messages are provided: 


Description 


number 6. 

Message 


A PRIORI ATTITUDE 
WAS USED TO RESOLVE 
AMBIGUITIES IN SPINAV 


NO A PRIORI ATTITUDE 
AVAILABLE —AMBIGUITIES 
CANNOT BE RESOLVED IN 
SPINAV 


SPINAV DIVERGED— NO 
AVERAGES AVAILABLE 


The a priori attitude provided was 
used to resolve ambiguities, either 
because the deterministic logic 
failed or because lAPIOR = 8 
(NAMELIST OPMANl) 

User response: Warning message 
only. The correctness of the re- 
sults may depend on the accuracy 
of the a priori attitude 

The search technique for resolving 
ambiguities failed to yield any set 
of selections which contained more 
than one useful attitude 

User response: Increase the block 
size and/or the time range of the 
block. Increase values of SPNSIG 
and SPNTOL. Try using lAPIOR 
= 4 and providing an a priori atti- 
tude (NAMELIST OPMANl) 

In the process of iterating to im- 
prove the selections for ambiguous 
pairs of attitudes, every attitude in 
the block was rejected as spurious 

User response: Increase rejection 
tolerances SPINSIG and SPNTOL. 
Try using lAPIOR = 8. Increase 
the block size and/or the time 
range of the block (NAMELIST 
OPMANl) 
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Message 


Description 


SPINAV FAILED TO 
CONVERGE— RESULTS 
ARE GIVEN FOR LAST 
ITERATION 


FITTING OF SUN DATA 
FAILED DUE TO ZERO 
DETERMINANT IN BLOCK 

XX 


FAULTY POINTER IN- 
FORMATION IN SUNPRO 


The process for selecting the cor- 
rect attitude from ambiguous pairs 
failed to converge within IT MAX 
iterations. The selections which 
are given for the last iteration are 
not necessarily consistent with the 
average obtained using those selec- 
tions 

User response: Increase IT MAX. 
Try using lAPIOR = 8 and specify- 
ing an initial attitude estimate. 
Increase the block size and/or the 
time range of the block (NAMELIST 
OPMANl) 

In attempting to perform a linear 
fit to the Sun angles as a function 
of time, a zero determinant was 
obtained 

User response; Set ISNPRO=l 
(NAMELIST OPMANl). Check 
Sun angles and Sun times for ex- 
treme erroneous values, and flag 
these frames 

The pointer information' in 
COMMON/POINT/ is invalid. 

This error should not occur 

User response: Correct the pro- 
gram so that the erx'or does not 
occur 


7.16.4 Printed Error Messages From OABIAS 

Error messages from the OABIAS Subsystem are written by subroutine BIASER 
on the FORTRAN unit number specified by ERUNIT. All messages start with 
the frame number and the name of the subroutine in which the error occurred. 
This information is followed by the error message. 
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In many cases, the same error message may be printed by several subroutines. 
Therefore the error messages are given below in the order in which they appear 
in the OABIAS error summary table. (See Section 7. 14. 16 for an example). 
Table 7-1 provides the correlation between the error summary table and the 
error message descriptions. 

The subroutine BIASER prints the appropriate error message based on the 
variable lERR which it receives via the calling sequence. lERR is a four-digit 
number. The first two digits indicate which subroutine has encountered an er- 
ror. The next two digits are the index for the error which was encountered 
within the subroutine, i.e. , the first, second, or third. The numbers which 
identify the subroutines are given in Table 7-1. The indices utilized by each 
subroutine are indicated in Table 7-1 by a reference to an error message num- 
ber. A blank space indicates no error message for that index. The values of 
lERR which can cause an error message to be written are listed in die message 
descriptions. 

In some cases, the message contains parameters which may be useful in de- 
termining the exact cause of the error. For example, 

***** ERROR IN FRAME 32 FROM DIAMOND - DATA REJECT. - REJOPT = 1, 
REJ = 8 

where REJOPT is the data rejection option in NAMELIST BIASNL, and REJ 
is the data rejection flag for frame 32 in COMMON/RE JCOM/. This message 
was caused by lERR = 0511. In the error messages given below, the characters 
'XXXX’ will be used to indicate the printing of a variable in that place in the mes 
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Table 7-1. OABIAS Error-Message Error-Table Correlation Matrix 


imber 

Subroutine 

Name 

1 

AMATRX 

2 

APARTS 

3 

BIASER 

4 

DIAFUN 

5 

DIAMOD 

6 

FRAPRO 

7 

INITIL 

8 

LCOMP 

9 

LNFUN 

10 

LNMOD 

11 

LPARTS 

12 

LRFUN 

13 

LRMOD 

14 

OABIAS 

15 

PRINTB 

16 

PSIPHA 

17 

RECURI 

18 

ROTATE 

19 

SANFUN 

20 

SANMOD 

21 

STMFUN 

22 

STMMOD 

23 

DMFUN 

24 

DHMOD 

25 

SCBFUN 

26 

SCBMOD 


Error Error Error 
1 2 3 


2 

3 

4 

8 

9 

10 


14 

15 


8 

11 

12 

18 

19 


16 

17 


8 

11 

12 


20 

9 10 12 

9 10 12 

21 

8 11 12 

16 22 
8 11 12 


Error Error Error 
4 5 6 


5 6 7 

11 12 13 


16 17 

13 


13 

17 
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The OABIAS error messages fall into two classes: warnings and information. 
The warning error messages are associated with computational failures. The 
program action is to bypass tlie observation wliich caused the error. If any 
user response is required, it is indicated after the description of the error 
message. The information messages record data rejection. If an observation 
has been flagged, it is not processed. No user response is required. 

Warning error messages are as follows: 

Message 


Number 

lERR 

Message 

Description 

2 

0401 

DOT PRODUCT OF SUN 
VECTOR AND SPIN 
VECTOR IS 1 or -1 

The Sun vector and spin 
axis vector are parallel 
or antiparallel (the com- 
puted Sun angle is zero or 
180 degrees) 

3 

0402 

SENSOR MOUNTING 
ANGLE =0.0 

Computed sensor mounting 
angle is zero 

4 

0403 

ARGUMENT OF SQRT, 
DELR, IS NEGATIVE 

Argument of square root 
is negative in computing 
the variable DELR in 
subroutine DIAFUN 

5 

0404 

ARGUMENT OF SQRT, 
DELP, IS NEGATIVE 

Argument of square root 
is negative in computing 
the variable DELP in 
subroutine DIAFUN 

6 

0405 

TOP**2/BOT** = -1, 
TOP**2 = XXXX 
BOT**2 = XXXX 

A denominator in comput- 
ing the single horizon 
dihedral angle model func- 
tion is zero 

7 

0406 

ARGUMENT OF ARCSIN, 
(TE, TH) .GT. 1.0 

The argument of arcsine 
of quantities in computing 
the single horizon dihedral 
angle model function is 
greater than 1. 0 
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Message 


Number 

lERR 

Message 

Description 

8 

0501 

EARTH SIGHTING TIME 

The observation is not 


1001 

IS NEGATWE. EINOUT = 1 

processed because the 


1301 

(0 = EARTH-OUT) 

central body sighting 


2401 

(1 = EARTH-m) 

time is negative (data 


2601 

INITIM = XXXX 

rejection flag). EINOUT 



EARTH-IN TIME = XX"XX 

indicates the type of 



EARTH-OUT ITME = XXXX 

triggering. INITIM is the 
block reference time in 
from September 1, 1957 

9 

0502 

SUNSIGHTING TIME IS 



2001 

NEGATIVE, SUNTIM = XXXX 

The observation is not 


2201 


processed because the 
sun sighting time is neg- 
ative (data rejection 
flag) 

10 

0503 

SUNANG (NOBS) = XXXX 

Sun angle observation is 


2002 

IS LT 0 OR GT PI 

either negative or greater 


2202 


than 7T . Error in the 
telemetry 

11 

0504 

DATA REJECT - 

Data is rejected because 


1002 

REJOPT = XXXX 

it was rejected by the 


1302 

REJ = XXXX 

corresponding model in 


2402 


OASYS 


2602 


« 

12 

0505 

BAD ESTIMATE, FIX- 

The sum of S^^^ and $ 2 ^ 


1003 

UP BY NORMALIZA- 

exceeds unity. S^ and 


1303 

TION 

S 2 „are adjusted so that 


2003 

SI = XXXX, S2 = XXXX 

Si^ + S2^ < 1 


2203 



2403 

2603 


User response: If this 
occurs frequently, use 
smaller PO 

13 

0506 

L VECTOR NOT OB- 

The current estimate of 



TAINED IN LCOMP 

the attitude and biases 


yields a scan which misses 
the central body. The 
model calling subroutine 
LCOMP is not processed 
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Message 

Number 


lERR 


Message 


Description 


14 0801 D = XXXX, XXXX, XXXX 

LT TOLDEN = XXXX 


15 0802 ARGUMENT OF SQRT = 

XXXX 

NEGATIVE IN QUAD- 
RATIC EQUATION 


16 0804 POSITION VECTOR IS 

1201 ZERO 
2501 


Denominator (D) in com- 
puting the horizon vector, 
Lg , is less than a speci- 
fied tolerance (TOLDEN 
in NAMEUST BLISNL) 

A complex solution for 
the horizon vector, Lr , 
was obtained. This im- 
plies that the current 
estimate of the attitude 
and biases does not yield 
a scan of the central body 

The position vector is 
zero. Error in the 
telemetry 


17 0805 ERROR IN EPHEm^ 

1202 
2504 


Observation is not proc- 
essed because error oc- 
curs in subroutine EPHEMV 


18 1101 THE DETERMINANT 

XXXX LT THE SING- 
ULARITY TOLER- 
ANCE = XXXX 


The matrix to be inverted 
in. computing the partial 
derivatives of the horizon 
vector, lijj , is singular 

User response: If this 
occurs frequently, de- 
crease the tolerance. 

This error is inevitable 
if the nadir angle is close 
to 180 degrees or the cen- 
tral body is extremely 
small. If this is the case, 
use the nadir vector pro- 
jection model and/or scan 
mid-time dihedral angle 
model only 
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Message 

Number 

lERR 

Message 

Description 

19 

1102 

ERROR IN MINVM 

Observation is not proc- 
essed because error 
occurs in subroutine 
MINVM 

20 

1901 

CHI GE 1, 

CHI = XXXX 
COS (X (10)) = XXXX 

Tlie cosine of the computed 
Sun angle exceeds 1 . 0 

21 

2301 

BOT = 0.0 

The denominator BOT in 
subroutine DHFUN is 
zero 

22 

2503 

ABS (RCHI) . GT. 1. 0 

Cosine of the computed 
mounting angle is greater 
than 1. 0 


7. 17 JOB CONTROL LANGUAGE 

Figure 7-44 shows the JCL required to load the required source modules from 
tape to disk, perform GESS table generation, compile the source modules, 
perform linkage editing with overlay, and execute tlie MSAD/OABIAS System. 
The source library on file 1 of the tape contains OADRIV and its associated 
modules and routines for the ODAP and OASYS Subsystems which were modi- 
fied for use with MSAD/OABIAS. File 2 contains the modules ib the OABIAS 
Subsystem. All other required routines come from the non-executable load 
module libraries referenced in the link step. 

Note that in the linlc step subroutine RJPLT, which reads the JPL Ltmar and 
Solar Ephemeris File, has been replaced by a dummy routine in order to save 
about 12K bjrtes of core storage. Thus the option to read the JPL Ephemeris 
File will not be available; subroutine SUNRD can be used to obtain lunar and 
solar ephemeris. If the user requires subroutine RJPLT then the cards with 
sequence numbers 810 and 820 must be deleted. 


7-162 

PRELIMINARY DRAFT 



OABIAS COMPILE 


//7RJRWOAR jne ir.Hft001852C,T,G00529,005008).PFF 
//•MAIN LINFS-lOO 
//• 

//• RFLOAD SOURCE CODE FROM TAPE TO DISK 

//POSURl EXFC Pr.«*POSUR ,REGI0N*080K 
RFlOAn F=T1 
RELOAD E=T?»T=01 

//Tl DO niSP=WPASSl,LAheL = ll*BLP» fUN ITa?400, OSN»T APEOl , VOL*SER*3671 5, 
// 0CR«IRFCEM.FH,LRFCL=HO,8LKSrze=7?RO,DEN=3J u ^ 

'<P*55I ,L*l>eL = l2,f>LP),vnL = HtF=«.I),DCB = ».Tl,r)SN«».Tl 
n ISP* (NEW, Pass I , space = icyl. is, i .si i ,uni t= 231 A,nsN*ccB, 

// DCfi=»,Tl 

D ISP. (NEW, PASS I .space *(CVl, 15,1,51 ), UNI T*231 A, nSN*Ct A, 

/ t HL n * • « I ? 

//SYSPRINT DD SYSOUT«A 
//• 

//• ASSEMBLE MSAO TABLES. 

//• 

//ASMG EXEC ASMG»PABM=*N0X,LINECNT=76* ,REGION*300K 
//SYSLIB DD 

// no 

.. A, ^0 OISP=SHR«OSN=ATT1T.GESSMAC.ASM 
//SYSLIN DD UNI T= (0 1SK,4 I ,SPACe = ICYL* (5 » 1) ) 

//SYSPRINT no UNIT = (D1 Sk, 4) 

//SYSUT2 DD UN |T= (0 ISK ,4 1 
//SYSUT3 DD UNIT=(01SK.4| 

//SYSIN DD niSP=(nLD,PASSJ ,OSN*CCA(TABLES) 

//J EXECUTE MSAD TABLE GENERATION PROGRAM. 

AB FXFC L INKGO, PARMs »NCAL »l I ST, MAP * ♦ CONO«< 16.LT ) , R6GI ON. G0= 1 50K 
//SYSLIB DD DISP = SHR,DSN=ATTlT.r.eSS.LOAO lt^coIun.ou 

//SYSLIN DD DISP* {DLO.PASS) 

//OBJECT DD • 

INCLUDE SYSLIB ITAflGEN) 

ENTRY TARGFN 

nSNseCNONR6S,UNITsDISKtOISP = INEW,PASS) , 

J ,f>CBs(RkEMsUT.HLKS 126 *5000 J 
//GO.X OD nSN*tUnOMOO,OlSP*inLr),OELETEI 

//I TO SeaJENTIAL 

//SFOl EXEC SEOL.DSN.'CCA' 

AECOPY 

blkplt 

flinCKDAT 

DISKS 

OADRIV 

DANA IN 

oaplot 

ODAP 
ODAPIN 
DDAP2 
OPIN 1 T 
DPMA IN 
PLDTOC 
SETELG 
TAPMED 
TC0N20 
TC0N40 
//• 

//• compile CCA 
//• 

//COMPLl EXEC F0RTRANH,PARM=»0PT=2,MAP.XR6F » ,REGI0N»300K 
//SYSLIN DD UN! T* <0 1SK.4 J ,SPACe«(CYL, 15, 1 ) ) 

//SYSPRINT DO UN1T= (01 SKjA) i. f i p , i # # 

//SYSIN DO OSN-CCSOURCe ,DISP«(0L0, DELETE) 

//• 

//•CONVERT CCB TO SEQUENTIAL 

//SF07 FXFC SFOL,OSNa*CCB» 

AMATRX 

APARTS 

6TASER 

BLOCKOAT 

DHFUN 

DHMOO 

OIAFUN 

OIAMOO 

EPHEMG 

EPMFMV 

FRAPRO 

GFTMOR 


oonoo 100 
00000200 
ooono‘^00 

O0DO0400 
00000500 
00000600 
00000700 
OOOOOflOO 
00000900 
OOOOIOOO 
00001 100 
00001 200 
00001 300 
00001400 
0000 I 500 
00001600 
00001700 
onnoiBOO 
OOOOIOOO 
00002000 
00002100 
00002200 
00002300 
00002400 
00002500 
00002600 
00002700 
00002800 
00002900 
00003000 
00003100 
00003200 
00003300 
00003400 
00003500 
00003600 
00003700 
00003B00 
00003900 
00004000 
0O005 100 
00004200 
00004300 
00004400 


000O4S00 

00004600 

00004700 

00004B00 

00004900 

00005000 

00005100 

00005200 

00005300 

00005400 

00005500 


Figure 7-44. 
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preliminaW draft 




HFMfTR 

INITIL 

tCOMP 

LNFON 

LNMnn 

LPARTS 

LRFUN 

LRMOn 

MINVM 

OABTAS 
ORfiGFN 
PR INTfi 
PR INTP 
PS IPHA 
RFCUK I 
rnTATE 

SANFUM 
SANMnn 
SCR FUN 
SCRHnn 
stmfun 
STMMD n 
VFCROT 
//• 

//* CnMPILE CCB 
//* 

//COMPL2 FxeC FORTRANHtPARM* 'nPTs?,HAP,XReF ' tREGI0N»300K 

//SYSLiN nn uNiT=(oiSK.A),sPAce=(CYL,(i,i)) 

//SYSPRiNT nn uNns(niSK,4) 

//SYSIN nn nsN=ccsnuRCF.niSP*(nLf>tnFLFTF) 

//* 

//• LINK FniT WITH OVFRLAY. 

//* 

L INK,CONn=< I6.LTI ♦NBLK=2SO,PARM=(LETtLIST,MAP,OVLY» 


//LINK FXFC 
//SYSLIB 
ft 
// 
ft 


L «CUrMi;=( 1 ^ t L I > t NO L K= £*)(>• KAKM- 

00 OlSPsSHK.DSNsATTIT.GESS.LOAD 
on D!SP=SHR,OSN«ATT I T , AOOA • OA SYSA, L I B. OBJ 
no niSP=SHR,DSN=ATT IT.OPRLIB.OBJ 
, on DISP=SHR,nSN=ATTIT,RAEB.LOAa 

ft no DSN=SYS 1 .FORTL IB ,n 1 SP=SHR 

if no osnssys 2 .gsfcl iB,nisp=SHR 
// on nsN=sYSi .sspAK,n !SP=SHR 
//SYSLMon nn 0 SN=££LnnMnD( 0 ABi 4 S) 

//SYSuTi nn uNn = (oiSK,3) 

//SYSLIN nn # 

INCLUOH AFl IH( IHCFRRM) 

CHANGE IHCFRRM(FnRTC-RRM) 

INCLUOF SYSl 1 H( IWCFHRM) 

RFPLACF NONRESl .NnNRES 2 
INCLUOF RFS 

INCLUOF SYSL IB IMSAOMA INJ 
INCLUOF SYSLIRIMSAnBLOKl 
INCLUnF SYSL IR {PLOTS ) 
lNCLt»DF SYSUR(GCOMM) 

INCLUOF SYSL IB (fancKon 
CHANGE BYPASS(RJPLT) 

INCLUnf SYSL 18 (BYPASS! 

CHANGE BYPASSIGCONVC) 

INCLUOF SYSLIBfBYPASS) 

INCLUOF AFLIH(ZBPABA) 

FNTRY MSAOMAIN 
OVERLAY A 

^^|NSERT^INIT,R0XSTP,R0RTCM,H$F0RM,FMTC0M 

INSERT OAORIVtOPMA IN 
OVERLAY B 

INSERT AFCnPY.WRTnuT.WRTRUF 

INSERT R I TS4.B ITSfl, GAUSS *NOlS£4,NOISEfl,RAN0U*TCnN20,TC0N40 

UVFRLAY R 

INSERT nOAP,ATOATA,FLAG* ICH AR « OOAP IN « 00 AP 2 tOHEGA «R£PORT 
OVERLAY R 

Double ♦OHCOM.OBLCOM, EXTRA, OUT INF 

INSERT niSKS*0P06TS,TERCMK«0BLCRS 
OVERLAY E 

^^^NS6RT^CH00SE,ERR0R1,0P!NIT,SP1NAV,$PNAV1,SUNPR0 

INSERT BLKPLT, 0 AGRAF,CURTIM 
overlay ft 

INSERT nABIAS»PRlNTB,BIAS 6 R,OISPLY*V 6 CROT 
INSERT EPHROT,RESIOU,F I LPAR , 6 RRCOM, VCOHON 

overlay C 

^V^NSERTcINITIL,PSIPHA,ROTATE,nAPLOT,PR|NTP 

INSERT AMATRX.APARTS,0UFUN,0IAM0D*FRAPR0,LNFUN 


00005600 
00005700 
OOOOSftOO 
00005P00 
00006000 
00006100 
00006200 
00006300 
00006400 
00006500 
00006600 
00006700 
00006BOO 
00006900 
00007000 
00007100 
00007200 
00007300 
00007400 
00007500 
00007600 
00007700 
00007BOO 
00007P00 
00008000 
00008100 
OOOOP200 
0OO0B300 
0000B400 
00008500 
00008600 
00008700 
OOOOHCOO 
00008900 
00009000 
00009100 
00009200 
00009300 
00009400 
00009500 
000096(10 
00009700 
00009POO 
00009900 
00010000 
00010100 
00010200 
00010300 
00010400 
000 10500 
00010600 
00010700 
OOOIOBOO 
00010900 
0001 1000 
00011 100 
00011200 
00011300 
00011400 
00011500 
00011600 


Figure 7-44. JCL to Compile, Link, and Execute 
MSAD/OABIAS (2 of 4) 
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INSERT LPRRTS.LREUN.LRNOD.RECURl 
INSERT AB.EPHPMV 
OVERLAY HSAOlREGinN) 

IKPert 

INSf RT *C PF XFC . CP IN I T, DA EXFC , DA I N I T , C P 1 BY I 

v^!=r?AjA ,nAC0Nv,n isPcnLS.niSPnATAtDiSPPARH 

INSERT Pl.TE XFC t PUCOK.GENGnP.VPTFRM 
INSERT EXPAND, GFTPUT ♦ IDF NT IF Y. SCOP LT 

nVFRUAY MSAD ’ 

INSERT nHMnD,SCBMnn,SANMOD,STMHOO, stmfun,lcomp,minvh 
^ INSERT nHFUN,SCHFUN,SANFUN,l.NMOO »LCunr,m,Nvn 

nVFRlAY MSAD 

INSERT graph, scale, SETFLG 
nVERLAY MSAD 
INSERT ATTDFT 


INSFRT ML ir,HT , L I GHT , DD A PER ,00 AP I .DP S I M 
INSERT DROfR, SEARCH, SPHERE .SUNVIS 


nn DISP=SHR,DSN=A7TIT.AOOA.tOAO 

LOAD MODULE, AND TABLES FROM DISK TO TAPE 


INSERT TF RMIN.uNCERT ,VPHASF 
INSERT DRLU,DHCnMP,UNCnH,UNCDBL 
//RES DO DSN=£/.nHjMOD»n ISP = inL0,06LFTE) 

//AFL I H *- 

//* 

//* UNLOAD SOURCE CODE 
//» 

//PDSUR2 FXFC PGMsPOSUR ,REGION=80K,C0N0»(8,LT J 
UNLOAD FrDl,T*n,L 
UNLOAD F=D?,T=T?,L 
UNLOAD F=D^,T«T^,L 
UNLOAD F=D4,T*TA 

//Dl DO DISPsIOLD, DELETE), DSN=C£A 
//D 2 no ni SP= (OLD.DFl FTE ) ,DSNsCC8 
//D3 DO DlSP=(OLO,PASS),USN*ttLOUHUn 

//DA no DISPsIOLO.PASSI.OSMsLCNnMRES 

> tP^SS) , LABEL -(2,BLP),VOLeREF=*,Tl,OCB=*.Tl,DSNs*,Tl 
//T? nn DlSP=t,PASS),LABeLs(3,RLP),vnLsREFs*.nJoCB»*.Ti,DSNs*,Tl 
//s^SPRiNPnn'l?snuUk‘"'^‘-''‘’’»‘-'’'*’''’‘-=« • dc«=*-ti.dsn=».ti 
// 

// 

//* 

//• e U CUTE OABIAS 
//• 

' ♦REG10N»4 50K,C0ND*EV6N 
nSNsSYSl.UNKLlH,niSPsSHR,DC8sBuFN0*l 
// DO DSNsSYS2,GTS,OISPsSHR.OCBsBUFNOs1 ‘ 

It nn DISPsSHR,nSNeC£LODMt )0 

<< ^ 0ISP*SMR,OSSsATTlT,OPRLlB.LOAO 

//TJSSP ni) nSN=GIS.JCLO,OISP = SHR 

//PILP nn UNIT = Z3l*,01SP = SHR,vnL=REF = SYSl.LINKUB 
//SYSMTl nn HN|T=J J14,SP4CE=(TRK, ( 10, 51 ) 

//SYSHTr nn hnit»?31y,space= (trkI ( lols n 
//TJSC4 nn nSN=GTS.LnGREC.OISP=SHR 
//lEPRnR nn nsN=GTS.jCLO,niSP=SHR 

//SYSPkINT on dummy MSaU dumps 

//FT05F001 DO DISP=.SHK,nSfJ=r.KSMS.nATAUB.D4TAinPHAN) ,L46EL = ( INI 
//!!■ SETS. ASSUMING IPRINT=ERUNIT<6. 

MSAD, OAORIV, OOAP , OASYS, AND OABIAS. 

//•^PRINTOUT FROM OABIAS. IIPRINT-MI. 

//• PRINinuT FROM OABIAS. IIPRINT*?!. 

//cTS5coi?^’nrt^?Gc« (Earth as central body) ~ ■ - 



0001 1700 
000 1 IfiOO 
000 1 IQOO 
00012000 
00012100 
00012200 
00012'»00 
00012400 
00012S00 
0001 2^00 
0001^700 
00012ROO 
00012000 
00013000 
00013100 
00013200 
00013300 
0001 3400 
OOOISSCJO 
00013600 
00013700 
0001 3R00 
000 13900 
00014000 
00014100 
00014200 
00014300 
00014400 
00014500 
00014600 
00014700 
00014B00 
00014900 
00015000 
00015100 
00015200 
000153O0 
0001 5400 
00015500 
00015600 
00015700 
OOOl 5800 
00015900 


00016000 
00016100 
00016200 
00016300 
00016400 
00016500 
00016600 
00016700 
00016800 
00016900 
00017000 
00017100 
00017200 
00017300 
00017400 
00017500 
00017600 
00017700 
00017800 
00017900 
00018000 
00018100 
00018200 
00018300 
00018400 
00018500 
00018600 
00018700 
0OOIB8O0 
00018900 
00019000 
000 19 100 
00019200 
00019300 
00019400 
00019500 
00019600 
00019700 
00019800 


Figure 7-44. 
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PRELir/il 


JPL LUNAR AND SOLAR ePHEMERlS FILE (NRJPLT). 
GTDS SEOUFNTIAL SPC. EPMEMERIS FILE (NGTDS>, 
OODS SPACECRAFT EPHEHERIS FILE (N0R81). 

GTDS DIRECT ACCESS SPC. EPHEN6RIS FILE 
OAtilAS DATA SET. 


// UNIT«(niSK.A> .SPACe*<CYL. 10.2) ) 

PRINTOUT FROM OAHIAS. (IPRINT+S). 

//FTllFOOl no SYSOUT=A.DCB*IrECFH*VBA,LRECL= 137,BLKSI2E*137A»BUFNO*1) t 
// UN iTalOISK.A I . SPACF=I CYL ♦ ( 0,2) » 

//FT13F001 00 OlSPsSHRtOSNsATTIT.SSS I .MSAP.gftROR.DAT A 

//FTIAFCOI no niSP*SHrt.OSN*ATT IT. SHSB.O A81AS.DAT A 
//FT23FOni no UNlT = AFF=r,PAKOO 
//FT2PF00! DO DUMMY 
//FT2PFOOI no DUMMY 

//FT30F001 on dummy 

//FT31F001 00 dummy 
//FT 50F001 no DUMMY 
//* 

//• SCRATCH DATA SET FOR SIMULATED OABIAS DATA. 
ft* 

//FT51F001 00 UNIT* (OISK.A) ,SPAC6*(CYL, (2,2) ) ♦ 

// 0CB*IRECFMaFB,LR£CL»SO,BLKSI2Es20O0,ttUFNOal) 

//r.D. FT90F00I f)D 0SN = f.CNONRES ,0 1 SP« ( 01 0 . P ASS) MSAO TABLES (IFTABL). 

//FT91F001 no UN1T«2250 GRAPHICS DEVICE (IFTUBE). 

//FT04FD01 no DUMMY UDAP REPORT PRINTOUT (E^kRTH AS CENTRAL BODY). 

//FT9SF001 on DUMMY QDAP RFPORT PRINTOUT (MOON AS CENTRAL BODY). 

//GO. SYSUOUMP DO SYSOUT=A,SPACF = (CYL, (0,2) ) ,UNtT=(OISK,A) 

//« 

//• CALCOMP PLOT TAPE. 

//* 

//PLOTAPF 00 UNlT=2AOD-7,OCB»OFNsI,LAB6L*(.BLP>,OISP=(N6W,KEeP) , 

// DSN»0,PinT ,VOL=SER=35033 


00019<?00 
00020000 
00020100 
00020200 
00020300 
OOO20AOO 
0OO20S00 
00020600 
00020700 
00020800 
00020900 
00021000 
00021100 
00021200 
00021300 
00021400 
00021 SOO 
00021600 
00021700 
00021000 
00021900 
00022000 
00022100 
00022200 
00022300 
00022400 
00022500 


Figure 7-44. 
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The following notes, which apply to the GO-step DD cards, describe alterna- 
tives available to the user. In all of the examples it is assumed that the 
FORTRAN unit numbers are as shown in the sample NAMELIST (default values 
for all unit numbers except IPRINT = 6 (BIASNL), ERUNIT = 6 (BIASNL), and 
IDISK = 51 (MAIN). See Section 7. 3. 8. ) 

9 The SYSPRINT data set may be directed to a printer if the user 
desires to obtain an MSAD SNAP dump if the system abends when 
operating in a graphic mode 

The following DD cards are required 

//SYSPRINT DO SYS0UT=A,0Ce=(RECFH=V8A,LRECL=137, 8LKSUE=7265, BUFNO=l ) , 000017A0 

// UN1T=(CI SK,4) .SPACE- (CYL ,(C,2) J 00001745 

The user should note that if these cards are provided, it may take 
several minutes of real time for thj dump to be generated, during 
which time the program will not respond to the operator. If the 
program abends in a nongraphic mode, if the GESS Executive is 
unable to intercept the abend in a graphic mode, or if the operator 
does not attempt to recover or is unable to recover -from the abend 
in a graphic mode, then a dump is always provided on the 
SYSUDUMP data set. 

• If a spacecraft attitude tape is to be used to simulate data for a 
time-varying attitude then the following DD cards are required: 

//FT12FC01 CD UNIT=24CO,DISP=(OLO,KEEPJ.OSN=2B2HSATT,LA8EL=(,ELP), 00001880 

// DCB=(BECFP*V8S.LRECL«50,eLKSIZE=3004,DEN=3,8UFNO=ll. VOL-SER-XXXXX 00001885 

where XXXXX is replaced by the library number of the attitude tape. 
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• If it is desired to obtain solar and lunar epliemeris data in equinox 
of 1950.0 rather than equinox of date, using the SUNRD routine, then 
the following DD card is required: 

//FTl<iFC01 CD UNIT=DI SK,0 1SP»SHR, V0L=SER=00DS01. DSN=CT0S.SLPIS50. JAN71 C0001900 

• If a 2260 graphics device is to be used, then the following DD card 
is required: 

//FT23FC01 CD UNIT=2260 00001905 

• If the run is to be made in a nongraphics mode, then the DD card 
with sequence number 1905 must be removed. 

• If it is desired to obtain limar and solar ephemeris data from a 
JPL Ephemeris File, then the following DD card is required: 

//FT20FOOI CO UN1T'=2ACO,DISP=(OLD,KEEP),OSN=ZB2M'‘.JPL.LABEL = < .ELP)» 00001910 

// DCB*<ReCFF=V8S tLRECL=e276 ,6LKS I ZE= 8280 , CeN=3 , BUFN0= 1 ) r VCL=SER= XXXXX 00001915 

where XXXXX is the library number of the JPL Ephemeris tape. 
(Note that the link step must be changed to include subroutine RJPLT 
if the user desires to access the JPL tape. See the notes on tlie linlc 
step above. ) 

• If it is desired to obtain spacecraft ephemeris data from a GTDS 
sequential file then the following DD cards are required: 


//FT29F001 CD UN I T=2900 ,0 ISP= (OLD , KEEP J , DSN=ZB2MSG TO , LABEL- ( , BLP 1 t 00001920 

// 0CB=(RECFH=V8S,LRECL= 1096 ,BLKS I ZE= UOO, 0EN=3 , BUFN0= I ) , VOL =SER= XXXXX 00001925 


where XXXXX is the library number of the GTDS ephemeris tape. 
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e If it is desired to obtain the spacecraft ephemeris from a DODS 
EPHEM tape, then the following DD cards are required; 

//FT30F001 CD UN I T = 2A 00 , D 1 S P= (OLD , KEE P > . DSN=ZB2MSD0D, L ABEL= ( , BLP ) , 00001930 

// DCB=(RECFM = VBS,LRECL=28C'i ,B LKS I ZE = 2808 , D£N = 3 , BUFND= 1 ),VCL*SER = XXXXX 000019 3 5 

where XXXXX is die library number of the DODS EPHEM tape. 

• If it is desired to obtain the spacecraft ephemeris from a GTDS 
direct-access file, then the following DD cards are required: 

//FT31F001 00 UNIT=0l SK,0ISP=SHR,V0L=SER=XXXXXX,DSN=YYYyYY 00001940 

where XXXKXX is the name of the disk pack on which the GTDS 
file resides, and YYYYYY is the data set name. 

® If it is desired to read an AE data set, then the following DD card 
is required: 

//FTA9F001 00 UNIT=DISK,0ISP-SHR,V0U=SER=XXXXXX,0SN=YYYYYY 00001950 

where XXXXXX is the name of the disk pack on which the AE data 
set resides, and YYYYYY is the data set name. 

® If it is desired to process real data from an OABIAS data set then 
the following DD card is required; 

//FT50FC01 CD UN I T = 0 I SK ,0 ISP= SHR , V0L = SER=XXXXXX, 0SN= YY YYYY, DC B=BUFN0= I 00001960 

where XXXXXX is the name of the disk pack on which the OABIAS 
data set resides, and YYYYYY is the data set name. 

NOTE: To read this data set, the parameter IDISK (NAMELIST MAIN) must 
be restored to its default value of 50. In order to avoid operator con- 
confusion, the convention has been established that imit 51 is used for 
a scratch OABIAS data set for simulating data or copying AE data; 
unit 50 is used for a real OABIAS data set, which should not be written 
over by the MSAD/OABIAS System. 
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If printed output of the summary reports from the simulator are de- 
sired, then the following DD cards are required: 


//FT94FC01 CD SYSCUT = A, CCP= ( RECFH = VBA,LRECL = U7, BUKS IZE= 1 37A , BUFNO=l ) . 00002040 

// UNIT=(DISK,4) ,SPACE=(CYL i(0,2) I 00002045 


//FT95FC01 CO SYSCUT=A,0CE*(RECFM=veA,LReCL=l37,BLKSUE=1374,BUFN0*l), 00002050 
// UNIT=(01SK,4),SPACE*(CYL,(0,2n 00002055 


(These emits are normally dummied to save core storage). 

• If it is not desired to generate CalComp plots of plot displays, then 
the DD cards with sequence numbers 2090 and 2100 must be deleted. 

Figure 7-45 shows the JCL required to execute the MSAD/OABIAS System as- 
suming that the executable load module resides on ATTIT. OPRLIB. LOAD and 
the GESS nonresident tables reside on ATTIT. AODA. OABIAS.NRT, DATA. 

The notes given above for the GO-step DD cards apply equally to this example, 
ample. 
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I 


//ZB2f'S00'i JOB IGH7001857A,T,G00^3l, 004010), FFF,MSGLEVEL=t2,0) 

//♦ 

//♦ EXECUTE FSAD/QABIAS FROM A LOAD HGOULE AND MSAD TABLES ON GISYS4. 
//♦ 

//GO EXEC PGM=0AB1 AS,REGICN=4C0K 

//STEPLie CC UNI T=CISK,CI <P=SFR.V0L=SER=G1SYS4, 

// DSN=Gl-ZBPEP.OABI AS.LCACHCD 
//G0.FT05FCC1 DO ♦ NAPELIST INPUT- 

EFAIN ISIR=l ,0MEGAI=10C. ,NFR /1ME = 25 , NFP AM2=25 , 

ICISK=51, 

INC1SE=2,ST0V=4*0.005, SCUANT=0.5, 

I EPHEM = 0, ITERH=0, 1CB = 0, EEND 
£LIST IATAPE=l,ICB=l,ICaLAT=C, 

ISUN=1,ISPC=1,IPCCN=0, 

ALPHA=3 5 1. ,0ELTA=-20. , 

CREITE=l974.,l,,4.,22,,51,.2£., 

A =24 54 8. ,E=0.73264,EYE=28. 3,£MO=0.,WO=180.,RANODE=260., 

IFLAG=2, ISKIP=1,TH£TAC=0., 

CANGE = -9C.,STEP = 0,,S1CPA = 86. , 

EPRBET=0.25,ERRGAP=0.1,EPRA=C. 1 , ERRAD=0. 5. ERRT I M= 30. , 

T0=l<574. , l.,4.,2 3. ,04. ,DELTAT = 15., CENO 

ECPPANl CUFLAG=l,EPSlLN=C. ,D£BLG=8, 1 SNPRO-l , 1CBLAT=0, 

DANGE=-9C, , 

ISUN=1,ISPC=1,IPCCK=0, 

TCRBIT=1974. , 1 - , 4. , 22. , 5 1 . , 2 6. , 

A=24548. ,E=0. 73264, EYE=2 8.3, EM0=0. ,W0= 180. ,RAN0DE=260., 
EPRBET=0.25,ERRGAP=0. 1 ,E PR A=C. 1 , ERRAC=0.5 , ERRT I M=30 . , 

I AP1CR=8,APRA=350.,APDEC=-19., CENO 

CelASNL ATT0PT=5,REJCPT= 1, BE SMCD=6*T , F , T , RCOEL= 6* T , F , T , 

IPRINT = 6,ERUMT = 6, IUUNPL=4, I LEVEL = 12 , 1 PLOI= 12*3 , REDUNT= 10*0 , 
CAGRAN=.5,STPGRN=.l ,SANGRN=3C. ,RH0GRN=2., 

XC=12+C., P0=3*1CC.,9*10., CENU 
//♦ 

//♦ SYSOUT CATA SETS. ASSURING 1 PR I NT = ERUN 1 T = 6. 

//♦ 

//* PRINTOUT FROM RSAO, OADRIV, ODAP, CASYS, ANO OABIAS. 

//GO.FTC6FC01 DO DCB= I REC FH= V EA , LRECL= 137 , BLKS I ZE=7265 , BUFNO* 1 ) , 

// SPACE=(CYL ,( 2,2) ) ,UNIT = (OI SK,4) ,SYSCUT=A 
//♦ PRINTCUT FROM OABIAS. (IFRINT+1). 

//FT07FC01 CD SY SCUT= A , CC 8= ( R EC FH=VBA, LRECL= 1 37, BLKS I ZE= 1 374 , 8UFN0= 1 ) , 
// UN1T = ID1 SK,4) ,SPACE= fCYL ,(0,2)) 

//SYSPRINT CO DUMMY MSAC SNAP DUMP DATA SET. 

//P PRINTCUT FROM OABIAS. (IFRINT*2). 


OOOOOOlO 

00001360 

00001370 

00001380 

00001385 

00001386 

00001367 

00001390 


C0001650 

00001660 

00001670 

00001680 

00001690 

00001700 

00001710 

00001720 

00001730 

00001740 

00001750 


Figure 7-45. JCL to Execute IVISAD/OABIAS (1 of 2) 
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I 

M 

to 


//FT08FC01 DD SY SOUT = A , DC6= ( R ECFM = VB A , LRECL= 137, BLKS I ZE= 1374, 8UFNa= 1 1 , 
// LNIT=(CISK,4 3,SPACE=<CVL ,(C,2)) 

//♦ PRINTCLT FROM CDAP (EARTH AS CENTRAL BODY), AND OABIAS (IPRlNr*3). 
//FT09FC01 CO SYSCUT=A,DCe=(RECFH=VBA,LRECL=137,BLKSIZE=1374,BUFN0=l), 
// UNIT=(DI SK,4) ,SPACE=(CYL,(C,2) ) 

//♦ PRINTCLT FROM OOAP (MCON AS CENTRAL BODY), AND OABIAS (IFRINT+4). 
//FTIOFCOI CO SYSCUT=A,DCE=(RECFH=VBA,LRECL=137,BLKSIZE=1374,BUFN0*1 ), 
// L'NIT=(OI SK,4 ) ,SPACE= (CYL ,( C,2) ) 

//♦ PRINTCLT FROM OABIAS- (IFRINT+5). 

//FT 1 IF CO 1 CO SYSOUT = A,DCB= (RECFH*VBA,LRECL=137,BLKSIZE=1374,EUFN0=1 ), 
// UNIT=(OISK,4),SPACE=(CYL,(C,2) ) 

//♦ 

//FT12F001 CD DUMMY SPACECRAFT ATTITUDE TAPE. 

//♦ SUNRD LUNAR AND SOLAR EPHEMERIS FILE. 

//FTI4FC01 CD UNIT=OISK,DISP=SHR,VCL=SER=DODSOl,DSN=GTDS.SLPTCO. JAN71 
//FT23FCOI CD UNIT=2250 GRAPHICS DEVICE (IFTUBE). 

//FT28FC01 CO DUMMY JFL LLNAR AND SOLAR EPHEMERIS FILE (NRJPLT). 

//FT29FC01 DD DUMMY GTOS SEQUENTIAL SPC. EPHEMERIS FILE (NGTCS). 
//FT30FC01 CD DUMMY DCDS SPACECRAFT EPHEMERIS FILE (NCRBl). 
//FT31FC01 CD DUMMY GTCS DIRECT ACCESS SPC. EPHEMERIS FILE 
//FT49FCOI DO DUMMY AE DATA SET. 

//FT50FC01 CD DUMMY OABIAS DATA SET. 

1 1 * 

//♦ SCRATCH DATA SET FOR SIMULATED OABIAS DATA. 

/t* 

//FT51F001 CD UNIT = (OISK,4) ,SPACE=ICYL,(2,2n , 

// DCB=(RECFM=FB,LRECL=5C,BLKSIZE=2000,BUFNO=1) 

//FT94FC01 CD DUMMY OCAP REPORT PRIivIOUT (EARTH AS CENTRAL BOOY). 

//FT95FCOI CD DUMMY OCAP REPORT PRINTOUT (MOON AS CENTRAL BODY). 

//FT96FC0I CD UN I T=D I SK ,D I SP= SHR , VCL=S ER=G1 S YS4, 

// OSN=Gl.ZBMEP. OABIAS. TABLES 
//FT97FC01 CD « MSAD NAMELIST 
GCONTRL CENC 
//• 

//♦ CALCCMP PLOT TAPE. 

//♦ 

//GO. PLOT APE CD UN I T = 7TR ACK , OCB = DEN= 1 , L ABEL= ( , BLP ) ,OSN=ZB2MSPLT, 

// DISP=(CLD,KEEP) ,VOL=SER*3l636H 

//GC.SYSUCUMP DD S YSOUT = A ,•$? ACE= ( C YL , ( C , 2 ) ) , UN IT = ( 0 1 SK , 4 ) 


00001760 

00001770 

00001780 

00001790 

00001800 

00001810 

00001820 

00001830 

C0001840 

00001850 

00001860 

00001870 

00001880 

00001890 

00001900 

00C01905 

00001910 

00001920 

00001930 

00001940 

00001950 

00001960 

00001970 

00001980 

00001990 

00002000 

00002010 

00002040 

00002050 

00002056 

00002057 

00002058 

00002060 

00002070 

00002080 

00002090 

00002100 

00002110 


Figure 7-45. JCL to Execute MSAD/OABIAS (2 of 2) 
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